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The aim of the work presented in this thesis is to produce f., 
direct method to design structures subject to deflection constr¬ 
aints at the working loads* The work carried out can be divided 
into four main parts. In the first part, a direct design procedure 
for plane steel frames subjected to sway limitations is proposed. 

The stiffness equations are modified so that the sway in each 
storey is equal to some specified values. The modified equations 
are then solved by iteration to calculate the cross-sectional pro¬ 
perties of the columns as well as the other joint displacements, -v 
The beam sections are selected initially and then altered in an 
effort to reduce the total material cost of the frame, A linear 
extrapolation technique is used to reduce this cost. In this* 
design, stability functions are used so that the effect of axial 
loads in the members are taken into consideration. The final 
reduced cost design is checked for strength requirements and the 
members are altered accordingly, * .. 

In the second part, the design method is,applied to the - 
design of reinforced concrete frames in which the sway in the 
columns play an active part in the design criteria. The second 
moment of area of each column is obtained by solving the modified 
stiffness equations and then used to calculate the minimum column ' 
depth required. Again the frame has to be checked for all the 
ultimate limit state load cases. 

In the third part, the method is generalised to design pih- 
jointed space frames for deflection limitations. In these the 
member areas are calculated so that the deflection at a specifi§4,*i 
joint is equal to its specified value. 

In the final part, the Lagrange multiplier technique is 
employed to obtain an optimum design for plane rigidly jointed 
steel frames. The iteration technique is used here to solve the 
modified stiffness equations as well as derivative equations 
obtained in accordance to the requirements of the optimisation 
method. 
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CHAPTER ONE. %■ ; = 

HISTORICAL REVIEW AND SCOPE OF THE PRESENT WORK 

1=1 INTRODUCTION 

In multistorsy buildings the wind load is usually resisted 
either by a stiff core or by specially designed bracing frames= 

In such structures the frames may be designed to resist vertical 
loads onlyo However, such construction may, in some cases, be 
found undesirable in order to satisfy other functional require¬ 
ments of the building. In these cases the transverse frames are 
designed as plane sway frames and the horizontal forces are taken- 
by rigid frame actiono 

There i^ a wide range of possibilities open to the designer by 
which the member properties of the frame can be obtained. This 
depends on the design criteria imposed. Fundamentally, the object 
of the design is to find the most economical structure consistent 
with the safety and serviceability requirements. Thus, the design 
of sway frames often includes the following;- 

1) Safety - based on the ultimate strength of the structure 
(strength requirements) 

2) Serviceability - the deflection of the structure should not 
adversely affect the appearance of the structure (deflection 
requirements) 

3) In addition the cost of the structure may be as small as 
possible. 

Methods of designing a sway frame may be classified into three 
categories. In the first category the strength requirements are 
considered to be of overriding importance. In the second type both 
the strength and deflection requirements are considered to be 










important; and in the third, th,e optimum design is aimed at and 
in which it is possible to satisfy all the design criteria simul-' 
taneousl/o Thus, in this chapter the design methods are reviewed 
according to the above classification. 

1.2 DESIGN METHODS FOR STRENGTH REQUIREMENTS 

1.2.1 Early Design Methods for Multistorey Frames 

The first detailed investigation into the design of multi¬ 
storey building frames was carried out by the Steel Structure 
Research Committee (SSRC) from 1929 to 1936. Before this period, 
the accepted method of design for frames braced against wind was 
that now known as 'simple design' (Allwood, et al, 1961). In this 
method the beam-column connections are considered to be pinned. 

Thus, the beams are assumed to be simply supported. The columns 
are designed for axial load, as well as for the end^^m^^ 
by the beam reactions acting at an assumed eccentricity with respect 
to the column axis. The design is based on linear elastic behaviour 
of the structural material. Sections are chosen so that the maxi¬ 
mum stresses in a member do not exceed certain permissible values. 

An extension of this method to the design of rigidly-jointed 
sway frames, in which there is no wind—bracing, is explained in the 
Steel Designers Manual by Allwood, et al (1972). The internal forces 
in the members due to vertical loading are calculated as above. 

Those due to wind loading are found by making the frame statically 
determinate. This is done by assuming points of contraflexure at 
both the mid-height of the columns, and the mid-span of the beams. 

stage of design connections are taken to be rigid. The 
total internal-member forces are then obtained by superposition. 

The Steel Designers Manual gives three simple design methods. 





all of which use the aboye assuipptiono The jiiost eonffidnl^ used- 
of the three methods is the portai method;*whieh**a 
each bay acts as a simple portal and^theCtotal horizontal Idad'tis 
divided between the bays in proportion to the spans of the bays 

The final report of the Steel Structures Research Committee 
and the work by Baker (1954) showed that practical joints made of 
rivets or high tensile black bolts do not act as being pinned under 
the effect of vertical loading. The rigidity of the joints pro¬ 
duces restraining moments at the ends of a beam reducing the 
maximum design moments by some 17% to 25%, and causing large end 
moments in the columns. They have also shown that all the columns 
in a frame do not necessarily bend in double-curvature under wind 
loading, particularly in frames with weak beams. This clearly 
means that the assumption of a hinge at the mid-height of a column 
is not always true. However, the safety of the above methods pro¬ 
bably arises from a judicious choice of numerical values for per¬ 
missible stresses, eccentricities, etc. 

It has long been realised that if the individual members of 
a frame are designed according to the requirement to satisfy the 
allowable stresses, then the multistorey frame has a greater load¬ 
carrying capacity than permitted by the elastic design method. 

This has led to the application of rigid-plastic theory in struc¬ 
tures, which enables the design of a frame to be based on its over¬ 
all ultimate strength, 

1,2,2 Rigid-Plastic Design Method 

This was first developed in Cambridge by Baker, Horne, and 
Heyman (1956), One assumption made in the analysis of structures 
by the plastic theory is that the strains and deformations in the 










structure, at the limit gf proportionality, are small compared to 
those at the onset of strain hardening, . The theory therefore 
neglects the elastic strains and deformations. It assumes that 
the structure remains undeformed until suddenly a number of fully 
plastic hinges develop at discrete sections in the structure which 
converts it into a collapsing mechanismo At the begining of 
collapse (the instant when the structure is just converted into a 
mechanism with the formation of the last hinge) the structure is 
isostatic and, the bending moments throughout the structure can be 
calculated by using the static equilibrium of forces alone. 

The assumption of negligible deflections in the pre-collapse- 
state is generally applicable in frames braced against sidesway 
and in very few cases of low-rise sway frames where the effect of 
wind load is very small. In general the horizontal deflections in 
sway frames are of major importance and can not be neglected. 

They cause rapid deterioration of the overall stiffness of a frame. 
In fact. Wood (1958) shows that it is possible for a frame to lose 
all its stiffness, and therefore collapse, before the occurrance 
of a full mechanism. This results in the frame collapsing at a 
load factor less than the plastic failure load. He has shown that 
in a four-storey one-bay frame, after the formation of only two 
hinges, the modified critical load of the remaining parts of the 
frame is less than half of the original load. Furthermore, the 
hinges formed did not correspond to those of the rigid-plastic 
mechanism either in their position or in the load levels at which 
these hinges formed. Thus, the overall effect of instability in 
multi-storey sway frames is significant. For this reason, in 
recent years several methods have been proposed to include the 
instability of the frame into the rigid-plastic design. 







lo2„3 Plastic Design Methods £ox Praxes Instability 

Heyman (1960), proposed an approach to theniltimate ioad^’ * 
design of multi-storey sway frames by assuming a system of plastic 
hinges which involved collapse in both beam and columns. Hinge 
moments of the beam and columns are calculated with the factored 
load. Heyman suggests the use of the beam hinge moments calcu¬ 
lated to design the beams according to full plastic moment values. 
For the columns, because of the possibility of instability effects, 
Heyman proposed that the columns should be designed to remain 
elastic under the combined action of the hinge moments and axial 
forces. As a safeguard against instability it was suggested that 
a frame should be elastic at the working load without allowing 
excessive sway deflections. A method was given for the estimation 
of these deflections which made no allowance for the reduced frame 
stiffness due to instability. 

Heyman's design method has the advantage of directness and 
simplicity. Frame instability is not considered, but there is a 
degree of conservatism in the design procedure introduced by the 
elastic design of the columns. No suggestion was made as to what 
action should be taken if the deflections were not satisfactory. 

Holmes and Gandhi (1965) proposed a design method in which an 
allowance was made for instability. This was restricted to regular 
rectangular frames in which collapse was assumed to occur by beam, 
combined and sway mechanisms in the upper, middle and lower regions 
respectively. An initial design, obtained by considering this 
pattern of hinges, was modified to allow for: 

i) axial load effects, considering the reduction of column stiff¬ 
nesses due to compressive axial forces by using the ’stability 








functions' m, n and 0;-. 

ii) the reduction of beam stiffnesses,, _ due. to, the formation of 
plastic hinges; and 

iii) points of contraflexure not occurring at the mid-height of 
columns. 

The design procedure was repeated until two successive designs 
were identical. 

More recently Holmes and Sinclair-Jones (1970) have used the 
same principle and developed a more accurate method for modifying 
the initial design, and also for dealing with boundary regions of 
the frame, such as the top most storey, the bottom storey and the 
external columns. 

All the above methods reviewed in this section made use of 
an assumed pattern of hinges, which might not correspond to that 
which would form in practice. In fact, Majid and Anderson (1968), 
and Anderson (1969) show'that the assumed pattern of hinges was 
unrealistic and far from the actual behaviour of the frame at 
collapse. Furthermore, there can be no certainty of either the 
safety or the economy of a design obtained by such methods. 

The above methods did not design for deflection although some 
of them pointed to the significance of excessive deflection in 
tall frames, 

l,2o3.1 Ultimate-Load Design of Multi-storey Sway Frames 

An empirical expression for the actual failure load of a frame 
is given by Merchant (1954) as. 

1/A. = 1/A^ + 1/^p (I.n 

where X and X are the elastic critical load and plastic failure 






load respectively. Approximate methods for calculating are. 
given by Wood C1974;), Horne (1975;), Bolton C1976.), Moy C1976a), 
Williams (1977) and Williams (1979), A rapid calculation of 
for use in equation (1,1) is given by Neal and Symonds (1952), 
The B/20 Draft Specification (1977) permits the use of 


equation (1,1) within certain limits. The 


se are: 


i) If the ratio of A /A ^ 10, then A. is limited to A 

c p f p 

ii) If 4, then the above equation cannot be used and 

a full non-linear elastic-plastic analysis is necessary. 

The reason is that for this category of frames equation (lol) will 
give higher values of A^ than those given by accurate elastic- 
plastic analysis. 

The first condition is stated above in general applicable to 
braced frames, A method for designing a braced frame for ultimate 
load is given by the Joint Committee’s second report (1971), In 
this method it was assumed that it is sufficiently accurate for 
design purposes to consider only that limited part of the frame to 
which the member is directly connected. The limited frames for a 
beam and a column design are shown in Figures 1,1 and 1,2, in 
which the neighbouring members are assumed to be fixed. However, 
this design method is approximate and it is only applicable to non¬ 
sway frames. 

For sway frames, when A^/A^ does not exceed the limits imposed 
by the above conditions, equation (1,1) can be used to calculate 
A^, This involves the calculation of A^, Wood (1974) uses a 
slstitute frame based on the work of Grinter (1937), Figure 1,3b 


shows the substitute 
shown in Figure loSa^ 


frame used by Wood for the original frame 
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THE LIMITED FRAME FOR BEAM DESIGN IN THE JOINT COMMITTEE' 
METHOD 



1 2 THE LIMIT FRAjME FOR COLUMN DESIGN IN THE JOINT 
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IGURE 1.3 THE EQUIVALENT GRINTER FR.AME IN WOOD'S DESIGN METHOD 








The stiffness of a beam in the substitute. fra®e at any 
level is given by = 3 where k^, is I/L for a beam in the 

real frame, I is the second moment of area, and L is the span of 
the real beam. The stiffness of the column of a Grinter frame 
at any storey level is the sum of the stiffness of the columns 
in the real frame. As the wind and vertical loads on the frame 
increase, the rotational stiffness of each joint in the substitute 
frame decreases gradually. Frame instability is said to have 
taken place when the rotational stiffness of one of the joints of 
the substitute frame becomes zero or negative. The load factor at 
this stage is the elastic critical load, 

Williams (1977) modified the substitute frame of figure 1,3b 
to obtain a close lower bound on which is supposed to produce 
a safe design by equation (1,1), More recently, Williams (1979) 
presents a unified, exact treatment of the substitute frames which 
was used by Wood to approximate unbraced multi-storey frames. The 
method is only exact for the substitute frame and can be used to 
calculate the exact sway deflection for this frame 
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FIGURE 1,4 FRAf'E 


SUBASSEMBLAGES IN MOY'S METHOD OF DESIGN 






Moy Cl976a) proposed a method for calculating when the 
frame is subjected to vertical loads only. In this miethod points 
l^fl®ction are considered to develop at thehnid'-lehg the 

members. With such an idealization, a storey is assumed to consist 
of a series of subassemblages composed of two columns and of one 
or two restraining beams (Figure 1 , 4 ) depending upon whether the 
subasserablage is an exterior or an interior one. The concept 
upon which the method was based is that the stability of a whole 
frame is a function of storey stiffnesses. As the gravity loads 
increase, the storey stiffnesses decrease gradually. Frame 
instability is said to have taken place when one storey stiffness- 
becomes zero or negative, Moy accounts for the sudden change in 
the storey stiffness, when a plastic hinge forms in a member. This 
was reflected by a sharp drop in the storey stiffness. However, 
his method ignores the wind loads on unbraced frames, which have 
a major effect on the frame instability. 

In all the above methods reviewed in this section the frame 
is idealized employing mathematical models to depict the character¬ 
istics and behaviour of the real structure. These models may not 
accurately represent the real structure. All the methods calculate 
X , but the design using equation (1,1) also involves the calcu¬ 
lation of rigid-plastic failure load X^, Finding the right mechanism 
for a large multistorey frame is a complex operation. Moreover, 
the design using this equation may still involve a large amount of 
iterative calculations to obtain a safe design, 

A solution obtained by using the Merchant-Rankine formula 
(equation 1,1) may not satisfy the specified limitations on the 
horizontal deflections. Even when the facility for calculating these 
deflections exists, as in the method by Williams (1979), no 







recommendations are made tofimproye the "design if the deflections 
are found to be excessive. Furthermore^ the final,* design obtained 
by any of these methods may not be the most economical one. 

These methods were intended to provide the engineer with a 
manual design method or an easy method to program on a desk or 
pocket calculator. However, the lack of accuracy in these methods 
and the rapid development of digital computers after 1956 made 
it possible to carry out more accurate analysis of a frame. 

All the above methods are in no way design methods. This is 
because they select the cross section of the members and spend a 
great deal of time to carry out the "Analysis" of the frame. 

Another weakness of these methods is that they do not design the 
frames to satisfy the deflection (or sway) requirements. Once an 
analysis shows that the deflection in a frame is excessive, the 
above methods do not suggest methods for reducing it. 

1,2,4 Computer Methods of Design 

The matrix methods of structural analysis were first used by 
Livesley (1956), who developed methods to analyse elastic frames 
with or without considering the instability effects due to axial 
forces in the members. This method was developed by, amongst 
others, Livesley (1959), and Jennings and Majid (1965), who prepared 
programs which can be used to trace the actual load deflection 
behaviour of a frame up to collapse, 

Livesley, Jennings and Majid, used the matrix displacement 
method of structural analysis in the above programs. In this 
method, the unknown joint displacements are obtained by solving the 
matrix equation 


L = K X 


(lo2) 


where is the external load vector, K.i.s the overall stiffness 
raatirx, and ^ is the vector of joint displacements» Member forces 
are then calculated using these joint displacements„ Livesley 
C1956j introduced the effect of axial loads into the above 
equations by using stability '(j) functions' to modify the member 
stiffnesses. Iterative calculation is necessary as the axial forces 
in the members, required to calculate the stability functions, are 
unknown initially. 

The elastic-plastic analysis program proposed by Jennings 
and Majid led Home and Majid (1966) to suggest a method for the 
design by computer, of sway frames under proportional loading, A 
satisfactory structure was obtained by repeated cycles of elastic- 
plastic analysis and redesign. 

The design criteria which were adopted permitted collapse to 
occur only when a minimum load factor was exceeded. To satisfy 
these criteria, restrictions were placed upon the formation of 
plastic hinges in the beam and the columns. These restrictions 

were: 

1) No plastic hinge should develop in a beam below the load 
factor of unity, and the frame should be entirely elastic 

under the working load, 

2) Under combined vertical and wind loading, no plastic hinge 
should develop in a column below a specified permissible 

load factor. 

The method required that each proposed design is subjected to an 
accurate analysis. As a result, points of weakness and overstress 
are revealed, and can be removed. Tne final design satisfies the 
above design criteria, with efficient distribution of material 


throughout the frameo 








However, because of the large storage required for the’overall 
stiffness matrix in all the above methods, Jennings developed 

a compact method for the storage and solution of the stiffness 
equations o 

This technique was used by Majid and Anderson C1968) to 
develop an efficient method for linear and non-linear elastic-plasti 
design. The criteria used are those adopted by Home and Majid, 
stated above. 

The elastic-plastic design method is based on very realistic 
design criteria. However, the provision for the hinges not to form 
in columns below a permissible load factor, and in beams below the 
working load does not prevent excessive deflections in the frame. 

In fact, the frames designed by Majid and Anderson showed that the 
sway deflection of the frames were excessive, although the frames 
themselves had a reasonable reserve of strength. 

For example, for a fifteen storey frame the sway deflection 
ratio for the top storey at working load was 1/91, This is more 
than three times that required by the B/20 Draft, Furthermore, 
the elastic-plastic method is yet another method of repeated 
analysis in which the members are modified because an analysis 
reveals that such modifications are necessary, 

1,2,5 De sign of Reinforced Concrete Frames 

Until recently two methods have been in general use for the 
design of structural concrete members: the permissible stress method 
and the load factor method. These have sometimes been used as 

alternatives and sometimes used in combination. 

The permissible stress method was used in England prior to 
1957 for all reinforced concrete members. In the subsequent 



Code of Practice (^CP 114, 1957) it was retained as an altematiye 
to the load factor method for flejcure, with and without compression, 
and remained the basis of design for shear and bond. 

In the permissible stress method, the moments and forces 
acting on a structure were calculated from the actual values of 
the applied loads, but the limiting permissible stresses in the con¬ 
crete and the reinforcement were restricted to only a fraction of 
their true strengths, in order to provide an adequate safety factor. 
For example, in CP 114 (1957), the permissible tensile stress in 
steel reinforcement was normally the yield stress divided by the 
factor 1.8, while the permissible compressive stress for concrete 
in bending was the cube strength divided by the factor 3. Tlie 
method assumed a linear stress-strain relationship and a constant 
modular ratio of steel to concrete. 

The permissible stress method models the behaviour of a section 
under service loads fairly well, but it gives an unsatisfactory 
indication of condition as failure approaches. This is because the 
assumption of a linear relationship between stress and strain in 
the concrete no longer remains true, and thus the distribution of 
stress in the concrete differs from that under service load. 

The load-factor method of design was introduced into CP 114 to 
overcome this disadvantage. Theoretically, this method involves 
the analysis of sections at failure, the actual strength of a section 
being related to the actual load causing failure, with the latter 
being determined by ‘factoring’ the design load. The 1957 code for 
reinforced concrete (1967 reprint) permitted the use of this method 
with a central load factor of 1.8. However, in order to allow for 
the fact that there is a greater variation in the strength of 
concrete than in the strength of steel, it was specified that in 







calculating the strength, the cube strengths of concrete should first 
be reduced to the ratio of 1.8/3o This brings the load factor 
method into line with the alternative permissible stress methodo 


In the above two methods only failure has so far been con- 
sideredo A member was designed to have an ultimate resistance lar-; ; 
ger ■ than the applied load. However, it has been noticed that 
ser\'’iceability, no less than safety, is a requirement in structural 
design. This led to the concept of a limit state design used in 
the present Code of Practice CP 110:1972, 

Many text books have been published to explain the principles 
of this method (Reynolds and Steedman, 1974), (Allen, 1977), 

(Bennett, 1973), (MacGinley, 1978) (Astill and Martin, 1975), 

With this method, the design of each individual member or section 
of a member must satisfy two separate criteria: the ultimate limit- 
state, which ensures that the probability of failure is acceptably 
low, and the limit-state of serviceability, which ensures satis¬ 
factory behaviour under service (i,e, working) loads. The prin¬ 
ciple criteria relating to serviceability are the prevention of 
excessive deflection, excessive cracking and excessive vibration. 

Although CP 110 outlines serviceability limit-state calculations 
to ensure the avoidance of excessive cracking or deflection, a full 
analysis of every section was considered to be time consuming. 

For this reason the new code specified certain limits relating to 
bar spacing, span over depth ratio, slenderness, etc,, and if 
these criteria are not exceeded, more detailed calculation was con¬ 
sidered to be unnecessary. As it stands the code does not impose 
any limit on sway deflection of frames. The only limit is that the 
depth of the cross section of a concrete column should not be less 
than (HyjO), where is the effective storey height. 






Io2o6 Design of pin-Jointed Space Frames* 


Pin-jointed space frames are a good solution.to the problem of 
spanning large areas without intermediate supports» Their relative 
lightness combined with relatively high stiffness makes them attrac¬ 
tive as an alternative to beam or grillage systems. 

It appears that the design of a safe and serviceable space 
truss is considered to be a straightforward process, which requires 
a standard linear analysis follovv^ed by member sizings Thus, research 
into the design of pin-jointed space frames can be classified into: 
i) Research into the carrying capacity of the whole structure» 
ii) Research into the behaviour of structural elements« 

Little attention seems to have been given to investigation 
into the behaviour of the complete space frameo However, very 
recently [Schmidt, et al, 1980) carried out full scale tests on 
9'^o60 m space trusses o The tests showed that premature collapse of 
the complete truss was obtained, below that predicted by a conven¬ 
tional elastic analysis» This was due to two factors: 

a) The significant magnitude [about 9% yield values) of the 
initial member forces due to the fabrication of the joints 
and the members of the truss, and not attributable to the 
weight of the structure itselfo 

b) The reduced collapse load of a compression member below its 
average peak load obtained from the element testso 

Therefore, Schmidt’s tests showed that a design based on con¬ 
ventional elastic analysis, modified by a load factor, may possess 
a smaller margin of safety at collapse than its commonly supposedo 
Schmidt’s paper show that the mid-span^deflection of the truss was 
somewhat greater than the limit of span/360 specified by the B/20 
draft for maximum allowable deflection at working loado 
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Schmidt did not suggest the design o£ the space frame for 
deflection limitation. Their proposals , to modify the 'conven¬ 
tional elastic analysis to match the test results were to make two 
modifications to the analysis; 

i) A set oj initial lack of fit and assembly forces should 
be added to the externally applied load before the 
analysis is carried out, 

ii) The carrying capacity of compression members should be 
reduced by some 20% from tire theort ical value 
Although these modifications improved the analysis considerably, 
it is impossible to estimate the initial forces in the members 
without actual measurements. This is because, the initial forces 
depend on the fabrication and erection procedures. Thus, Schmidt’s 
analysis cannot be carried out before actually building the truss. 

The desigia of pin-jointed space-frame elements in tension or 
in compression is outlined in B/20 draft. To reduce the effect of 
buckling of the compression elements, smaller compression stresses 
are used. Design stresses are given in Table (6,2,2) of B/20 
draft for different values of the Robertson constant and slender¬ 
ness ratio, 

1,3 DESIGN TO DEFLECTION REQUIREMENTS 

The tendency for deflection criteria to dominate the design is 
necessarily increased by the introduction of more refined design 
methods, and the introduction of higher-strength steels. High- 
yield steel provides the frame with more strength than mild steel, 
and with less expense, Needham (1977) has demonstrated that steel 
Grade 50 is more economical to use than (mild steel) Grade 43, 
However, deflection becomes more critical when a higher grade steel 


is used. 







In the design of concrete elements, when the permissible stress 

design was being used with traditional conservative design methods, 

problems due to excessive deflection were practically unknown„ 

However, cracking in service has recently become more commono This 

is due to the use of the limit state design which, because of a more 

effective use of the concrete compression zone, have resulted in 

shallow members. More significant, is the introduction of high- 

strength materials with their correspondingly higher service stresses. 

For example, in British practice, there has been a progressive 

increase in permissible stresses with a consequent reduction in 

overall safety factors. For a concrete having a mean cube strength 
2 

of 40 N/mm at 28 days, the maximum permissible stresses according 
to the British codes were as follows (Clarke, 1974)I 


Date 

CP114:1957 
CP114:1965 
Draft CPI10:1969 
CP110:1972 


Maximum permissible stresses 
10 N/mm^ 

12,2 N/ram^ 

13,66 N/mm^ 

18 N/mm^ 

(This figure is for ultimate design stress) 


This increase in stress together with a reduction in flexural 
rigidity (El) due to cracking, with no corresponding increase in 
modulus of elasticity, results in higher deflections, since deflec¬ 
tion in concrete is directly proportional to stress and inversely 
proportional to flexural rigidity. 

Figure 1,5 curve (a) shows how the margin of safety (defined as 
ultimate stress-working stress \ j. ^ „ 

- - working ' st?S^ -^ reinforcement in tension 

has reduced successively over the years (Rowe, et al, 1965) and 
(Clarke, 1974). Curve (b) of Figure 1,5 represents the same 
information for concrete in bending. 
















From what has been said above one could conclude that the 
need for design for deflection requirement has become essential» 
However, up to now few methods exist for such a design, 

Stevens (1960) suggested that the real basis for the design 
of sway frames should be the prevention of unacceptable deformations 
under the working loads, A design to satisfy such a criterion 
was obtained by selecting a curvature pattern that would produce 
specified deformations, and was compatible with a bending moment 
distribution in equilibrium with the external loads. Sections 
were then selected by using moment-curvature charts. 

No direct allowance was, however, made for instability. A 
further defect of the method was that initial assumptions of 
relative flexural rigidities had to be made for statically indet¬ 
erminate frames. 

A second design method due to Stevens was based upon the 
collapse state. Maximum overload deformations were specified 
and used in formulating virtual work equations corresponding to 
collapse mechanisms in the deformation frame. Member sections were 
then selected, and the resulting design analysed by approximate 
methods. If it was found that the specified deformations were 
exceeded, then a further cycle of design and analysis was necess¬ 
ary. 

Although in this method the effect of the specified deflections 
was included in the virtual work equations, instability effects 
were ignored, and premature collapse could therefore occur. 

Moy (1974 and 1976b) proposed a method for sway frames in 
which an initial design with adequate strength but inadequate stiff¬ 
ness was corrected to satisfy permissible limits on horizontal 
deflections. The frame was divided into subassemblages on the basis 




of the portal method. These subassemblages were shown earlier in 
this chapter in Figure 1.4, for exterior and interior bays. 

Expressions for the storey stiffnesses were obtained in terms of 
the member second moments of area of the storey, and in terms of 
the permissible sway deflections at working load. Two assumptions 
were made: 

i) The vertical loads were assumed to have a negligible 
effect on horizontal deflection. 

ii) A point of contraflexure was taken to exist under hori¬ 
zontal loading at the mid-height of each column (except 
in the bottom storey), and at the mid-span of each beam. - 
Thus, the frame was made statically determinate above 
the bottom storey, and each storey was considered in 
isolation. 

The method is based on hand calculation but it requires a pre¬ 
liminary design and repeated calculation to arrive at a feasible 
solution. 

More recently Anderson and Islam (1979), and Islam (1978) 
suggested a method for the design of multi-storey frames to sway 
deflection limitations. The assumptions used are those adopted by 
Moy, stated above. Expressions relating the sway deflection over 
a storey height to the second moment of area of the corresponding 
column and surrounding beam were derived. These expressions were 
linked to a cost function to obtain an economical solution. A 
typical internal subassemblage for an intermediate storey in Anderson 
and Islam's method is shown in Figure 1.6. The complete frame is 
shown in Figure 1.7. The expression for the second moment of area 
of an internal column in this subassemblage was given as: 
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FIGURE 1.7 HIE COMPLETE FRAME IN ANDERSON AND ISLAM'S METHOD 
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where q^ ^ «..., q^ are the cost factors of the upper beams, 

^2 1^ ^2 2^ ®°5 ^2 m cost factors of the lower beams, and 


^3,1’ ^3,2’ ^3,ni+l 


are the cost factors of the columns» In 


equation (1„3), is the total horizontal shear in the columns 

of the storey under consideration, while and are those in 

the storey above and below respectively„ The quantities 5,^, £> 2 , 

oooo, a are the bay widths, where m is the number of bays, 
3m' 

LC=J!.i +£,„ + £_+ . 0.0 + £ ) is the total width of the frame, A 
is the horizontal deflection of the storey, h 2 is the height of 
the storey under consideration, while h^ and h^ are those for the 
storey above and below respectively, E is Young's modulus of 
elasticity, I „ refers to the second moment of area of the first 

o , Z 

left internal column. The expression for the second moment of area 
of the second left bay beam in the storey under consideration was 


given by: 


_ CPjhj + ^2^2^ ^2^2 ^5,2 
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(1,5) 


when the values of I, - and I, ^ were calculated from equation 

O y2 1 5> ^ 

(1,3) and (1.5) the second moment of area for the remaining columns 






and beams in the storey under consideration were calculated by 
proportion using the following two expressions respectively 
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Equation (1,3), Cl.5), Cl.6) and Cl.7) were repeated for each 
storeyo 

The values of the cost factors were given in table 2,1 and 
2,2 of Islam's thesis for universal column and beam sections. To 
obtain a quick conservative design, Anderson and Islam have shown 
that the values of the cost factors can be taken as being equal 
to one. In fact, their design of the 6 storey, 4-bay fixed base 
frame shows no change in weight when the real values of the cost 
factors were used, 

Anderson and Islam derived similar expressions to those given 
by equation Cl.3) and Cl.5) to deal with the boundary regions of 
the frame, such as the top-most storey, the bottom storey and the 
external columns. However, they suggested that a quick design could 
be obtained if the design of the top storey was considered as the 
storey below it, and the design of the bottom storey was taken as 
the one above it. 

In equation Cl.3) Anderson and Islam took the sway of each 

storey A as being equal to the maximum allowable deflection of 

h./300, where h. is the height of the storey under consideration, 

1 ’ 1 

However, such a ratio is difficult to be maintained for the ground 
floor columns, because these are connected to the foundations and 






therefore deflect less than the others» 

The methods of Moy and Anderson and Islam are approximate, as 
hinges were assumed in the mid-length of each member to avoid con¬ 
sidering the frame design in its entirety, and in order to separate 
It into storey slices. Furthermore, the vertical loads have a 
considerable effect on the sway especially in non-symmetrical 
frames. Thus, the effect cannot be ignored. 

1<,4 OPTIMUM DESIGN 

There are two types of optimisation: the classical type (which 
uses calculus), and the programming type. 

The first type is efficient only when it deals with equality 
constraints. One method of this type is the Lagrange multiplier 
technique (Hadley, 1970 edn). For instance, it is required to 
optimise the function 

f(x^,X 2 , .o,, x^) = f(x) (1„8) 

subjected to m equality constraints: 

g.(Xj,X2, xj = b., i=l, m (1.9) 

The procedure is first to convert each of the equations (1,9) to 
the form: 

Gf = A^[b^-g^(x)], i=l, ,,, , m (1,10) 

where the new variable is known as the Lagrange multipler. The 
constraints are than added to the objective function Z to form 

F(x, X) = f(x) + A. [b.-g.(x)] (1,11) 

Here F(x,X) is called Lagrangian function. The difference between 
F(x,X) and f(x) is that F(x,X) has more variables than f(x) arid the 
turning points of F(x,X) also include those of f(x). The turning 
point for the unconstrained function F(x,X) is found by setting to 


zero the partial derivatives of F(x,x) with respect to each of the 
m+n variables , j = l, ,,,, n and i=l, o,,, m. Thus 

3F _ 3f m 

3x. 3x. " i^l ^i 3^^ j = l, ..«,n (1„12) 

3 1 j 

8F 

- g^Cx) = 0, i=l, , m (lol3) 
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An advantage of Lagrange's method, from a practical point of 
view, is that the optimum solution can be obtained with relatively 
simple numerical techniques. However, the method is only useful 
when the constraints are equalities. 

The programming methods, are often more efficient than the 
classical type. This is because they can deal with inequalities 
and with relatively large probelms. For these reasons, these 
methods are often used to design rigidly jointed framed structures. 

The principles of some of these methods are given in some of the 
recently published text books (Majid, 1972), (Majid 1974), and 
(Gallagher and Zienkiewicz, 1973), 

Both the rigid-plastic and the elastic design method can 
be used to formulate the optimisation problems. The former method, 
however restricts the design constraints to cover strength require¬ 
ments only and cannot deal with deflection constraints. 

In the elastic design method, the design criteria used are 
that the stresses in the members and the deflections at the joint 
of a structure should not exceed certain permissible values. It 
then becomes necessary to express the deflections and the stresses 
in the structure in terms of the design variables. This can be 
carried out by employing either of two matrix methods. 




In the case of the matrix displacement method the design 
problem becomes one of finding the sectional properties of the 
member so that three constraints are satisfied. These are the 
stiffness equalities and the deflection and the stress inequalities. 

An alternative formulation uses the matrix force method to 
find the same member properties while satisfying the compatibility 
constraints which are equalities, as well as the deflection and 
the stress inequalities. 

Whichever method is employed the optimum design of a frame 
is a non-linear problem even if the instability effects due to 
the presence of high axial forces in the members are ignored. In 
this case, the linear programming techniques cannot be directly used 
to solve such problems. Several non-linear programming methods 
are available to overcome this difficulty but these may lead to 
complications and increase the computer time when designing large 
frames. Linear programming methods may, however, be used in the 
solution of non-linear problems by using an iterative procedure. 

Each iteration consists of a small linearised step of the problem. 
Different methods of solving a non-linear problem by such linearis¬ 
ation and iteration, e,g, the cutting plane method, the piecewise 
linearisation method etc, have been discussed in detail by Majid 
(1974), 

Saka (1975) presented a general computer program for the 
automatic optimum design of realistic, rigidly jointed multi-storey 
frames, A non-linear programming algorithm was proposed for the 
minimum weight design of structures in which the joint displacements 
and stress constraints were included as design variables. The 
design problem was formulated using the matrix displacements 
method. The problem was then linearised using the cutting plane 



method, and the optimum design was obtained using the simplex 
method„ 

Although, this method achieved a design in which all the 
design criteria were satisfied, it was not applicable to large 
frames, as it was essentially a procedure requiring extensive stor¬ 
age space and consumed a considerable amount of computer time.. These 
difficulties can be overcome, however, by making use of the backing 
store of a computer (Al-Pasha 1981-82)o 

1.5 SCOPE OF PRESENT WORK 

The design of multistorey sway frames with realistic wind 
loading can be governed by the limitations specified by the codes 
on the permissible horizontal sway of each storey. These are now 
specified in terms of the storey-heights at the specified unfactored 
loads (8/20 Draft Specification), It has also been recommended that 
at a load factor of unity the frame should remain elastic in con¬ 
sideration of both strength and serviceability of structure, 

(Majid and Anderson (1968),(t1oy, 1976) and(Majid, Stojanovski and 
Saka, 1980). 

In general a rational approach to designing sway frame for 
deflection limitations is not adopted in existing methods, practiced 
in industry. Furthermore, economy in design is sacrificed in 
order to obtain a fast manual procedure which is either restricted 
in its application or inaccurate in its approach. 

Many of the recent proposals for improved methods of design 
(Moy) and (Anderson and Islam 1979) have been based upon assumed 
rather than actual behaviour. For example, the assumptions of 
hinges at the mid-point of each member, and the neglect of the 
vertical loads on the frame etc are proposed to simplify the design 
methods. Accurate design methods for deflection limitation has not 





been generally available. 

The work presented in this thesis began with- an attempt to 
find an accurate direct design method to limiting horizontal 
deflections at the unfactored loads. Chapter 2 presents a method 
for the design of plane rigidly-jointed multistorey steel sway 
frames to withstand specified values of sway deflections. Using 
these values, the overall stiffness equations ]L = modified 

and then solved by iteration not to calculate deflections but to 
calculate the unknown cross-sectional properties of the columns. 

Thus the method proposed is that of design as opposed to that of 
an analysis. Simple linear extrapolation is then proposed to 
progressively reduce the cost of the material of the frame. 

Chapter 3 gives suitable criteria to cover both the strength 
and the deflection requirements in sway frames. Using stability 
functions, the axial loads effects are catered for and a refined, 
more accurate design is obtained. The procedure is then computer¬ 
ized, and an explanation of the computer programming is given in 
detail in Chapter 4, 

Chapter 5 contains the examples solved by this program. The 
effect of the axial force and a comparison of the results with 
other methods are pointed out in this chapter. 

In Chapter 6 the design procedure is applied to the design of 
rigidly-jointed multistorey reinforced concrete sway frames. 

The depth of each column and beam in the frame is calculated first, 
to withstand specified sway in the columns. The reinforcement of 
the columns and of the beams are then calculated so that the strength 

requirements are also satisfied, 

Tne method of Chapter 2 can be used to design pin-jointed 
space frames, for deflection limitation, by formulating the problem 




in terms of the unknown cross-sectional areas« This is done in 
Chapter 7 where a computer program for"this purpose is also 
presentedo 

In Chapter 2 an approximate extrapolation technique is used 
to change the beam sections in an attempt to reduce the material 
cost of the frameo In Chapter 8 the results of Chapter 2 are then 
used in conjunction with the method of Lagrange multipliers in 
order to obtain an optimum solution for the design problemo This 
new extended method is also iterativOo The computer program for the 
optimisation method is also given in this chaptero 





- CHAPTER TWO . , , 
DESIGN OF STEEL FRAMES 


2.1 INTRODUCTION 

The design of a rigidly jointed sway frames requires that the 
overall stiffness equations =^^must be satisfied while, in 
addition, the stresses and the sway in each member are within upper 
bounds specified by some limit state design code. In the case of 
sway frames, it has been found (Anderson, 1969 and Majid and Anderson, 
1968) that the deflection or sway requirements are more severe than 
the stress requirements. For this reason it is possible to satisfy 
the sway requirements first. If the sway in each column is equated 
to its upper bound, then the solution of _L = _KX becomes easier. In 
this chapter it is assumed that the sway in each member is known and 
the equations L - _K^ are solved to calculate the cross sectional 
properties of the columns. An iterative method is suggested to solve 
these equations. This chapter only deals with the design of two 
dimensional rectangular steel sway frames. 

For a general joint in a sway frame, expressions are obtained 
to calculate the deflections in terms of those at the neighbouring 
joints as.well as the second moment of area of the connecting members. 
The value of the horizontal sway at each storey-level is first 
specified, and the sway equations are used to evaluate the unknown 
sectional properties of the columns. The stiffness equations for 
each joint are solved by an iterative technique. This produces 
values for the other deflections in the joint as well as the second 
moment of area of the columns. The beam sections are selected by a 
preliminary method such as that suggested by (Islam, 1978, and 
Anderson and Islam 1979). These are then altered in an effort to 








reduce the total material cost of the frame. 

Simple linear extrapolation is used to find the cost of each 
storey. The beam which produces the least cost is used in the next 
round of an iteration process, which is repeated until no further 
reduction in material cost of the frame can be gained. As the 
Iteration process continues, the beam sections, the column sections 

as well as the ratio between the external and the internal columns 
all change. 

2 »2 CONTRIBUTIONS OF A MEMBER TO THE OVERALL STIFFNESS MATRIX 

The contributions of a member to the overall stiffness matrix 


[Kj of a frame is:- 

First end Second end 



where 







A = a£^ + bm^ 

P P 

B = (a-b) £ m 
P P 

C = -d m 

F = am + b£'^ 

P P 

T = d£p (2,3 

e = 4EI/L 
f = 2EI/L 

a = EA/L 

b = 12 EI/l3 

d = -6 EI/l2 

In equation (2,3), E is the Young's modulus of elasticity, I is 
the second moment of area for the member, L is the length of that 
member, £p and m^ are the direction cosines for the longitudinal P 
axis of the member, and A is the member area. 

For rectangular skeletal frames it is convenient to simplify 
the K matrix and construct separate matrices for beams and columns. 

It is assumed that the direction of the P axis for a beam is from 
left to right, while that for a column is vertically downwards. The 
sign convention adopted for forces and displacements is in accordance 
with the right hand screw rule, and is shown for a horizontal and 

a vertical member in Figure 2,1, 

For a beam, the stiffness equations are simplified to become: 







Column 



FIGURE 2.1: 


THE SIGN CONVENTION AND THE POSITIVE DIRECTION OF 
THE LOCAL AND GLOBAL AXES 
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where L = (H. V. M. 

— Ill 


(x. y. r. X. y. r.} 
1 1 1 J 3 3 


a column: 


, 0.0 Vj } is the load matrix, 
is the joint displacement vector„ 


and {X} = 
Similarly for 
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(2.5) 


If the axial stiffness (EA/L) of a member is neglected, equations 
(2,4) are simplified further by removing the first and the fourth 
columns and rows from the stiffness matrix. Equations (2,5) are 
also simplified by removing the second and the fifth columns and 
rows of K, In this manner, the contributions for a beam and a 
column, to the stiffness matrix become the same. This is provided 
that the sign convention and the direction of the axes of the members 
are kept as specified. Thus the contributions of a member to 






{L} = K {X} becomes: 
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( 2 , 6 ) 


Notice that suffixes 


i and j refer to the ends of the member. 
First end Second end 
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2,3 MEMBER AND JOINT SPECIFICATION 

In this thesis, only orthogonal frames with horizontal beams 
and vertical columns are considered. This simplifies the numbering 
of the joints and the members and, thus the data preparation can be 
reduced to a minimum. In Figure 2,2 a frame is shown with joints 
and members numbered. Joints in the beams are numbered first. This 
is done starting at the left end of the top storey. The numbers are 
consecutive and are from left to right. The columns joints are then 
numbered in the same way, A member is given the same number as the 
joint at its first end, which is specified by the tail of an arrow 
placed on the member. This arrow also specifies the positive dire¬ 
ction of the longitudinal p axis of the member. Thus on the second 
floor, left bay of the frame shown in Figure 2,2, the beam connecting 












the two columns carries a load along its span„ This is at joint 4o.j 
The joint to the left of this is numbered 13. In our numbering, 
system, the beam to the left of joint 4 is beam because joint 

13 is at the first end of this beam. The beam to the right of joint 
4 is beam because joint 4 is at its first end. The left side 
outer column at the second storey is because joint 13 is at its 
first end. All supports are numbered as zeroes. 

2.4 STIFFNESS EQUATIONS FOR A GENERAL JOINT IN THE FRAME 

A column joint is at the end of a column while a beam joint 
lies somewhere along a beam. Figure 2.3 shows a general beam which 
may be carrying a load at joint j and other loads at i or k. Alter¬ 
natively, i or k may be joints at the head of columns on either end 
of the beam. In our specification, joint i is to the left of j, 
and k is to its right. Tlie beam joining i to j is beam with its 

second moment of area I .. The beam to the right of j joining j to 

d1 

k is beam B^ with its second moment of area 1^^. For any column 
joint, there are six possible configurations. These are shown in 
Figure 2.4. Of these, type 6 is the general case where two columns 
and two beams meet. This type is also shown in Figure 2.5 where 
joint j is shown together with its neighbouring joints. The second 
moments of area for the members connected to j are also specified 
in this figure. All the other types 1 to 5 can be obtained from 
this configuration. For this reason joint type 6 is used to derive 
the relevant equations for the design of a column. 

2.4.1 Stiffness Equations for a Beam Joint 

The relationship between the external loads ^ and the vectori- 
ally equivalent joint displacements in terms of the system coor¬ 
dinates XY, is expressed by the equation 









CONFIGURATIONS OF THE MEMBERS AT A COLUMN JOINT 
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L = K X 
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This is used to evaluate the deflections X, thus 


X = K" 


C2o9) 


For a beam joint j, the vector X is {x., y., r.}„ Here x. is the 

- 1 3 3 3 

horizontal displacement; y^. is the vertical displacement; and r^ is 
the rotation of the jointo If the axial stiffness (EA/L) in a beam 
is neglected, the displacement vector becomes X= {y^, thus 

two equations are needed for defining y^ and r^ o These are obtained 
from the contribution of the beams and connected to joint j, 
(Figure 2o3)<, The stiffness coefficients of beam B^ are: 




d. 

3 




( 2 . 10 ) 




where is the length of the beam. Similarly for B^ the stiffness 
coefficients are; 


b. 

1 




d. 

1 


e. 

1 






( 2 . 11 ) 


f. 

1 




where L . is the length of B., 

1 

2,4,1,1 The Vertical Deflection 


The vertical deflection y^ for a joint in a beam is obtained 




from the stiffness equations as follows: 

i) The contributions of the beam B. to the vertical deflec- 

1 

tion of joint j, which is at its first end, are calcu¬ 
lated from the first row of equations (2„6), thus 


V. 

J 



d.r. - b.y, - d.r 
J J J k 3 


Rearrange: 


= (V. + d.r. 
3 11 


+ b .y + d.r, )/b . 
ik : k" : 


( 2 . 12 ) 


(2,12a) 


ii) The contributions of the beam to joint j, which is 
at its second end of the beam, are calculated from the 
third row of equations (2.6), thus 


0 = -b.y. + d.r. + b.y. + d.r. (2.13) 

11 11 ij ij 

The vertical load is considered to be acting at end j 
of Bjo This is why in equation (2.13), the left hand 
side is zero indicating that no load is acting at the 
second end of B^. Rearranging equation (2.13) for y. 
gives: 

y. = (+b.y. - d.r. - d.r.)/b (2.13a) 

^3 ^i"^! 11 X 3 i 

iii) .A.dding the contributions of B^ and B^, given by equations 
(2.12a) and (2.13a) results in the following equation for 


V. = (V. + d.r. + b.y, + d.r + 
'1^1 11 1 k 1 k 


b.y. 

1 1 


d.r. 


1 1 


d.r.)/(b. + b.) 
1 ] J 1 


(2ol4) 




2o4,l„2 The Rotation r. 

_1 


Similarly the rotation r^ for a beam joint is obtained as follows: 

i) The second row of equations C2o6) gives the contributions 

of the beam B. to the rotation of joint j as: 

J 


M. = -d.y. + e.r. + d.y, + f.r, 
■J J J J J 1 k J k 


(2,15) 


Thus r. = (M. + d.y. - d.y, - f.r, )/e. 

J J J J 1 k J k' J 


(2,15a) 


ii) The fourth row of equations (2,6) gives the contributions 


of beam B. as: 
1 


0 = d.y. + f.r. + d.y. + e.r. 
1-^1 11 1 J 1 j 


(2,16) 


r. = (+ d.y. - f.r. - d.y.)/e. 

J I'l 11 1 1 


(2.16a) 


It should be noted that is taken into account when 
considering the contributions of B^, As a result, in 
equation (2,16) and (2,16a), is not reconsidered, 

.ii) The rotation of joint j is obtained by adding equations 
(2,15a) and (2,16a), This gives: 


r. = (M. + d.y. - d.y, - f.r, 
1^1 11 1 k J k 


(2,17) 


+ d.y. - f.r. - d.y.)/(e. + e.) 
11 11 1 y ^ j 1 


2,4,2 Stiffness Equations for a Column Joint 


The displacements vector for a column joint, is similar to that 
)f a beam, except that the axial deformation of the column must be 


considered. Thus, three equations are needed to find , y^ and 

r.. These equations are obtained from the contributions of the beams 

B and B., and the columns C. and C which are connected to joint j, 
j 1 1 m 







see Figure 2,5<. The column C. which connects joint j to n is 

1 

defined by the joint number j at its first end<. Its second moment 

of area is I„.o Members B., B. and C are defined similarly by 
Cj i’ 3 m , 

the joint at their first end„ The second moments of area of these 

members are I„., . and „ 

Bi’ Bj Cm 

The member stiffness coefficients a, b, d, f and e defined in 

equations (2„3) are found for each member connected to joint jo For 

column C. these coefficients are b^., d„., e^. and f„., while for 

3 C3 ’ C3 ’ C3 C3 ’ 

column C they are b„ , d^ , e„ and f^ , These coefficients for 
m Cm Cm Cm Cm 

the beams B^ and B^ have been given by equations (2.10) and (2oll)o 

The axial stiffness for column C. is: 

J 




C2ol8) 


where A . is the area of the column while the axial stiffness for 

the column C is: 

m 


^Cm ^^m'^^Cm 


(2,19) 


2o4o2,l The Horizontal Deflection x^ 


The expression for x^ (the horizontal displacement of joint j) 
is obtained from the stiffness equations as follows: 

i) The contributions of the column C^ to the sway are taken 
from the first row of equations (2,7), which is 


H. = b„.x. - d„.r. - b X - d„.r 
3 C3 3 C3 3 C 3 n C 3 n 

Thus X. = (H. -I- d„.r. + b X + d„.r )/b„. 

3 J C 3 3 C 3 n C 3 n^ C 3 


(2,18) 

(2,18a) 


iil For column C , joint j is at its second end, so its con- 
^ m 

tributions to the horizontal deflection x. are obtained 

J 

from the third row of equations (2,7), which gives: 





(2,19) 


0 = -b X + r + b_ X. + d- r. 
Cm m Cm m Cm j Cm j 


and thus: 


X. = (+ b^ X - d„ r - d^ r.)/b^ 
j Cmm Cmm Cmj Cm 


(2,19a) 


iii) There are no contributions from the beams to the 

horizontal deflection x. at j because the axial stiff- 

J 

ness of these is ignored. 


iv) The total deflection x^. at j is thus obtained by adding 


equations (2,18a) and (2,19a), Hence 


X. = (H. + d^.r. + b^.x + d^.r 
J 1 Cj J Cj n Cj n 


+ b^ X - d„ r - d_ r.)/(b_. + b_ ) 
Cm ra Cm m Cm j Cj Cm'^ 


( 2 , 20 ) 


2,4,2,2 The Vertical Deflection y. 


The vertical displacement of joint j is obtained from the con¬ 
tributions of all the members connected to it as follows: 


i) The contributions of the beams and to y. are given 
in equations (2,12a) and (2,13a) for a beam joint. These 
remain the same for a column joint, 

ii) The contributions of a column C. to the vertical deflection 
y^ are obtained by considering the axial deformation of 
that column. This is given by the second row of equations 
(2,5), Thus 


0 = a_.y. - a„.y 
Cj^: Cj^n 


( 2 , 21 ) 


That is 




(2,21a) 


in') Joint i is at the second end of column C , Thus the 

contribution of this column to the vertical deflection 






at j is given by the fifth row of equations (2<,5), thus: 


0 = -a^y + a_ y. 
Cm m Cm 3 

or y. = ay /a 

j Cm m Cm 


( 2 , 22 ) 


( 2 , 22 a) 


iv) The total value of y^ is then obtained by adding equations 
(2„12a), (2„13a), (2.21a) and (2,22a) to give: 


y. = (V. + d.r. + b.y, + d.r, + b.y. - d.r. - d.r. 

J J 3 3 1 K J K 11 11 13 


a„.y + a„ y )/(b. + b. + a„. + a. ) 
C 3 n Cm m 3 1 C 3 Cm 


(2,23) 


Notice that the external load applied at a joint is only 
considered once with the contribution of any member 
connected to that joint. In the case of y^, the external 
load at j was considered when y^. for the beam was 

derived. 


2,4,2,3 The Rotation r. 


The rotation r. of a column joint j is obtained as follows; 


i) The contributions of the beams and B^ to r^ are 

already derived and are given by equations (2.15a) and 
(2,16a), 

ii) The contributions of columns C^ to r^ are obtained from 
the second row of equations (2,7), Joint j is at the 
first end of this column. Thus: 




C2.24) 


i,e, r. = (+d^.x. - d„.x - f„.r )/e . 

' 3 ^ C3 3 C 3 n C 3 n^' C 3 


(2.24a) 


iii) Joint j is at the second end of column C^,. Hence the 

fourth row of equations (2,7) gives the contribution of 






this member to r.= Thus: 

0 = -d X + f r + d„ X. + e_ r. 
Cm m Cm m Cm j Cm j 


(2,25) 


i.e. r. = (+ d„ X - f r - d^x.)/e^ 
J Cmm Cmm CmyCm 


(2.25a) 


iv) Adding equations (2„15a), (2ol6a), (2o24a) and (2„25a) 
together, gives the final rotation of joint j as; 


r. = (Mj + d.y. - d.y^ - f.r^ + d.y. - f.r^ - d.y. » 


d_.x. - d_.x - f„.r + d_ X - f_ r - d„ x./ 
Cj j Cj n Cj n Cm m Cm m Cm y 


(e. . a. * Cj,. . (2.26) 

Equations (2o20), (2o23') and C2o26) express the joint 
displacements for a member with configuration 6o Figure 
2c4 shows this configuration^, which is that of a general 
case^ To obtain similar equations for configurations 1 
to 5, terms contributed by the missing members should be 
omittedo Eogo for configuration 3^ terms contributed by 
column C to all three equations should be omittedo 

In the orthodox stiffness method of analysis, to cal¬ 
culate the sway of a complete storey, one adds the rows 
and the columns in the stiffness matrix corresponding to 
Xj for all the joints in the storey^ Thus equation (2o20), 
which calculates the sway of a column joint, can be extended 
to calculate the sway of a complete storey by adding up 
the terms in the equation to include all the joints in 
the storeyo The numerator SN^. of equation (2,20) is 






Figure (2,6) shows an intermediate storey of height h in a frame of 
non-uniform bay width<. The simple "portal method" of design assumes 
that the horizontal shear at any level is divided between the bays 





in proportion to their relative width CAllwood, et al, 1966)„ . 
(Anderson and Islam 1979), Hence, for the columns. 


I 

L 


4 ^ 



L- + 

I j 


o o o o o o o 



C2» 


where nb is the number of bays. The ratio l2'^^l referred to 
as the internal to external column ratio, and it is denoted by p. 
For the first column joint in the storey (that is the one on the 
immediate left) is equal to 1, Hence for the rest of the joints 
of the columns p is defined as: 


' Vh' 


= lyq, 


= h'h- 


- lyh 


ooo, P 1^1 

n (nb +1) 1 


( 2 . 


Using equations (2,32): 


p^ = 1,0, p^ = (L^ + p^ = (L^ + 


O Q O O O ^ 


% = ‘■nb/q 


( 2 . 


If the sway of the complete storey is specified, then equation 
(2,29) can be rearranged to calculate I^, in terms of the sway x^, 
and the ratios p^. The second moments of area for the remaining 
columns in the storey are then calculated from equation (2,33), 

It should be stressed that the ratios p for the second moment 
of area of internal and external column has to be specified before 
the above equations are used for the design of the columns. How¬ 
ever, it should be stressed that the pontal method described above 
is not a necessary part of the design method. The poatal method 







is only used to begin with to specify some initial patios ,yj^||ate^ 
as the design of the frame proceeds, the value of p changes con¬ 
tinuously so that the final design has the most economic value of 
p for each storey. 


2,4o3„l Equation for Second Moment of Area of the External 
Left Column 

i) In the design procedure the contributions of each column 
in the storey to should be found first. These are 
calculated from equation (2,31), which is rewritten in 


terms of p. and I,. Thus, the contributions of columns 
1 1 

C. to I, are: 

1 1 


I-, = (H. - b^x. + b„ X - d„ r - d_ r.)/ 

1 ^1 Cmj Cmm Cram Cmj 

(2.35) 

(p. (b^.,x. - d^.,r. - b„. X - d_.,r )] 

^ 1 Cjl j Cjl j Cjl n Cjl n-'^ 

ii) Terms of equation (2.35) are added for all the columns 

in the storey. The numerator of equation (2.35) is SN! 

given by: 

SN! = (H. - b„ X. + b- X -dr - d r.) (2,36) 

j '■ j Cm j Cm m Cm Cm Cm j 


while the denominator is SDj 

= Oj ('>cnN - "cji'j - - "cnV 


(2.37) 


is thus given by: 


j~'is 

I, = SNI/.Z. SD! 

1 3 l=Jl 3 


(2.38) 


IVhere J, and J are as given before in equation (2.29). 
1 s 






2o 4<,4 Aii Example oii the Derivation of the Stiffness EqijatiQns 


The above equations are best explained by deriving them for 
a single-storey single-bay frame, shown in Figure 2o7. To keep the 
problem as simple as possible the axial deflection of the columns 
is disregarded at this stage. The stiffness equation L = K X for 
the frame in Figure 2.7 are; , —• 
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where r^, r and r are the rotations of joints B, F and C respec- 

B F C 

tively, X is the horizontal sway, and y is the vertical sway of 

r 

joint F. The second moment of area for the column and beam are 
and 

For the beam joint F, equation (2,14), or the third row of 
equation (2,39), gives: 


(L/2)2 
12EI, 


(L/2)2 


(L/2) 


(L/2) 


(2,40) 











while equation (2,17) or the fourth row of equation (2,39) gives 


4El2 4Eq 

—“—rrz 


(2,41) 


For the column joint B, the sway x can be obtained from equation 
(2.29), or from the first row of equation (2,39); thus 


6EI^ 6E1^ 

- r + - r 

h^ ^ h^ ^ 

12EI^ 12EI^ 


(2,42) 


If, X is specified (e.g, x = h/300) can be obtained from equation 
(2,42), thus 


1 ,12E 121 6E “SI 

(- + -3 X-r-r 

h3 h3 h2 ® h^ ^ 


(2,43) 


Application of equation (2,38) gives the same result. 

Tlie rotation r^ is found from equation (2,26), which is the 

D 

same as the second row of the stiffness matrix. Both give: 
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4EI^ 4EI2 

- + - 
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(2.44) 


similarly the rotation r^ is given by: 

6 EI 6EI 2] 

0 + —^ X— ^ y. - - 
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(2,45) 
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2,5 PRELIMINARY BEAM SECTIONS 

In the design process the beam sections should be selected to 
satisfy the stress constraints. First of all, a lower bound is 
obtained for a beam by considering the vertical load on it and 
deriving the equation of beam mechanism. In most cases, the section 
obtained in this manner is not adequate for resisting the side-sway 
except for the top one or two storeys. Nevertheless, no beam section 
is allowed to be reduced belovNr that obtained by beam mechanism. 

Under wind-loads, preliminary beam sections could be obtained to 
resist the sway deflection by a method due to Anderson and Islam 
1979) as described in Chapter 1, In some cases, sections obtained 
in this way are found to be totally unsatisfactory. However, 
Anderson and Islam's approach can still be utilized provided that, 
in this preliminary calculation, the allowable sway is severely 
restricted for the time being. It has been found that equations 
given by Anderson and Islam underestimate the beam sections con¬ 
siderably. 

Another method, which could be used to select an initial 
section for a beam, is to use the simple design methods of The 
Steel Designers' Manual (Allwood, et al, 1972), In these methods 
a point of contraflexure is taken to exist at the midheight of each 
column and at the raidlength of each beam. In most cases, sections 
selected by this method to satisfy strength requirements are found 
to be unsatisfactory for resisting the sway deflection. However, 
for a preliminary trial, these sections can be magnified by a cer¬ 
tain factor and obtain an upper bound for each beam. This factor 
might change from the beams at the top storey to those at the bottom 
and it can have any valuo between 1,5 and 2 at the top and 3 to 5 
at the bottom, depending on the height of the frame. 






2.6 SOLUTION OF THE MODIFIED STIFFNESS EQUATIONS BY 

INERATIVE TECHNIQUE 

The modified stiffness equations given for a general joint 
section 2.4.2 and 2.4.3) are applied to each joint of the 

frame to obtain a sufficient number of equations for the evaluation 
of the joint displacement and the column second moments of area. 

This is provided that sections have already been selected for the 
beams. 

Although these equations are non-linear, an iterative tech¬ 
nique, which is a significant part of this thesis, proved to be a 
powerful method for solving them. This technique starts by assuming 
estimated values for the unknoivns which, in the absence of a better 
value, may be taken as zero or one. A second revised estimate 
may then be found using the modified stiffness equations one at a 
time. The current value of a variable, at each stage of the iter¬ 
ation is used as soon as it is available to calculate a new value 
for the next variable in the iteration chain. 

It must be emphasised here, that only an unknown related to 
the leading diagonal element of the stiffness matrix is calculated 
from each equation. The obvious reason is that the coefficients 
of the leading diagonal element is much larger than the others 
(Lennox, et al, 1974, Cohen, et al, 1973). This prevents divergence 
and also accelerates the convergence of the process. 

2.7 AN EXAMPLE ON TriE SOLUTION OF THE MODIFIED STIFFNESS EQUATIONS 

The solution method is best explained by applying it to the 
single bay frame dealt with in section 2.4.4. Assuming r^ = r^ = 
r = 0, y = -1, x/h = 1/360, and is calculated from the beam 
mechanism, which gives = 21345 x 10*+ ram^^. Equation (2.40) gives 






Xp = -11092 jnnio For the rest of this cxcle, this new value of 
Xp is used for the calculation of all the other variables. Equation 
(2.41) gives: Vp = 0,0, and equation (2,43) gives the first value 
for the column section as = 3106 x lO^ With these new values 

for Xp» Ip, equation (2,44) gives r^ = -0,00395„ Thus r^ 

is calculated from equation (2.45) as 0,00267, This procedure is 
repeated until all the variables converge to stable results. To 
test this convergence, a tolerance is introduced, and a variable 
is said to have become stable if its: 


new value - old value 
new value 


< The tolerance 


Table 2,1 shows the results of the iteration process, with a tol¬ 
erance of 0,01, It can be noticed that 13 iterations are needed 
to converge to the final design values, which include the second 
moments of area. 


2,8 AN OUTLINE OF DESIGN PROCEDURE 


At this stage it is useful to give an outline of the design 
procedure. This consists of the following steps. 

Step 1 - Select an initial set of beam sections bx the beam 

mechanism method under vertical loads. 

Step 2 - Specifx the horizontal deflection at each storex level 

Step 3 - Define the ratio of the second moments of area of an 

internal and an external column, Anx ratio is acceptable 
but it was decided to start with the ratios calculated 
above in section 2,4,3, 

Step 4 - Assume the second moment of area of each column, the 

deflection and the rotation of all the joints. Infea¬ 
sible values such as I = 1 mm^*, x = 1 and r = 0,0 are 





sufficient for the purpose of starting the iteration • 
processo 

Step 5 - Using the modified stiffness equations given above, 

calculate the deflections of the beam jointso 

Step 6 - Using these values of deflections calculate the first 

unknown sectional property of one of the colunuio 
Using this value of calculate the vertical deflection 
and the rotation of that column joint» 

Step 7 - Continue with iteration until a new value is obtained for 

each unknown and all the modified stiffness equations 
for all the joints are use<l 

Step 8 - Repeat steps 5, 6 and 7 until convergence is achieved. 

At this stage one design cycle is completed and a set 
of sections become'available! for columns. All the joint 
deflections and rotations also become available. Notice 
that the beam sections and the horizontal deflection of 
the joints are specified but the other unknowns are cal¬ 
culated, The column sections and the joint deflections 
calculated fully satisfy the stiffness equations and 
the design obtained is feasible and most structures 
designed in this manner are ready to be constructed. 
However, the design process continues for the following 
reasons: 

a - To include the effect of axial forces in the columns. 

b - To achieve economy by altering the ratios between 
internal and external columns, 

c - To achieve economy by repeating the operation with 






a different set of beam sections. This is provided‘no 
beam section is reduced below those given by the beam 
mechanism, 

d - To check that while deflection and sway requirements 
are satisfied, the stress (plastic hinge moments) are 
satisfactory in all the members. 

Step 9 - Modify the stiffness equations by including stability 

functions, see later. 

Step 10 - Alter some of the beam sections if this is economical, 
see later. 

Step 11 - Repeat the iteration process, i,e, steps, 5, 6 and 7 

until convergence is obtained once again. Notice that 
when entering the iteration this time, the current 
values of column sections and joint rotations are emp¬ 
loyed. These values are more realistic than the initial 
infeasible (I = 1, y = -1 and r = 0) values. For this 
reason the second round of iteration converges very 
quickly. 

Step 12 - Alter the ratio between the internal and external col¬ 
umns, if this is economical, see later. 

Step 13 - Repeat the iteration, i,e, steps 5, 6 and 7 until con¬ 
vergence is achieved. Once again the iteration is 
initiated with current available values and convergence 
is very fast. 

Step 14 - Repeat steps 9 to 13 until no economy is achieved either 
for altering the beam sections or for altering the 
ratios of internal to external columns. At the end of 



this step the frame stability is satisfied and the 
deflections and sectional properties calculated fully 
satisfy the stiffness equations„ 

Step 15 - Check the stress requirements in the frame. This is 
done by; 

a - checking that the plastic modulus of the beams 
under combined vertical and wind loads are satisfactory. 
With the unit load factor X = 1 no plastic hinge should 
develop in any beam and the entire frame must be elastic 
under such a load factor. See details later, 
b - Checking that under factored loads, X > 1 (say 
X = 1.29 as in the recent codes), no hinge develops in 
any column. See details later. 

If the plastic modulus of all the sections are satisfactory, 
which is often the case, the design process is complete. If not, 
select new beam and column sections for the members and repeat 
steps 5, 6 and 7. Notice that if I^ of the columns are known the 
iteration at this stage calculates the horizontal joint deflections 
instead. Thus during each round of iteration either column sections 
are unknown and calculated with given specified values of the hori¬ 
zontal storey deflection or these deflections are calculated using 
column sections selected to satisfy strength requirements. 

It is evident that the values of the horizontal deflection at 
each storey can be specified in more than one way. In the following 
sections a method is proposed for this purpose so that the recent 
suggestions made by the British Standard documents are satisfied. 

2.9 GENERAL EQUATION FOR THE SWAY DEFLECTION CONSTRAINTS 

The draft of the British Standards Institution B/20 document 
77/13908 DC restricts the differential deflection between the two 





ends of any column to not more than a constant value. This constant 
may be taken as aH = H/300, where H is the length of the column. 
Excessive differential deflection causes cracks in the walls, 
cracks in the cladding and may also damage windows. If the sway 
in each storey is the maximum of H/300 then the deflection profile 
of a multistorey frame is a straight line such as OA in Figure 
2,8, In practical frames, the sway in the lower columns is less 
than those in more flexible upper columns and thus a linear profile 
cannot take place. It is therefore necessary to derive a non¬ 
linear deflection profile for the frame which represents the defor¬ 
mation of the frame more realistically. Measuring x from the top 
of the frame, let the deflection y at a distance x from the top be 
given by the deflection function, 

y = a^ + a^ X + a^ + a^ x^ + a^ x^ (2,46) 

where a^, a^, a^, a^ and a^ are constants and must satisfy the 
boundary conditions at the top and the bottom of the frame and n 
is a variable which for economic reasons should be numerically as 
high as possible. As n increases, the frame becomes more flexible 
and the non-linear profile of the frame approaches the linear 
upper bound OA, The constants a^ to a^ are found as follows; 

1 - Since in a practical frame each floor is required to be 

stiffer than or as stiff as the floor above, it follows that 
the sway in each storey must not be more than the sway in the 
storey above. In this thesis■it is proposed that a sway.of 
H/300 should be allowed only in the top storeys. Thus the 
slope of the deflection profile, curve OB in Figure 2,8, is 
made equal to a i,e, when x = 0 dy/dx = a. 
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dy/dx = + 2a^ x + 3a^ x^ + na^ x^”^ 

and when x = 0, dy/dx = a 


a = a^ + 0 


and a^ = a = 1/300 


2 - The top end of the frame is free and considering the whole 
frame as a single element, the bending moment at this free 
end will be zero. Thus wlien x = 0 d^y/dx^ = 0. Now ' 

d^y/dx^ = 2a^ + 6a^x + n(n-l) a^x^ ^ 

and when x = 0 d^y/dx^ = 0 



3 - At the free end, at the top of the frame, the shear force is 

also equal to zero. Thus when x = 0 d^y/dx^ = 0 and thus 

aq = 0. 

4 - Because of the fixed end condition at the ground level, we have 

\\/hen X = L, y = 0.0. Thus equation (2.46) gives: 

0.0 = a^ + aL + a^L^ (2.47) 

where L is the total height of the frame. 


5 - Furthermore, for the condition of fixity at the base, when 
X = L, dy/dx = 0.0. This leads to 

0.0 = a + na^L^ ^ (2.48) 

Solving equations (2.47), and (2.48) gives the value of a^ 
and a^ as: 


(2,49) 





' ' ' - '5'= V ‘ -- * 
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n L 


n-^1 


C2,50) 


Equation (2.46) now becomes: 


a 


(k+3) 


k+4 


- L 


Ck+3) k+3 .0^+4) 


X + 


k+4 


] 


(2.51) 


where n = k + 4. As k tends to infinity, equation (2.51) tends 

to the ultimate straight line deflection, OA in Figure 2.8. In 

lc+4- 

practice, k is taken as large as possible. Provided that x 
does not lead to an exponential overflow. During the calculation 
of the deflection for each storey using equation (2.51). The 
value of k at which the overflow occurs changes from one frame 
to another. For each frame therefore the value of k is taken 
as a number which is smaller than that which causes an exponential 
overflow. Equation (2.51) ensures that the sway in each storey 
will be less than the sway in the storey above and this sway 
reaches H/300 only at the top floor. For a large value of k, the 
profile OB is nearly parallel to OA except for the bottom part of 
OB. 

In the special case, where m storeys of the frame are of equal 
height h, equation (2.51) can be written in terms of the storey 
level r Figure 2.8. This gives: 


y 


r 


ah 


(k+3) 


m 


(m-r) 

k+4 


k+4 


k+3 

m 


(m-r) 


+ 


(k+3) k+4 

7T— ro 
(k+4) 


(2.52) 


where y is the sway of the storey at level r. 
r 

2.10 UNIVERSAL BEAM AND COLUMN SECTIONS USED IN THE DESIGN 


In multi-storey steel frames it is normal practice to use the 
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available rolled sections unless, in special cases, built-up 
sections are found to be Tnore economical. The computer program 
written for the design of sivay frames selects universal column 
sections for column members, and universal beam sections for beam 
members. Because of local instability problems, certain sections, 
depending on the grade of the steel used are not suitable,(Horne, 
1964) and are not included in the final selection. For mild steel 
(grade 43) and high yield steel grade (50) the universal column 
and beam sections (BCSA, CONSTRADO, 1971) are reorganized in the 
order of their cost. The high yield steel (grade 55) is not used, 
because of the high probability of local buckling. Furthermore, 
in multi-storey work. Grade 50 is frequently the most economic 
(Needham, 1977). Expensive sections are omitted from the list. 

The prices of the universal column sections and beam sections are 
taken from the price list of steel sections (British Steel Cor¬ 
poration, 1977). A section is considered expensive if there is 
another less costly section with an equal or larger second moment 
of area. The sections with the minimum price-inertia ratio are 
selected. Tables 2.2 and 2.3 show the economical beam and column 
sections respectively, for mild steel grade 43; and Tables 2.4 and 
2.5 are for high yield steel grade 50. 

For high-rise multi-storey frames, universal sections are not 
adequate. For this reason a few built-up sections are added. These 
are fabricated from the largest universal section with two rein¬ 
forced steel plates. 

To satisfy the deflection requirements the beam sections (cal¬ 
culated in section 2.5) are selected from the list before designing 
the column sections. Thus the discrete available beam sections are 
made full use of during^ the process of satisfying the deflection 


requirements. 





2oll ALTERING THE BEAM SECTIONS 

The selected beam sections (Section 2,5) do not necessarily 
result in the most economic design for the complete frame. For 
this reason, a linear extrapolation is used to alter the beam sec¬ 
tions in an effort to reduce the cost of the frame. To apply 
this extrapolation, two preliminary designs are needed, and initi¬ 
ally two sets of beam sections must be provided for these two 
designs. The beam sections selected in Section (2,5) upgraded by 
one provide the first set of beams. The second set of beam sec¬ 
tions are provided by those obtained in Section (2,5), The 
iteration technique is entered with the first set of beams to obt¬ 
ain the sectional properties of the column, deflections and 
rotations, so producing the first cycle in the design process. 

For the second design the iteration technique is entered again with 
the second set of beams. The cost of the material for the frame 
in the first and second cycles is found, and used in an extra¬ 
polation process to reduce cost. In this process 1^^ and I^^ are 
the second moment of area for the beam and the column in a given 
storey for the first cycle, 1^2 and 1^2 those for the second 
cycle. If Ig 2 is smaller than 1^2 by section, then Figure 

2,9 shows that the second moment of area I^^ ibe column section 
is given by: 


^C3 " ^C1 "" ^^C2 " ^Cl^ ^^Bl ” ^B3^^‘^^B1 “ ^B2^ 


(2,53) 


(2,54) 


Similarly, this equation can also be used if the beam section for 
each storey is upgraded by one section to find the expected coliamn 


section. 






It must be emphasised that the extrapolation is introduced to 
avoid solving the modified stiffness equations when’beam section 
is increased or decreased„ The value obtained for the second 
moment of area of a column is approximate and is used only for 
costing purposes„ However, it was found that the value of the 
second moment of area for a column obtained using equation (2„54) 
is very near the value obtained by solving the modified stiffness 
equations. 

The approximate cost of the material at each storey (that is 
the cost of the beam and the columns supporting it) is calculated 
by extrapolation for the following: 

1 - for the first design cycle; 

2 - for the second design cycle; 

3 - for the case when the beam is downgraded by one section; 

4 - for the case when the beam is upgraded by one sectioni 

The beam section which gives the lowest cost is chosen for the 
third design cycle. This process is continued for each storey 
until successive design cycles will show that a decrease or an 
increase in any beam section would not lead to a cheaper design of 
the frame. 

To avoid reverse column taper, a check is made in each design 
cycle to ensure that each column section is less than or equal to 
the column below it. If a current design cycle results in reverse 
column taper in a given storey, the beam section of the previous 
cycle is adopted, and no reduction on this beam size is allowed, 

2,11,1 Example on Reducing Beam Sections 

The four storey single bay frame shown in Figure 2,10a is 
designed as an example. This frame has been designed by (Majid, 
1972) using the elastic plastic design method for sway frames. 






This showed that the sway deflection of the frame was somewhat 
excessive. This is in spite of the fact that the frame itself 
had a reasonable reserve of strength over and above that required 
by (BS 449,1969), The result of this design is summarized in 
Figure 2„10b which shows the second moment of area for each member 
after selecting sections from the list of rolled sections, and 
in Figure 2,10c which shows the order of the formation of hinges 
under combined loading. It can be seen that collapse takes place 
at a load factor of 1,49, which is above an adopted permissible 
value of 1,4, However, the maximum deflection for this design is 
53,8 mm which is more than allowable, ITiis frame is used here for 
an extensive investigation into the design method proposed in this 
thesis. Although the reduction in the cost of the material is con¬ 
sidered in the actual computer program, a reduction in weight is 
used for this example to compare the results with those given by 
an optimum design program based on weight optimisation (Saka, 1975) „ 
The axial force effect and the axial stiffness are not considered 
in this example for similar reasons. 

Application of equation 2,52 for various values of k is 
shown in Figure 2,11. The value of k just before an overflow 
occurred was 56; the deflection profile for this value of k is used 
to specify the horizontal deflection at each beam level. It can be 
seen in Figure 2,11 that for k = 56 the non-linear and the linear 
profile are nearly the same. Any difference between the adopted 
deflections and those obtained by the straight line are in the 
bottom storeys, 

i) First Design cycle (Figure 2,12) 


The second moment of area of the beam sections for the first 





design cycle are found as follows: - : ; 

1) The second moment of area for the sections given by 
beam mechanism are found with a load factor of 
1,75. This gives I = 21345 cm^, 

2) The values obtained in (1) above are multiplied by 
a magnification factor. This factor can take any 
value between 1,5 and 2, In this example this 
factor is taken equal to 1,75 for all beams (i,e, 

1.75 X 21345 = 37354), 

3) Select sections from the universal beam sections 

4 

this gives I = 40414 cm , 

4) The sections obtained in (3) above are upgraded 

4 

by one section which gives I = 55779 cm 
These values of I are used in the iteration using the modified 
stiffness equations. The result is summarized in Figure 2,12, 

The second moment of area for the beams (chosen from the list) 
is shown on the top of each beam. The calculated and the selected 
second moments of area for the columns are shown next to the columns, 

ii) Second Design Cycle (Figure 2.13 and 2,14) 

The second moment of area of the beam sections for this design 
cycle are those of the first cycle downgraded by one section (i.e. 

I = 40414 cm^). The result of the iteration using the modified 
stiffness equations is shown in Figure 2.13, This shows that 
columns number 9 and 10 have a section larger than that of columns 

II and 12 in the bottom storey. For this reason this design is 

unacceptable. To improve it, beam 3 was increased to the next lar- 

4 

aer size in the table with I = 55779 cm , 'Fhe iteration gives 
columns 9 and 10 again larger than the columns below. Therefore, 



beam 3 is increased once more to the next larger size which has 

4 

I = 75549 cm o Figure 2,14 shows the results with beam 3 having 

4 4 '■ 

I - 75549 cm and all other beams having I = 40414 cm „ This 

modified design replaces the one in Figure 2,13, It can be seen 

that the design in Figure 2,14 is lighter than the design in Figure 

2.12. That in Figure 2.14 is used to start the extrapolation for 
the third design cycle, 

iii) Third Design Cycle (Figure 2,15) 

To obtain the second moments of area for the beam sections 
the procedure of section 2,11 is applied. Each beam section (except 
beam 3) is upgraded and downgraded by one section and the weight 
of each storey is found with the column sections obtained from 
equation 2,54, In this equation 1^^ and are those of Figure 

2.12, while 1^2 ^.nd 1^2 taken from Figure 2,14, The beam 
which gives the least weight in each storey is used in the iter¬ 
ation, The result is shown in Figure 2,15, The design of the 
top storey is altered to avoid beam mechanism collapse, due to the 
formation of hinges at the top ends of the columns. The top columns 
are given the same sections as the columns below, and .the beam is 
designed to satisfy strength requirements. This gives I = 25464 
cm"^, no further reduction on this beam section is possible, 

vi) Fourth Design Cycle (Figures 2,16 and 2,17) 

Further reduction on the beam sections of the third cycle 
gives the design shown in Figure 2.16, This is unacceptable, bec¬ 
ause the I values for columns 9 and 10 are larger than those for 
columns 11 and 12, To avoid reverse column taper the I value for 
beam 4 in the fourth design cycle is kept as it is in the third. 

The fourth design cycle is shown in Figure 2,17, 








This is the final design cycle, the result isshown in Figure 

2„18„ No further reduction in beam sections could be made. The 

reason is that any reduction in the bottom two beams will result 

in reverse taper column, while any change in the top two beams 

will violate the strength requirements for vertical load. 

Figure 2,19 shows the result obtained by Saka (1975), using 

an optimisation method. The second moments of area for the beam 

sections obtained by Saka are shown above each beam, while the 

available second moment of area for the sections selected from 

the list is given below the beam. To calculate the weight of the 

frame the length L of each member is multiplied by the second 

moment of area I, and by the mass per unit length m„ 

The sums EIL and ZmL are calculated. It was found that ZIL 

for the final design proposed in this thesis is 1,9032 x 10^^ mm^, 
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while for the optimum design is 1,8835 x lo mm , The percentage 
difference being 1,03%, However, using the list of available 
sections for both designs it was found that the total mass for the 
author's final design was 5349,0 kg which is less than that for the 
optimum design (5497,3 kg) by 2,77%, This difference is solely 
due to the discrete nature of the beam sections available. The 
computer time needed to complete the design by the present method 
was 1/20 of that needed for the optimum design. The storage 
needed was 1/10 of that needed by the optimum design. 

Finally, to summarize the design cycles given in Figure 2,12 
to 2,18 a graph is drawn for the total mass of each design cycle 
in Figure 2.21, This shows the reduction in mass in each step. 

The increase in the total mass in the third design cycle is due 
to the increase in the top column sections to avoid the formation 




of plastic hinges at the top ends of these columns. 


2ol2 • SELECTING THE INTERNAL TO EXTERNAL COLUMNS RATIO 

Equation (2.34) gives the ratio of the second moments of area 
of the internal and external columns as might be obtained using a 
simple design method. This ratio is not the most economical one, 

A procedure to obtain a more economical distribution of inertia 
between internal and external columns is now considered. 

In the first two preliminary design cycles, equation 2,34 
gives the ratio between the internal and the external second mo¬ 
ment of area of the columns. After calculating the columns iner¬ 
tia and then selecting a section for each column from the list, 
the selected sections will have a ratio different from the one 
obtained from equation 2,34, The new ratio between the selected 
sections might be more economical to use. Therefore, a new design 
is examined for cost, the design having the same set of beams, 
but with the new ratio. The flow diagram shown in Figure 2,21 
summarizes the procedure for economy in selecting beam section, and 
internal to external column ratio. 
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TABLE 2.1: SHOWS THE CONVERGENCE OF THE DEFLECTIONS AND SECOND MOMENTS OF AREA 
OF THE COLUMNS TO THE DESIGN VALUES, WITH A TOLERANCE OF 0.01 














EAR AND NON-LINEAR DEFLECTION PROFILES 
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Second moment 

Area 

Price £ per m 


UB sections 

of area cm^ 

cm^ 

914 

X 

419 

X 

388 

1911319 

996„7 

185.54 

+ 2 

*“ 

oU mm 

plates 


914 

X 

419 

X 

388 

1480449 

829,1 

156.00 

+ 2 

“ 

40 mm 

plates 


914 

X 

419 

X 

388 

1082870 

661.5 

126,45 

+ 2 

- 

20 mm 

plates 



914 

X 

419 

X 

388 

717325 

493.9 

86,91 

914 

X 

419 

X 

343 

623866 

436.9 

76.83 

914 

X 

305 

X 

289 

503781 

368.5 

64,30 

915 

X 

305 

X 

253 

435796 

322,5 

56.29 

914 

X 

305 

X 

224 

375111 

284.9 

49.84 

914 

X 

305 

X 

201 

324715 

256.1 

44.72 

838 

X 

292 

X 

194 

278833 

246.9 

42.68 

838 

X 

292 

X 

176 

254412 

223.8 

38.72 

762 

X 

267 

X 

173 

204747 

220.2 

37,80 

762 

X 

267 

X 

147 

168535 

187.8 

32.12 

686 

X 

254 

X 

140 

135972 

178.4 

30,59 

686 

X 

254 

X 

125 

117700 

159,4 

27.31 

610 

X 

229 

X 

113 

87260 

144.3 

24.46 

610 

X 

229 

X 

101 

75549 

129.0 

21.87 

533 

X 

210 

X 

92 

55225 

117.6 

19.64 

533 

X 

210 

X 

82 

47363 

104.3 

17.51 

457 

X 

191 

X 

82 

37039 

104.3 

16.56 

457 

X 

191 

X 

74 

33324 

94.9 

14.95 

457 

X 

191 

X 

67 

29337 

85.4 

13.53 

457 

X 

152 

X 

60 

25464 

75.9 

12.93 

457 

X 

152 

X 

52 

21345 

66.5 

11.21 

406 

X 

178 

X 

54 

18576 

68.3 

10.99 

406 

X 

140 

X 

46 

15603 

58.9 

9.98 

406 

X 

140 

X 

39 

12408 

49.3 

8.46 

356 

X 

127 

X 

33 

8167 

41.7 

7.26 

305 

X 

102 

X 

28 

5415 

36.3 

6.22 

305 

X 

102 

X 

25 

4381 

31.4 

5.55 

254 

X 

102 

X 

22 

2863 

28.4 

4.86 


TABLE 2„2 


UNIVERSAL BEAM SECTIONS FOR GRADE 43 STEEL 







UC sections 

Second moment 
of axea cm^ 

Area 

cm^ 

Price £ per m 

356 
+ 2 

X 

406 X 
60 mm 

634 

plates 

639385 

1316,9 

241,74 

356 
+ 2 

X 

406 X 
40 mm 

634 

plates 

499630 

1147,3 

211,86 

356 
+ 2 

X 

406 X 

20 mm 

634 

plates 

378668 

977,7 

181,95 

356 

X 

406 

X 

634 

275140 

808,1 

142,02 

356 

X 

406 

X 

551 

227023 

701,8 

123,42 

356 

X 

406 

X 

467 

183118 

595,5 

104,61 

356 

X 

406 

X 

393 

146765 

500,9 

88,03 

356 

X 

406 

X 

340 

122474 - 

432,7 

76,16 

356 

X 

406 

X 

287 

99994 

366,0 

64,29 

356 

X 

406 

X 

235 

79110 

299.8 

52,64 

356 

X 

368 

X 

202 

66307 

257,9 

44,95 

356 

X 

368 

X 

177 

57153 

225,7 

39,38 

356 

X 

368 

X 

153 

48525 

195,2 

34,04 

305 

X 

305 

X 

158 

38740 

201,2 

32,86 

305 

X 

305 

X 

137 

32838 

174,6 

28,50 

305 

X 

305 

X 

118 

27601 

149,8 

24,54 

254 

X 

254 

X 

107 

17510 

136,6 

22,26 

254 

X 

254 

X 

89 

14307 

114,0 

18,51 

254 

X 

254 

X 

73 

11360 

92,9 

15,18 

203 

X 

203 

X 

71 

7647 

91,1 

14,91 

203 

X 

203 

X 

60 

6088 

75,8 

12,60 

203 

X 

203 

X 

52 

5263 

66,4 

10,92 

152 

X 

152 

X 

37 

2218 

47,4 

7,97 

152 

X 

152 

X 

30 

1742 

38,2 

6,47 


TABLE 2.3: UNIVERSAL COLUMN SECTIONS (GRADE 43) 






UB Sections 


Second moment 
of area cm^ 


Price £ per m 


As in Table 2„2 for steel grade 43 



838 

X 

292 

X 

194 

278833 

246,9 


762 

X 

267 

X 

173 

204747 

220,2 


762 

X 

267 

X 

147 

168535 

187,8 


686 

X 

254 

X 

140 

135972 

178,4 


610 

X 

229 

X 

140 

111673 

178,2 


610 

X 

229 

X 

125 

98408 

159,4 


610 

X 

229 

X 

113 

87260 

144,3 


533 

X 

210 

X 

109 

66610 

138,4 


533 

X 

210 

X 

101 

61530 

129.1 


533 

X 

210 

X 

92 

55225 

117,6 


457 

X 

191 

X 

89 

40956 

113,8 


457 

X 

191 

X 

82 

37039 

104,4 


457 

X 

191 

X 

74 

33324 

94,9 


457 

X 

191 

X 

67 

29337 

85,4 


457 

X 

152 

X 

60 

25464 

75,9 


457 

X 

152 

X 

52 

21345 

66,5 


406 

X 

140 

X 

46 

15603 

58,9 


356 

X 

127 

X 

39 

10054 

49,3 


356 

X 

127 

X 

33 

8167 

41,7 


305 

X 

102 

X 

28 

5415 

36,3 


305 

X 

102 

X 

25 

4381 

31,4 


305 

X 

102 

X 

22 

2863 

28,4 


TABLE 2 A: 


UNIVERSAL BEAM SECTIONS FOR GRADE 50 STEEL 









UC Sections 


Price £ per m 


Second moment 

Area 

of area cm'^ 

CFi^ 


As in Table 2„3 for steel grade 43 


356 X 368 X 202 

66307 

257„9 : 

44.95 

305 X 305 X 198 

50832 

252.3 

41,18 

305 X 305 X 158 

38740 

201.2 

32.86 

305 X 305 X 137 

32838 

174.6 

28,50 

254 X 254 X 132 

22416 

167.7 

27,46 

254 X 254 X 107 

17510 

136,6 

22.26 

254 X 254 X 89 

14307 

114.0 

18.51 

203 X 203 X 86 

9462 

110.1 

18,06 

203 X 203 X 71 

7647 

91.1 

14.91 

203 X 203 X 60 

6088 

75,8 

12.60 


TABLE 2o5: UNIVERSAL COLUMN SECTIONS FOR GRADE 50 STEEL 
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CURE 2.14: MODIFIED SECOND DESIGN FIGURE 2.15: TMIR: 

CYCLE 



























L = 1.9032x10 ‘ min^ (using I as ETL - 1.8835 x 10 ' mm 

leu1at e d) calculated) 

ss of beams 2552.8 kg Mass of beams 2854.8 k 

ss of columns 2796.2 kg Mass of columns 2642.5 

tal mass 5349.0 kg Total mass 5497.3 kg 

























BEGIN 

V y 



Choose the cheaper preliminary deign to start 


Change internal to external column ratio keeping the previous 
beam sections 


the design with the new ratio 
cheaper? 


NO Keep the 
old ratio 


Adopt the new ratios 


Change beam sections keeping the previous column ratio 


Adopt the 


new oeam 


sections 


the design with the new beam 
section cheaper? 


FIGURE 2.21 THE PROCEDURE FOR ECONOMY 















CHAPTER THREE 

REFINING THE DESIGN PROCESS 

3„1 REDUCTION IN THE BENDING STIFFNESS OF THE COLUMNS 

In the presence of high member axial forces the bending stiff¬ 
ness of a column is reduced by including the "stability functions" 
in the stiffness equations of the columno If Livesley's ()> stability 
functions are used, equations (2.5) are altered to become: 
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Accordingly, the modified stiffness equation for a joint in a 
column (given by equations 2.20, 2.23, 2.26, 2.29, 2.31 and 2.38) 
are altered. For example, equation (2.20) becomes: 


X =rH.+d.(i>^-i‘.+b.(b^.x + 

j ^ j C3 ^2 c3 3 C 3 ^ 5 c 3 n 

d . d)^ . r + b (pr- X - d r (■ 

cj ^2 c 3 n cm ’^ 5 cm m cm 2cm m 

-d d) r.)/(b - ^ ) 

cm ^2cm 3 ^ ^ c3 ^5c3 cm ^5cm 

where 3 and m refer to member numbers. 

Obviously, the values of the (j) functions are dependent on the 


second moment of area 


of the column but this is not known before 










88 


the iteration starts. Furthermore, during the first design cycle 
the values of I for the column change rapidly as the iteration is 
started with unacceptably small values of I, For these reasons, 
and since the first design cycle is merely a procedure to obtain 
approximate sections, the effects of the axial forces are neglected 
in this cycle. Tlius (j)^ = (((^ = (j)^ = (|)^ = ^^ = 1„ 

At the beginning of a subsequent design cycle an approximate 
value of I is available for each column. This is used for cal¬ 
culating the (}) function in the first iteration. At the end of each 
iteration, a new I is obtained for each column which is then used 
in the next iteration to calculate a new and revised set of stabi¬ 
lity functions, 

3.2 CALCULATION OF THE AXIAL FORCES 

The values of the stability functions are dependent on the 
axial force in the column which is also not known before the 
iteration starts. However, to assess the axial forces it is found 
that it is sufficiently accurate to proportion the vertical loads 
between the columns so that the external columns equally carry 
half the load on the external bays above. On the other hand, 
internal columns proportionally share the rest. The use of this 
approximation in assessing the axial forces instead of calculating 
the exact ones simplify the design procedure and reduce the computer 
time. Alternatively, it is possible, using a computer, to cal¬ 
culate the axial forces exactly from: 

{p} = k . A . {x} (3.3) 

where {p} is the vector of member forces, is the member stiffness 
matrix, A is the displacement transformation matrix of the member. 








and {x} is the vector of joint displacements. 

To obtain the exact axial forces for each design cycle the 
iteration must be used first to find the joint displacements. The 
internal forces are then calculated using equations C-^.3), This 
process must be repeated until the difference between two succes¬ 
sive sets of axial loads is smaller than a specific tolerance. 

Such a calculation is lengthy and it can be shown that the differ¬ 
ence between the approximate and the exact values of the axial 
forces is small. Their effect on the outcome of the design is even 
smaller. 

3,3 EXAMPLE ON THE AXIAL LOAD EFFECT 

The five-storey, three unequal-bay, frame shown in Figure 3,1 
was designed using the deflection profile shown in the figure. 

This frame is used here for an extensive investigation into the 
effect of the axial load in the design of sway frames. Two designs 
were obtained. The first design procedure initially neglected the 
effect of axial forces in the columns, until a reduced cost 
design was obtained, Tlie beams of this design were kept unaltered 
but the columns were subsequently altered to allow for the effect 
of axial forces. These are calculated exactly using equations (3.3), 
with a tolerance of 0.001, Second column in Table 3,1 shows the 
reduced cost design with the axial forces neglected, wnile the third 
column in the table shows the design with the axial load effect 
included. In both designs the second moments of area after selecting 
sections from the list are shown. It can be seen that some column 
sections are increased due to the effect of the application of the 
axial load. These sections are marked by an asterisk in Table 3,1. 
Including the axial load effects increases the material cost by 4.09%, 













Member 
number 
(taken as 
joint 
number 

Final linear 
design axial 
load effect 
is neglected 

Final Non-linear 
design. The lin¬ 
ear design is 
modified by con¬ 
sidering the 
exact axial load 

I in cm"^ 

Final Non-linear 
design, The_effect 
of the approximate 
ax;ial load is con¬ 
sidered from first 
design cycle 

I in cm4 

1-3 

47363 

47363 

33324 

4-6 

47363 

47363 

47363 

7-9 

21345 

21345 

29337 

10 - 12 

21345 

21345 

18576 

13 - 15 

18576 

18576 

18576 

16 

11360 

11360 

11360 

17 

14307 

17510* 

17510 

18 

11360 

14307* 

14307 

19 

6088 

7647* 

7647 

20 

11360 

11360 

11360 

21 

14307 

17510* 

17510 

22 

11360 

14307* 

14307 

23 

6088 

6088 

7647 

24 

5263 

6088* 

7647 

25 

11360 

14307* 

11360 

26 

5263 

6088* 

11360 

27 

5263 

6088* 

5263 

28 

5263 

6088* 

5263 

29 

7647 

7647 

7647 

30 

5263 

6088* 

6088 

31 

5263 

5263 

5263 

32 - 35 

5263 

5263 

5263 

Total cost 
of mater- 
ial in £ 

£2485„95 

£2592,03 

£2591.13 

Total mass 
in kg 

1 

11778 kg 

12286 kg 

12519 kg 


* marks sections change due to the effect of exact axial force 


TABLE 3ol: THE DESIGNS OF THE FIVE STOREY FRAME 













In the second design procedure the axial load effect is 
introduced from second design cycle as described in section 3,2, 
The axial forces for this design are calculated approximately 
as explained in section 3,3. The design obtained in this manner 
is shown in column four of Table 3,1. Although this design is 
different from that obtained by considering the exact axial force, 
the cost of the material for both designs is almost exactly the 
same. However the total mass is different by 1,9%, This shows 
that a consideration of a cost reduction in the design procedure 
is much more realistic than a consideration of a weight reduction. 

When considering the design with approximate axial forces, 
column 25 and 27 are found to have smaller sections than those 
obtained when the axial forces are calculated exactly. On the 
other hand beams 7, 8 and 9 have larger sections, see Table 3.1. 

From this example one could conclude the following two points 

1 - Taking axial load effect into consideration increases the 

initially selected column sections 

2 - Calculation of the axial loads approximately, simplifies the 

design, procedure but produces feasible designs which are not 
much different from those obtained using an exact procedure, 

3.4 economy in computer TIME 

Normally between five and twelve design cycles are needed to 
obtain the final set of sections for a frame. Although the iter- 
process in each cycle converges puicxly it is poscsible to 
reduce the computer tim.e in two different ways. These are: 

(a) Once the iteration in a design cycle is complete, the 
joint deflections, the axial forces and the second moment of area 
of each column becomes known. Before entering the next design 





cycle, a set of heasi sections are selected in an effort to reduce 
costo This set is only slightly different from the one used in the 
previous cycle as each beam is only up or down graded by one sec¬ 
tion „ Such a change does not alter the member axial forces and 
joint deflections drastically. Rather than starting the iteration, 
in the new design cycle, from unrealistic initial values for all 
the variables in the modified stiffness equations, it was decided 
to use the old values of these variables as given in the previous 
design cycle. Starting the iteration using these realistic ini¬ 
tial values, speeded up the subsequent design cycles considerably. 

For example, eight design cycles were needed to obtain the 
final reduced cost design of Table 3,1, column four. The number 
of iterations for each trial is shown in Figure 3,2, where it can 
be seen that the number of iterations in each cycle is reduced 
drastically after the second design trial. The reason is that, 
the results of one design cycle are used as the starting point for 
the next design. The large number of iterations in the second 
design trial is due to the introduction of the stability functions 
at this stage. The total number of iterations for the last six 
trials is less than that of the first or the second trial, 

(b) As an iterative technique is being used for solving the 
modified stiffness equation, the tolerance control of this technique 
has a vital influence on the speed of convergence, A small toler¬ 
ance increases the number of iterations and results in a consider¬ 
able increase in computer time. This is not only unnecessary but 
also can consume all the available computer time before convergence 
is achieved, A large number of tests showed that a tolerance of 
0.01 was adequate for relatively small frames (up to twelve storeys), 
and also for medium size frames (up to twenty-four storeys) which 










had equal or nearly' equal hay's. For raediujii irregular frames, with 
unequal bay's, it was found that a tolerance of 0,001 will be adequ¬ 
ate, This is also sufficient for larger frames. To demonstrate 
the effect of the tolerance on the design outcome, the five-storey 
three unequal bay frame of Figure 3.1 was first designed with a 
tolerance of 0.001. The results of this design are shown in 
Table 3,1, column four. This same frame was then designed using 
a tolerance of 0,01 and resulted in the same final sections. To 
compare the displacements of the joint, this design was analysed 
by an existing computer program which solves the stiffness equations 
explicitly (Celik, 1977). Comparison between sway deflections 
is not possible because these are constants in the design whatever 
the tolerance. For this reason the vertical displacements are 
compared . For joints 1, 2 and 3, see Figure 3,1, these are listed 

in Table 3o2o 


Kinds of Solutions 

y for joint 1 

y for joint 2 

y for joint 3 

Direct Solution 

-12,42 mm 

-2.61 mm 

-3,78 mm 

Iterative Solution with 
tolerance of 0.001 

-12,42 mm 

-2.61 mm 

-3,78 mm 

Iterative solution with 
tolerance of 0.01 

-12.40 mm 

-2.64 mm 

-3o79 mm 


TABLE 3 2 VERTICAL DISPLACEMENTS OF JOINTS 1, 2 AND 3 

OF FIGUTE 3.1 

It can be seen from the table that the direct solution of the stiff, 
ness equations leads to the same answer as that obtained by the 









iterative technique with, a tolerance of■*0,001, For the tolerance 
of 0,01, the displacements are nearlv the same, and the difference 
does not exceed 0,05 mm. 

The computer time needed, with a tolerance of 0,01, using 
Aston Computer (ICL 1904s) to complete the eight design trials was 
only 26 seconds. This included the time to check the strength of 
the final design. The efficiency of the design method is realised 
when it is noticed that a single explicit non-linear analysis of 
this frame, using a compact storage scheme, needed as much as 22 
seconds. The dotted graph in Figure 3,2 shows the number of iter¬ 
ations for each trial with a tolerance of 0,001, As expected, the 
number of iterations for each trial increased, and the computer 
time increased to 46 seconds. However, the final reduced cost 
design is the same. Thus, in this example, no benefit is gained 
using a small tolerance. Even then the computer cost of a com¬ 
plete design is only twice as much as an explicit analysis solving 
L = IK X directly. 

.3.5 ECONOMIC USE OF COMPUTER STORAGE 

As stated in Chapter 2, sections 2,4, 2,4,1 and 2.4,2 the 
modified stiffness equations are formulated and solved iteratively 
one joint at a time. Thus no matter how large the frame is, the 
storage requirements is that needed by a single joint. The com¬ 
puter program written for the design method does not therefore 
require the use of any backing store as is often the case in the 
analysis of large frames. Five one-dimensional arrays are needed 
in the core store. These keep the current values of areas and the 
second moments of area of the sections, the current ratio of inter¬ 
nal to external columns, the current rotations of the joints, and 



■ ■ ■■■■■ ^ ”|p|||^ 

their vertical displacements, A further, five similar arrays are 
also needed to store the previous values of these variables. Thus 
the storage requirement of the design method is very small indeed. 
For a total of n joints in a frame, the iteration process only 
required 10 n locations. For instance, the five storey, three bay 
frame shown in Figure 3,1 used only 8050 words (approximately 
8 k). This included storing the computer program itself, 

3,6 THE DESIGN CRITERIA 

Two categories of design criteria are adopted. The first 
category includes the deflection requirements, while the second is 
concerned with the strength requirements. For deflection require¬ 
ments, the existing BS 449 does not recommend any specific value 
as the limit for the horizontal sway in the columns of multi-storey 
frames. On the other hand the B/20 Draft (Document 77/13908DC, 
1977) which will replace BS 449 recommends that: 

**The sway of the column should not exceed the limit of 
the height of the storey/300". 

This design criterion dominates the design of the sway frame when 
the wind load is high. In reality the wind load is calculated 
according to CP3, Chapter five (British Standard 1972), If the 
frame is subjected to such a wind load, and is made from grade 50 
steel, then the horizontal sway in the columns usually governs its 
design. Frames made from grade 43 steel are also dominated by the 
column sway in many cases. Thus one criterion is 

i) The sway of the storey does not exceed its height/300. 

The second category of the design criteria caters for strength 
requirements. For this B/20 Draft Specification in Chapter 12, 
clause 12,5,4,1 recommends the following: 



"In ^ multi-storey rigid jointed frame which is not effectively 
braced against sidesway in its own plane, plastic design shall 
be used only when proper allowance is made for frame-insta¬ 
bility effects. This may be done either by carrying out a full 
elastic-plastic sway analysis using methods which make proper 
allowance for frame instability or alternatively, in appro¬ 
priate cases by the simple method specified in 12„5,4„3"o 

The method of Majid (1972), (Majid and Anderson 1969), and 
(Anderson, 1969) is one of these methods which makes proper allow¬ 
ance for frame instability. The overall instability of the frame 
due to gradual deterioration of its stiffness is fully considered. 
For this reason, most of their design criteria are adopted; these 
ensure that the frame has adequate load factors against collapse 
under various loading cases. The aim is to ensure that the two 
following conditions are satisfied, 

(ii) Under combined vertical and wind load from either side, 
the frame does not collapse until the load factor exceeds 

(taken as 1,29 in this thesis), 

(iii) Under dead load and vertical load the frame does not 
collapse until the load factor exceeds (taken as 1,75), 

Analysis of a large number of frames by Majid (1972), (Majid 
and Anderson 1969) and (Anderson, 1969) shows that if certain res¬ 
trictions are made on the formation of the plastic hinges, then 
in most cases the above two conditions are not violated. These 
restrictions are: 

(iv) No plastic hinge should develop in a beam below the 
load factor of unity and the frame should be entirely elastic under 
the specified unfactored loads, 

(V) No plastic hinge shall develop in a column until the load 
factor exceeds a permissible value A^ under combined loading or A 2 
under vertical loading. 

The condition (iii) and the restrictions on the formation of 
(iv) and (v) are used as design criteria for 


the plastic hinges 







this work. Condition (ii) is excluded because.it is necessary 
to perform an elastic-plastic analysis, which follows the successive 
development of plastic hinges. This kind of analysis needs consider¬ 
able computer storage and time. However, the restriction on the 
formation of plastic hinges is a safeguard for this criterion. 
Furthermore, the instability of the frame increases rapidly with 
the increase of the sway of the frame. For this reason the restric¬ 
tion on sway imposed by criterion (i) will provide the proper 
allowance for frame instability proposed in Clauses 12,5.4,1 of B,20, 

To ensure that the above criteria are all satisfactory, the 
final design is analysed at a load factor = 1,29 under combined 
loading and taking the effect of axial loads into consideration. 

The maximum bending moment in each column obtained by this analysis, 
is compared with the reduced plastic hinge moment of the member 
cross-section. If the section is unsatisfactory, a new satisfactory 
section is selected for that member. 

Table 5.2.1 of B/20 specifies the factors as 1.2 and X^ as 
1,6, This is equivalent to X^ = 1,29 and X^ = 1.75 used in this 
thesis if the minimum yield stress is used instead of the design 
strength, B/20 clause 5.7,1 specifies the design strength as 
follows: 

"P may be taken as 0o93Y^ ..,, where is the minimum 

yield stress specified in BS 4360". 

In this design method the minimum yield stress is used as the design 
strength to account for the increased in the design strength up to 
the yield stress. The values of X^ and X 2 are therefore taken as 
1.29 and 1.75 respectively. The reason for taking these values is 
that if M is the total plastic moment of the section, and is 



the plastic modulus, the design strength, is given by =; 1,2 
Z or alternatively 0,93 y = 1,2 M /Z , This gives y =1,29 

■IT P P ' ^ 

Mp/Zp, For similar reasons a load factor of 1,75, taken in this 
thesis, is equivalent to 1,6 specified in Table 5,2,1 of B/20, 

It should be noted that the design program has been written 
in general terms, and the permissible load factors may therefore 
be assigned any values. However, for the examples iji this thesis, 
the values given above are used. 

Notes: 

(1) Sections must be selected for columns so that reverse 
column taper does not occur. If it does, then the beam section 
of the storey is increased to the next larger size until the same 
column section could be used for two successive storeys. It should 
be noted that a redesign is needed for each increase in a beam 
section, 

(2) Certain sections which are not suitable for plastic 
actions are omitted from the sections list. This is because fail¬ 
ure due to local buckling for these sections might occur before 

they reach their full plastic moments, 

$ 

(3) This thesis only deals with plane rigidly jointed frames 
and assumes that satisfactory restraints are available to prevent 
buckling out of the plane of the frame, 

3„7 STRENGTH REQUIREMENTS 

Once a set of sections are selected which satisfy the sway 
requirements, the member forces are checked and the task becomes 
that of satisfying the strength requirements if these are violated, 

3 ^ 7^1 Forces in a Member 


The axial force, shear force and bending moment of a member 





are obtained from the equation 


{p} = k A {x} C3.4) 

where {p}, A and {x} are as defined before in Section 3o2„ For 
a member linking joint i and j the member forces are calculated 


from: 


first end 


-al 

P 

-am 

P 

0 

al 

P 

am 

P 

b(t)_m 

5 p 

-b(j. 1 

5 p 

d4.2 

-bd)_m 

5 p 

b<^3l 

‘'♦2”p 

-d<)) 1 

2 p 

e+3 



d*2l»p 

-d<)).l 

2 p 

*♦4 




second end 


(3,5) 


Notice that joint i is the first end of the member, and joint j 

is its second end„ For a beam member with axes as shown in Figure 

211 is equal to 1 and m is 0„ If the axial stiffness, EA/L, 

° P P 

is ignored, then the first row of equation (3,5) should be dis— 

regarded. Furthermore, if the axial load effect is not considered, 

then <?,<!)><(>, <(> 3 .nd <()_ should be equal to 1. Thus the equa- 
1^ 2 3 4 b 

tion which calculates the member force in a beam is given by: 

First end Second end 


-b d 

-d e 

-d f 


(3,6) 


For a column member with axes as shown in Figure 2,1, 1^ = 0 and 

m = -1 Because the axial load is estimated as in Section 3,2, 
'■p ” 

the first row of equations (3,5) is ignored. Hence the column 


forces are given by: 
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The joint displacements x^, and 0^ at the first end of the 
member, and x., y., 9. at the second end are calculated by the 
iteration method given earlier. This also gives the second moment 
of area of the columns. These are obtained under unfactored load. 
Thus, the shear force and bending moments at both ends of a 
member are found for working load conditions. To obtain the inter¬ 
nal forces under factorized loading, the values of these forces are 


multiplied by the appropriate factor, 

3o 7,2 Checking the Beam Sections 

It was stated in Chapter 2, section 2.5 that under vertical 
loading each beam is prevented from collapsing by a beam mechanism 
at a load factor below X 2 » ^h^^ restriction on beam sections is 

applied when selecting the beam section even before considering the 
deflection requirements. 

Later at the end of the design cycles under vertical and wind 
loading, the bending moments at the ends of a beam and at mid-span 
are calculated from equation (3.6) under specified loads. For a 
unit load factor none of these is allowed to exceed the full plastic 
moment of the section as selected for the reduced cost design. If 
it is, then a larger section is selected, for the beam. This is 


done using; 


Z = M ,^^/Y 
p max s 


wliGTG M i s tn G maxiniuiTi 
wnexG 


bending moment in the member, Y is the 



mininiiim yield stress, and Z is th.e plastic modulus of the sectibno 

P 

If we consider the deflection requirements, an increase in the 
beam section will reduce the column sectionSa For this reason, if 
a beam section is changed, the iterative technique is required to 
be applied again to obtain new column sectionso These may turn 
out to be less than the previous ones„ The joint displacements 
will also change. These new values are used to calculate new bending 
moments. This process is repeated until no change in beam sections 
is required. 

3,7,3 Checking the Column Section for Strength 


At the end of the design cycles, the bending moment at each 
end of a column is calculated from equation C3o7) at a unit load 
factor for combined loading. The member forces including the 
axial forces are then increased by the load factor which is 
taken as 1.29, Tlie factorized maximum bending moment ^ is 

then compared with the reduced plastic moment of the column. If 
the reduced plastic moment is less than the in the column, 

then a larger column section is needed. The reduced plastic modu¬ 
lus of this is given by: 


7* = X M /Y 
■^p '^l max'-'s 


(3,8) 


Anew universal column is then selected from the tables so that it 
has a reduced plastic section modulus larger than Z^, The iteration 
process is then repeated but this time the object is to calculate 
the actual sway in each storey after the column sections are modified. 
Naturally, because the strength requirements dominate, the sway in 
the storeys will be less than those permitted. The process of 
altering the column section and finding the internal forces is 
repeated until no change in any column section is required. 
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CHAPTER FOUR 
COMPUTER PROGRAM 

4a INTRODUCTION 

This chapter presents a computer program for the design of two 
dimensional, rectangular, steel, sway-frames„ The program makes 
use of the design procedure described in the previous chaters. It 
is written in Fortran and runs on the ICL 1905 computer at the 
University of Aston, 

The program consists of a master segment calling a number of 
subroutines, the functions of which are illustrated by the flow 
diagram shown in Figures 5,1, 5,2 and 5,3, There are nine subrou¬ 
tines in the program. These are: 

I - Subroutine ITERATE: This constructs the modified stiffness 
equations and solves them, 

II - Subroutine REDUCE: This calculates the cost of the material 
for each storey of a frame. It finds the best set of beams which 
should be used in the next design cycle. It also decided the ratio 
for the second moment of area of an itemal column to that for an 
external column, 

III - The Subroutine INTERNAL FORCES: This calculates the bending 
moments and shear forces in each member of the frame from the joint 
displacements obtained by subroutine ITERATE, 

IV - Subroutine COLUMN-I: This chooses a section for a column, 
from the universal column sections, after the value of I is cal¬ 
culated by the iteration process, 

V - Subroutine BEAM-I: This selects a section for a beam from 
the universal beam sections. 
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yi - Subroutine COLUIIN-'ZP; This chooses a section for a column 
from the universal column sections, once the plastic modulus of 
the column is knowuo 

VII - Subroutine BEAM-ZP: This chooses a section for a beam 
from the universal beam sections, once the plastic modulus of the 
beam is known„ 

VIII - Subroutine REDUCED-ZP; Calculates the reduced plastic 
modulus for a coliamno 

IX - The subroutine STABILITY-F: Calculates the stability 
functions of a column membero 

The flowchart of the master segment shown in Figures Sol, 

5o2 and 5o3 is divided into 16 steps. The number of each step 
appears in brackets to the right of the chart. These steps are 
explained in some detail in the next three sections of this 
chapter. 

4.2 INPUT DATA (Steps 1 and 2) 

STEP 1: The method of numbering the joints was described in 
section 2.3. The program is written to expect at least one joint 
in each beam, whether there is a load applied on the span or not. 
The data consists of four groups these are: 
a - preliminary data - one card 
b - Beam-joints data - one card per joint 
c - Column-joints data - two cards per joint 
d - Frame data - two cards. 

An example of the data preparation is given in Appendix A of this 
thesis. 

The data preparation can be tedious especially for large 


frames, and might involve human errors 


For these reasons, a 
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PRINT RESULTS OF FINAL DESIGN BEFORE AND AFTER 
SELECTING COLUMN SECTIONS FROM THE LIST 


CALL ITERATE TO CALCULATE THE SWAY AFTER SELECTING 
COLUMN SECTIONS 


PRINT THE SWAY AND OTHER 
DEFLECTIONS 


MASTER PROGRAM [Continued) 








reduction in th,e data ia always useful. For regular frames a 
program i? written to prepare data for the frame. The program 
creates a data file for the original design program. The only 
data needed in this case is the preliminary data, the frame data, 
the height of the storey, and the loads applied on the frame, 

STEP 2: The properties of the universal columns and beams 
must be given to the program as data. To avoid punching these 
for each frame designed, they are kept in a data file on the back¬ 
ing store. This file is called each time a frame is to be designed. 
The properties kept in this file are the sizes of the section, the 
current price, the area, the second moment of area, the weight, 
and the plastic modulus of the section, 

4,3 DESIGN FOR DEFLECTION REQUIREMENTS (Steps 5 to 11) 

STEP 3: -A set of lower-bound sections is selected for the 
beams. To avoid collapse under vertical loading at a load factor 
^2 (1,75), the lower-bound beam sections are selected so as to 
avoid collapse by a rigid-plastic beam type mechanism. The sec¬ 
tions obtained in this manner are retained, and each time a beam 
is redesigned, the section is compared with its lower-bound beam 

section, 

STEP 4 1 The axial loads are calculated approximately as 
described in section 3,3, 

STEP 5: Preliminary second moments of area for beam sections 
are calculated according to section 2,5, either by using the method 
of Anderson and Islame, or by using the portal method, 

STEP 6; The simple portal method of design is applied to 
calculate the preliminary ratios of internal to external column 
second moments of area. Equations (2,34) give these ratios using 
the portal frame methodo 
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STEP 7; For each, beam section the second'’mbment of *axea is 
calculated in Step 5, a universal beam is chosen from the list by 
calling subroutine BEAM>'I„ 

STEP 8: The second moments of area of the columns, and the 
joint displacements are calculated for the first set of beams. The 
subroutine ITERATE for doing this will be discussed later in this 
chapter. The second moments of area of the columns resulting from 
this subroutine might form a reverse column taper. In this case 
the appropriate beam section is increased to the next larger size 
on the list. Subroutine ITERATE is called again to calculate the 
new second moments of area of the columns and the joint displace¬ 
ments, This process is repeated until no reverse column taper 
occurs, 

STEP 9: The cost of the material and the number of iterations 

needed for a design cycle are printed, 

STEP 10: The beam sections are reduced to the next smaller 
size on the list, and the procedure of steps 8 and 9 is repeated 
for the second design cycle, 

STEP 11: The results of the two most recent design cycles 
are kept in the memory to be used in the extrapolation for reducing 
the cost of the material of the frame. This is done by subroutine 
REDUCE, which calculates the cost of the material for each storey 
for the most recent two designs. The subroutine reduces the 
sections of the beams to the next smaller sizes in the list of 
universal beam sections, finds the cost of the material, and decides 
which beam section is most suitable to use in the next design 
cycle. It also decides the best internal to external column ratio. 
This process of reducing the material cost is repeated until any 
chance in a beam section results in a higher cost, or in a beam 







. 

• \ 

section sjualler than that s.pecified as a lower bound in step 3, 

The sections of the final reduced cost design and their properties 
are printed,, Also the joint displacements are printed in addition 
to the cost of the material of the frame, the total weight, and 
the number of iterations. 

4,4 DESIGN FOR STRENGTH REQUIREMENTS (Steps 12 to 17) 

STEP 12: The final reduced cost design should be checked for 
strength. To do so, the joint displacements, obtained for the last 
design cycle are used in calculating the internal forces in the 
frame. This is done by subroutine INTERNAL FORCES, The resulting 
bending moment at both ends of a member, the shear force, arid the 
approximate axial force are printed in IS units, 

STEP 13: The procedure of paragraph 3,5,2 is applied to 
check the beam sections. If the strength requirements in one of 
the beams are not satisfied, then that beam section must be incre¬ 
ased. A new plastic modulus for the section is obtained by using 

7 = M /Y , A universal beam section, in turn, is obtained by 

p max s 

calling subroutine BEAM-ZP, An increase in a beam section may 
result in smaller column sections. For this reason, subroutine 
ITERATE is called to obtain new column sections and joint dis- 
placements o 

STEP 14: The procedure of paragraph 3,5,3 is applied to check 

the strength of the column sections. The reduced plastic modulus 

of a column section should be more than that required by = 

i M /Y The reduced plastic modulus is obtained by calling 
1 max s° 

subroutine REDUCED- ZP, If a column section is not satisfactory, 
then the section is increased and a new universal column section 
is found by calling subroutine COLUHK-ZP. This increase in the 
colum section results in a decrease of the s»ay below the 






allowable. For this reason §ubroutint ITERATE is called. This 
time to obtain the sway and the rest of the joint displacements. 

The procedure described in steps 13 to 14 above is repeated 
until all the sections satisfy the strength requirements, 

STEP 15; The results of the final design which satisfy 
deflection, strength, and reduced cost requirements are printed, 
before and after selecting universal sections for the columns. 

STEP 16: Selecting a universal section for a column reduces 
the sway in a storey below the allowable. To examine the new 
sway deflections, the subrouting ITERATE is called to calculate 
the resulting sway and other joint displacements. These are due 
to the difference between the actual second moment of area of a 
column and that of the selected universal section, 

4,5 SUBROUTINE ITERATE 

subroutine is used to calculate the second moments of area 
of the columns if the sway is specified, or to calculate the joint 
deflections if the second moments of area of the columns are known. 
In both cases the vertical displacements and the rotations of the 
joints are also calculated. A comprehensive flowchart for this 
subroutine is given in Figures 5,5, 5.6, 5,7 and 5,8, in which the 
integer NIR is used to count the number of iterations as they are 
carried out. JB and JC are used to count the joints of beams and 
columns as they are handled. 

The subroutine starts by setting the values of the joint dis¬ 
placements X, y and R, and the second moment of area of the column 
I to their values in the previous design. Obviously for the first 
design cycle these are set to their preliminary values (i,e, I = 1, 
V = -1 and R = 0), Level 600 of Figure 5,5 in the chart marks the 
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FIGURE 5o8: SUBROUTINE ITERATE (continued] 
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start of one iteration. The operations, between this leyel and 
level 300 in Figure 5,8 represents one iteration. The process of 
calculation for an iteration is divided into nine steps. The number 
of each step appears in brackets at the right of the chart, 

STEP 1: The number of iteration NIR is set to zero and is 
increased by one for each iteration, 

STEP 2: For each beam joint the contributions of the beams 
connected to the joint are calculated. Then the vertical displace¬ 
ment of the joint y and the rotation R „ are found, using 

J B JB 

equations (2,14) and (2.17) respectively, 

STEP 3: For each column joint the (j) functions are calculated 

for columns C. and C connected to joint JC, However, for the 
1 m 

first design cycle the stability functions are set to 1, as the 
axial load effect is neglected in this cycle, 

STEP 4: The second moment of area of a column needs to be 
calculated once in each storey. The rest of the column second 
moments of area are then found by proportion, Tlie calculation is 
carried out for the extreme left joint of the storey, and the value 
of I is obtained from equation (2,38), However, if the second 
moment of area of the column is known and the sway is wanted, 
then equation (2,29) is used instead to calculate the sway for the 
complete storey, 

STEP S* The second moments of area of the other columns in 
a storey are found from the internal to external columns ratios by 
multiplying I by the appropriate ratio p„ 

STEP 6‘ The vertical displacement of joint JC and the joint 
rotation must be found for each column joint. This is done by 
solving equations (2,23) and (2,26) respectively, 

STEP 7: To determine whether a desired degree of accuracy (Z) 
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is acquired or not. The vertical displacement y and the joint 
rotation R are examined using the inequality 

The current value - The Previous Value ^ ^ 

The Previous Value 

where Z is the required accuracy„ If the above condition is not 
satisfied, then another iteration will be needed starting from level 

600 o 

STEP 8: For each column joint, the above inequality is used 
to establish whether the value of I or x achieves the required 
accuracy or not« 

STEP 9: For each column joint the test is also applied on the 
value of the vertical displacement y and the rotation of the joint 

Ro 

In both steps 8 and 9 the control is transferred to level 
600 to start a new iteration if the above inequality is not satis¬ 


fied. 


CHAPTER FIVE 


DESIGN EXAMPLES 

5„1 A 12 STOREY-SINGLE BAY FRAME 

Most of the existing design methods for deflection limitation 
do not consider the axial stiffness and the axial load effect (see 
Chapter One). In order to show these effects, the 12 storey single 
bay frame of Figure 5.1 was first designed disregarding the axial 
stiffness EA/L and the stability functions. The deflection pro¬ 
file used in the design is shown in Figure 5.2. The value of k 
used in the deflection curve was 30. 

The second column in Table 5.1 shows the second moments of 
area of the initial beam sections. These are obtained as follows: 

1 - The beam sections are designed using the simple portal frame 

method of the steel designer's manual (Allwood, et al, 1966). 

It should be noted that no beam section is allowed to be 
below that obtained by a beam mechanism. 

2 - The second moment of area of the beam sections obtained in 

step (1) are magnified by a factor. This is taken as 2 at 
the top storey, and 4 at the bottom. From this initial 
design, eleven cycles are needed to obtain the final linear 
reduced cost design. 

Column three of Table 5.1 shows the second moments of area 
for this design after sections have been selected from the list. 

The number of iterations per design cycle is represented in Figure 
5.3. This shows the reduction in the number of iterations after 
the first design cycle. Figure 5.4 is for the reduction in material 
cost in each design cycle. Trial 6 shows an increase in the 





material cost, due to the increase in some,beam sections to prevent 
reverse column taper. The computer time needed to complete the 
design, with a tolerance of 0,01, was 27 second on Aston 1904s 
Computer (at a cost of £2,37), This is less than the 29 second 
needed for one explicit analysis solving L = ^ )C directly. 

In the second stage, the beams for the final linear design 
were kept unaltered, but the columns were subsequently altered to 
allow for the effect of axial forces. These v\fere calculated 
exactly using equations (3,3), with a tolerance of 0,01, The 
design with the axial load effect included is shown in the fourth 
column in Table 5,1 This gives the second moments of area of the 
sections selected from a list for universal beams and columns. 

The time needed to obtain the final non-linear design was 82 
seconds. This includes the time to check the strength of the 
final non-linear design, Tlius, when considering the axial load 
effect the computer time increases considerably, but even then the 
computer cost of a complete design is only £7,20. 

The effect of the axial stiffness and the stability functions 
on the design of the column sections is obvious, A glance at the 
third and fourth columns in Table 5,1 shows that, except for the 
top two columns, all the sections are increased when the effect of 
axial load is considered. This will cause the total material cost 
of the frame to increase by 12,24% from £5608,18 to £6390,69, 

The total mass increases by 11,18% from 25832,28 kg to 29082,36 kg. 
To check the safety of the final design, a non-linear elastic- 
plastic analysis was carried out for this frame by an existing 
computer program (Majid and Anderson, 1967), This showed that 
under combined loads the frame remained elastic at a load factor of 
1,29, No plastic hinge occurred in any member of the frame at this 


load factor. 



5.2 A 24 STOREY FRAME WITH THREE EQUAL BAYS 

The computer program described in Chapter Four was used to 
design the 24 storey frame shown in Figure 5„5„ This frame belongs 
to a structure consisting of parallel frames which are 4„5 m apart. 
Th'e vertical and the wind loads are calculated according to CP3 
Parts 1 and 2„ The wind load is calculated for a frame to be 
built on a surface with large and frequent obstructions. Each 
frame is symmetrical and is loaded symmetrically. 

The deflection curve used in the design and the maximum allow¬ 
able sway in each storey are shown in Figure 5,6, The value of k 
used in the deflection curve is 26, The results of the design in 
which the axial load effect is included are shown in Table 5,2, 

This gives the second moments of area of the selected sections, in 
which steel grade 50 is used. 

The design was carried out with a tolerance of 0,001, How¬ 
ever, the same frame was designed with a tolerance of 0,01 and the 
same sections were obtained with considerable saving in computer 
time. The time needed for designing the frame with a tolerance 
of 0,001 is 1020 seconds, which means that the cost of the design 
for a user from outside the university would be £88, The cost 
reduces to £23,76, if a tolerance of 0,01 is used. The computer 

time needed in this case is 267 seconds. 

To start with, the beam sections were taken from Anderson 
and Islam’s design described in Chapter One, using the linear 
deflection curve, Tne sway between the two ends of a column in 
their design was taken to be H/300, where H is the column's length. 
The results of the design using their method is shown in Table 
5.3, which gives the second moments of area of the selected sections 





Anderson and Islam's design does not consider ^the effect of the 
vertical load on the sway. Therefore, to compare the results, a 
third design should be obtained. This is done by using the 
computer program described in Chapter Four, but neglecting the 
axial stiffness of the columns and the axial load effects. Table 
5,4 shows the second moments of area for the selected sections. 

The price of the material and the total mass of the frames 
are given in Table 5,5 for the three designs. It can be seen that 
the design obtained by this method (neglecting the axial stiffness 
and the axial load effect) is 1,22% cheaper than that obtained by 
Anderson and Islam's method. It is also noticed that if the axial 
load effects are included in the design of this frame then the 
material cost will increase by 11,45%, 

In order to justify this increase and to examine the efficiency 
of Anderson and Islam's design, each frame was analysed using an 
existing computer program. The differential sway between the ends 
of each column, obtained by each analysis, is shown in Figure 5,7, 
The vertical line shows the h/300 limit on the differential sway. 
Curve 1 shows sway in each column for the final design. None of 
the points of this curve is to the right of the vertical limiting 
line. Curves 2 and 3 show the differential sway between the ends 
of each column for the frame designed by Anderson and Islam's 
method, given in Table 5,3, Curve 2 is for wind load only, exclu¬ 
ding axial load effect. Curve 3 is for wind and vertical loads, 
including axial load effects. Because this frame is symmetrical 
and loaded symmetrically, the difference between curves 2 and 3 
is solely due to the effect of axial loads. It can be seen that 
is si^[nificant and cannot be neglected. In fact 
an increase of 11,45% in the material cost is needed if the axial 




load effects are considered in the final design. 

Twelve design trials are needed to reduce the cost of the 
material from the initial beams, obtained from Anderson and Islam's 
design. The reduction in the material cost in each trial is shown 
in Figure 5,8, Trials 8 and 9 show an increase in the material 
cost, due to the increase in some beam sections to prevent reverse 
colvunn taper. The number of iterations used in each design trial 
for a tolerance of 0,001 is plotted in Figure 5,9, This shows a 
considerable reduction in the number of iterations after the second 
trial. This is due to the use of the previous design results as 
initial values for the current ones. The relatively large number 
of iterations in the second design trial is due to the introduction 
of the axial load effect at this stage. 

The final design shown in Table 5,2 was obtained by applying 
the deflection requirements only. However, the design was checked 
for strength requirements as described in Chapter Three sections 
(3,5, 3.5,1, 3. 5.2 and 3,5,3). This showed that all the sections 
were satisfactory. Furthermore, a non-linear elastic-plastic 
analysis of the frame with the final sections was carried out using 
an existing computer program (Majid and Anderson, 1967), This 
analysis showed that under combined loads the frame remained elastic 
at a load factor of 1,29. It should be noted that, for strength 
requirements a beam was allowed to develop a plastic hinge below 
a load factor of 1,29, but because of the deflection requirements 
governing the design, it was found that all the beams remained 
elastic at this load factor 
5^3 A 35 STOREY - TtVO UNEQUAL BAY FRAME 

The 35 storey - two unequal bay frame shown in Figure 5,10 
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was designed from steel grade 50. In the out-of-plane direction, 
the frames were assumed to be 4.5 m apart. The vertical and the 
wind loads were calculated according to CP3 Parts 1 and 2„ The 
wind load was calculated from a frame to be built on a surface 
with large and frequent obstruction. Because the frame is assym- 
metrical, both vertical and wind loads affect the deflection. The 
linear and the non-linear profiles of the frame are shown in 
Figure 5.11. For practical reasons, the left and the right spans 

in each storey are to be built out of the same beam section. 

The frame was first designed using Anderson and Islam's 

equations, Tne sway between the ends of the column was taken as 
height/300, The results of the design using their equations with 
the linear deflection curve are shown in Table 5.6. This 

shows the second moment of area of the members after selecting 
sections from the list of universal beams and columns. The 
design represented in the table was analysed using an existing 
computer program (Celik, 1977), Four different load cases were 
considered, and the sway between the ends of each column is 
shovvTi in Figures 5,12 and 5,13, 

(1) Curve 1 of Figure 5,12 is for wind load only with the wind 
acting from the leftj and neglecting the axinl load effects, 

(2) Curve 2 in Figure 5.12 is for the combined loading, with the 
wind acting from left; and considering the axial load effects, 

(3) Curve 1 of Figure 5.13 is for wind load only, with the wind 
acting from the right; and neglecting the axial load effects, 

(4) Curve 2 of Figure 5,13 is for combined loading, with the 
wind acting from the right; and considering the axial load effects. 
In both Figure 5,12 and 5,13 the vertical line shows the permiss¬ 
ible limit on sway, A point to the left of this vertical line 


means that the sway at a storey is acceptable, while a point to 
the right means it is not. It is obvious from Curves 1 and 2 in 
both figures that the sway exceed, the limit of H/300 by a con¬ 
siderable margin. In fact for the combined load case with the 
wind acting from the right, the sway of the 33rd floor is equal 
to H/120. This is 2,5 times larger than the limit. 

The design obtained from Anderson and Islam's equations was 
therefore found to be totally unsatisfactory. The beams obtained 
from their equations with the sway restriction of H/300 cannot be 
used as a starting point for the iteration technique. However, 
Anderson and Islam's approach can still be utilized to select 
trial initial beams provided that the allowable sway is severely 
restricted. It must be emphasized that this restriction on sway 
only applies when selecting beam sections from Anderson and Islam's 
equations. But when using the design method proposed by the author, 
the code limitation of H/300 is observed. In other words Anderson 
and Islam's approach is only used to select an initial trial set 
of sections. 

In the design procedure described in Chapter Two, economy in 
the cost of the material of the frame was achieved by reducing each 
beam section to the next available in the list of universal sections. 
For large frames such a policy can increase the number of design 
cycles. The design procedure can be speeded up by drawing a curve 
relating the cost to the beam sections as follows: 

(i) Anderson and Islam’s equation are used to calculate several 
sets of beam sections under different restrictions of the sway 
between the two ends of the column. Some of these, however, may 
be too weak to resist failure by a beam type mechanism. To avoid 
such a failure, these initial trial sections are Increased 


substantially. 


(ii) For each initial set of beams thus obtained, the columns are 
designed using the method developed in this thesis to satisfy the 
deflection profile given by curve b of Figure 5„11„ The axial load 
effects are included when designing these columns, A graph is 
drawn to show the cost of material for different designs as obtained 
from a number of trial sets of beams. Figure 5.14 shows the cost/ 
beam-section curve for the 35 storey frame. It can be seen that 
six trials are needed to draw the curve. This is considerably 
less than the 12 trials, which was needed to obtain the reduced 
cost design of the 24 storey frame shown in Figure 5,5 Figure 5,14 
shows that the most economic design is obtained when the sway in 
Anderson and Islam's equations are restricted to H/800, A sub¬ 
sequent analysis of this design showed that plastic hinges deve¬ 
loped in some of the columns below a load factor of 1,29, This 
is perfectly acceptable provided that failure takes place above 
such a load factor. However to avoid the formation of such 
hinges two modifications can be made: 

(1) The column sections where hinges occur could be increased, 
until an analysis shows that no hinge is developing below a load 
factor of 1.29, The beam sections are kept unchanged. 

(2) A smaller set of beams could be used to satisfy sway require¬ 
ments so as to obtain larger columns from the iteration. 

In fact the second solution was found to be more economical. 

In order to produce a design which satisfies both sway and 
strength requirements, the second approach described above was 
therefore used. First, the iteration technique using a set of 
beams obtained from Anderson and Islam's equations with a restric¬ 
tion of H/700 was analysed. This was found to be similar to that 



of H/800, Then, the design using a set of beams with sway restribS- 
tions of H/600 was analysed. The beams in this case were relatively 
small and the resulting columns were relatively large. The analysis 
of the frame designed in this manner under combined loading 
showed that no hinge developed below a load factor of 1,29, and the 
frame was completely elastic at this load factor. Although hinges 
were allowed to develop in the beams below a load factor of 1,29 
(for combined loading), it was found that these too remained 
elastic. 

The designs of the frame obtained using the computer program 
described in Chapter Four are summarized in Table 5,7 and 5,8, 

These show the second moment of area for the selected sections from 
combined loading, with the wind acting from either side. The 
second moment of area obtained by using Anderson and Islam's 
equations, with the restriction of H/600 on sway, plus the second 
moment of area for a section that satisfies a beam mechanism provides 
the actual beam section used in the iteration technique. The 
deflection profile used is that of Figure 5,11 (the non-linear 
curve b), The design was carried out with a tolerance of 0,001, 
and when a tolerance of 0,01 was used a premature convergence was 
noticed. It is obvious from Table 5,7 and 5,8 that the critical 
design is that for wind acting from the right. This design is 
considered as the final design because the column sections obtained 
are larger than those resulting from the wind acting from the left. 

It should be noticed that in this example, only two design 
cycles were sufficient for each initial set of beams. In the first 
design cycle the axial load effect and the axial stiffness were 
neglected. In the second design cycle the results of the first 
cycle were used to introduce these effects. 




The number of iterations needed for the first design cycle is 
565, while 493 iterations are needed for the second design cycle. 
Comparing these with the number of iterations needed for other 
frames in the first two design cycles shows that the number of 
iterations increases rapidly for large frames. However, the computer 
time needed, to obtain the design of Table 5.7 , 1076 second. 

This includes the time to check the strength of the final design, 
which means that the cost of the design fop a user from outside 
the university would be £92,05, The store needed was 3120 words 
(approximately 32 k), This included storing the computer program 
itself. 

To ensure that the sway is under the limit and to examine the 
accuracy of the program, the design of Table 5,7 was analysed 
considering the axial load effects. The sway between the ends of 
each column is shown by curve 3 of Figure 5,12, for combined 
loading with the wind acting from the left. This shows that all 
the points are within the permissable limit on sway represented 
by the vertical line. For the case of wind acting from the right, 
the sway of each storey is shown by curve 3 of Figure 5,13, It 
can be seen that for the case of combined loading with wind 
acting from the right, some of the point are slightly out of the 
limit. The reasons for this are; 

(1) The approximation made in calculating the axial loads in the 
iteration method, 

(7) The change of the second moment of area of the column due to 
a selection of a column section from the list of universal 

columns, 

5.4 CONCLUSIONS 


From the examples given above, it is possible to conclude 
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the following: 

1 - The examples indicate that the iteration technique can be 

used. Successfully for designing different sizes of frames. 
In all the frames designed, the iteration converged to a 
set of acceptable column sections. 

2 - For larger frames a small tolerance such as 0.001 is required 

to obtain a complete convergence. The 35 storey frame showed 
that a premature design was obtained when a tolerance of 
0,01 was used. However, the latter was used successfully 
for all the other frames. Obviously, a larger number of 
iterations and design cycles are needed for a larger frame. 

3 - A frame which is designed to satisfy deflection requirements 

results in the frame having an overall strength which is 
greater than that required by the elastic-plastic design 
methods. In fact, for all the above designs the non-linear 
elastic plastic analysis under combined loading shows that 
the frames remain elastic at a load factor of 1.29, 

4 - The cost of the computer time and storage is negligible 

compared to the material cost of the frame. 






















Initial 

Linear 

Non-Linear 

Level 

design 

I in cm'^ 

design 

I in cm^ 

design 

I in cm^ 


BEM4S 



Level 1 

503781 

324715 


Level 2 

435796 

204747 


Level 3 

375110 

245412 


Level 4 

375110 

168535 


Level 5 

245412 

168535 


Level 6 

245412 

135972 


Level 7 

168535 

117700 


Leve1 8 

135973 

75549 


Level 9 

75549 

55225 


Level 10 

75549 

47363 


Level 11-12 

75549 

21345 


324715 

204747 

245412 

168535 

168535 

135972 

117700 

75549 

55225 

47363 

21345 


COLUllNS 


Level 1-2 
Level 3 
Level 4-5 
Level 6 
Level 7-8 
Level 9 
Level 10 
Level 11 
Level 12 


66307 

57153 

48525 

38740 

32838 

27601 

17510 

11360 

5263 


79iiw 

79110 

66307 

57153 

57153 

57153 

32838 

1751C 

5263 


Total Mass 


25832,28 kg 29082,36 kg 


Total cost of material 


15608,18 


16390,69 


TABLE 5,1: DESIGNS OF THE 12 STOREY FRAME 
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Second moments of area x lo'^ 

Storey 

Beam 

External 

Internal 

Level 

4 

mm 

Column iran^ 

Column mm^ 

Level 1 

87260 

-i 

99994 

183118 

Level 2 

87260 

66307 

122474 

Level 3 

87260 

50832 

99994 

Level 4 

87260 

38740 

79110 

Level 5 

87260 

32838 

66307 

Level 6 

87260 

32838 

66307 

Level 7 

87260 

32838 

66307 

Level 8 

87260 

32838 

50832 

Level 9 

87260 

22416 

50832 

Level 10 

87260 

22416 

50832 

Level 11 

66610 

22416 

50832 

Level 12 

66610 

22416 

50832 

Level 13 

55225 

22416 

50832 

Level 14 

55225 

22416 

50832 

Level 15 

40956 

22416 

50832 

Level 16 

40956 

22416 

50832 

Level 17 

33324 

22416 

38740 

Level 18 

29337 

17510 

38740 

Level 19 

21345 

17510 

32838 

Level 20 

21345 

14307 

32838 

Level 21 

15603 

9462 

17510 

Level 22 

10054 

6088 

14307 

Level 23 and 24 

10054 

6088 

6088 

Total mass 

111597,0 kg 

Total cost 

£23817,19 


TABLE 5,2: DESIGN OF THE 24 STOREY - 3 EQUAL BAY FRAFE 




Second moments of area x 10'^ 


Storey 

Level 

Beam 

mm^ 

External 
Column mm'^ 

Internal 
Column mm'^ 

Leve1 1 

87260 

38740 

79110 

Level 2 

87260 

38740 

79110 

Level 3 

87260 

38740 

79110 

Level 4 

87260 

38740 

79110 

Leve1 5 

66610 

38740 

66307 

Level 6 

66610 

32838 

66307 

Level 7 

61530 

32838 

66307 

Level 8 

61530 

32838 

66307 

Level 9 

55225 

32838 

66307 

Level 10 

55225 

32838 

66307 

Leve1 11 

55225 

32838 

50832 

Level 12 

55225 

32838 

50832 

Level 13 

40956 

22416 

50832 

Level 14 

40956 

22416 

50832 

Level 15 i 

1 

37039 

22416 

38740 

i 

Level 16 

33324 

17510 

32838 

Level 17 

29337 

14307 

32838 

Level 18 

25464 

14307 

32838 

Leve1 19 

21345 

14307 

22416 

Level 20 

15603 

9462 

17510 

Leve121' 

15603 

7647 

14307 

Level 22 

8167 

6088 

14307 

Level 23 and 24 

4381 

6088 

6088 

Total mass 

100398o0 kg 

Total cost 

£21347ol8 


ANDERSON AND ISLAM’S DESIGN OF THE 24 STOREY 
- 3 EQUAL BAY FRAME 


TABLE 5o3: 









Second moments of area x 10"^ 

Storey 

Beam 

External 

Internal 

Level 

mm^ 

Column mm^ 

Column mm^ 

Level 1 

87260 

79110 

146765 

Level 2 

98408 

38740 

79110 

Level 3 

61530 

38740 

79110 

Level 4 

87260 

32838 

66307 

Level 5 

66610 

32838 

50832 

Level 6 

87260 

32838 

50832 

Level 7 

87260 

22416 

38740 

Level 8 

87260 

17510 

32838 

Level 9 

87260 

14307 

32838 

Level 10 

87260 

14307 

32838 

Level 11 

66610 

14307 

32838 

Level 12 

55225 

14307 

32838 

Level 13 

55225 

14307 

32838 

Level 14 

66610 

14307 

22416 

Level 15 

1 

37039 

14307 

22416 

Level 16 

55225 

9462 

22416 

Level 17 

29337 

7647 

17510 

Level 18 

37039 

7647 

17510 

Level 19 

21345 

7647 

14307 

Level 20 

29337 

7647 

9462 

Level 21 

15603 

6088 

6088 

Level 22 

15603 

6088 

6088 

Level 23 and 24 

10054 

6088 

6088 

Total mass 

99080„0 kg 

Total cost 

£21088„36 


TABLE 5„4: DESIGN OF 24 STOREY - 3 BAY FRAI4E NEGLECTING 
AXIAL STIFFNESS .AND AXIAL LOAD EFFECT 







Type of Design 

Cost of the 
material 

Total mass 
of the frame 
in Kg 

The final design includes 

the axial stiffness and 

the axial load effect 

23817,19 

111597,0 

Design neglecting axial 

stiffness and axial load 

1 

effect 

21088,36 

99080.0 

Anderson and Islam's Design 

21347,18 

100398.0 


TABLE 5.5: COST AND MASS OF THE 24 STOREY 
3 EQUAL BAY FRAME DESIGNED BY VARIOUS METHODS 






























Second moments of area for selected sections 
in oclO^ mm^ 

Storey 

Left 

Right 

Left 

Internal 

Right 

Level 

Beam 

Beam 

Column 

Column 

Column 

Level 1-4 

278833 

87260 

122474 

183118 

66307 

Level 5-6 

278833 

87260 

99994 

183118 

66303 

Level 7 

278833 

66610 

99994 

146765 

66303 

Level 8 

278833 

66610 

99994 

146765 

50832 

Level 9 

204747 

61530 

99994 

146765 

50832 

Level 10-12 

204747 

55225 

99994 

146765 

50832 

Level 13 

204747 

55225 

79110 

122474 

50832 

Level 14-15 

168535 

55225 

79110 

122474 

50832 

Level 16 

168535 

55225 

79110 

122474 

38740 

Level 17 

168535 

55225 

66307 

122474 

38740 

Level 18 

168535 

55225 

66307 

99994 

38740 

Level 19-20 

135972 

40956 

66307 

99994 

32838 

Level 21 

135972 

37039 

66307 

99994 

32838 

Level 22-23 

111673 

33324 

50832 

79110 

32838 

Level 24 

98408 

29337 

50832 

66307 

32838 

Level 25 

1 87260 

29337 

38740 

66307 

22416 

Level 26 

87260 

25464 

38740 

66307 

22416 

Level 27 

87260 

21345 

32838 

50832 

17510 

Level 28 

1 61530 

21345 

32838 

50832 

17510 

Level 29 

1 55225 

21345 

32838 

38740 

14307 

Leve1 30 

1 55225 

15603 

22416 

32838 

14307 

Level 31 

! 37039 

j 

15603 

17510 

32838 

9462 

Level 32 

1 29337 

8167 

14307 

22416 

7647 

Level 33 

^ 21345 

5415 

9462 

17510 

6088 

Level 34-35 

10054 

i 

4381 

6088 

9462 

6088 

Total mass 138773^50 kg 

Total cost £30253o82 


TABLE 5.6: TRIAL INITIAL SECTIONS FOR A 35 STOREY - 2 


UNEQUAL BAY FRAME OBTAINED BY ANDERSON AND ISLAM’S 

EQUATIONS USING THE LINEAR DEFLECTION CURVE WITH a=l/300 
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Second moments of area of the selected sections 

storey 

Levels 

Beam 

mm4 X 104 

right 
column 
mm4 X 10^ 

Internal 
column 
mm4 X 104 

1 ef t 
column 
mm4 X 104 

Level 1 

503781 

146765 

378668 

275140 

Level 2-4 

503781 

99994 

227023 

146765 

Level 5-8 

503781 

79110 

227023 

146765 

Level 9-12 

435796 

79110 

227023 

146765 

Level 13-16 

375110 

79110 

227023 

146765 

Level 17-19 

324715 

79110 

227023 

146765 

Level 20-21 

278833 

79110 

227023 

146765 

Level 22 

278833 

66307 

227023 

146765 

Level 23 

278833 

50832 

227023 

146765 

Level 24 

278833 

50832 

183118 

122474 

Level 25 

204747 

38740 

146765 

99994 

Level 26 

204747 

32838 i 

146765 

99994 

Level 27 

168535 

32838 

122474 

79110 

Level 28 

168535 

32838 

99994 

66307 

Level 29 

135972 

22416 

99994 

66307 

Level 30 

135972 

17510 

79110 

50832 

Level 31 

98408 

14307 

50832 

38740 

Level 32 

87260 

9462 

50832 

32838 

Level 33 

61530 

6088 

32838 

22416 

Level 34 

55225 

6088 

14307 

7647 

Level 35 

25464 

6088 

6088 

6088 

Total mass 221 

865.17 kg 

Total cost 

£49275.87 



TABLE 5o7: DESIGN OF THE 35 STOREY - 2 UNEQUAL 
BAY FRAME, COMBINED LOADING, WIND FROM RIGHT 




Second moments of area of the selected sections 


Storey 

Level 

Beam 

mm4 X 10^ 

Level 1 

503781 

Level 2-8 

503781 

Level 9-12 

435796 

Level 13-16 

375111 

Level 17-19 

324715 

Level 20-21 

278833 

Level 22 

278833 

Level 23 

278833 

Level 24 

278833 

Level 25 j 

204747 

Level 26 

204747 

Level 27 

168535 

Level 28 

168535 

Level 29 

135972 

Level 30 

135972 

Level 31 

98408 

Level 32 

87260 - 

Level 33 

61530 

Level 34 

55225 

Level 35 

25464 ' 

_— 

Total mass 2 1 54 14.50 kg 


right 


Internal 


146765 
99994 
7 9110 
79110 
79110 
79110 
66307 
50832 
50832 
38740 
32838 
32838 
32838 
22416 
17510 
14307 
9462 
6088 
6088 
6088 


378668 

227023 

227023 

227023 

227023 

227023 

183118 

146765 

122474 

122474 

99994 

79307 

66307 

66307 

50832 

32838 

32838 

17510 

9462 

6088 


left 

colramn 


275140 

146765 

146765 

146765 

146765 

146765 

122474 

99994 

79110 

66307 

66307 

50832 

50832 

38740 

32838 

22416 

17510 

14307 

6088 

6088 


TABLE 5,8: 


DESIGN OF THE 35 STOREY - 2 UNEQUAL 


BAY FRAME FOR ^ C OM BINE D LO A D ING WIN D F ROM 
LEFT 
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CHAPTER SIX 

DESIGN OF REINFORCED CONCRETE FRAMES 

6ol INTRODUCTION 

The design procedure given in Chapters 2-5 for the design of 
steel frames is used in this chapter for the design of reinforced 
concrete frames. The CP110:1972 Code of Practice for The Struc¬ 
tural Use of Concrete is based on the principles of Limit State 
Design, The object of limit state design is to achieve an accept¬ 
able probability that a structure will not become unservicable in 
its lifetime. It, therefore, sets out to examine all the ways in 
which a structure may become unfit for use. 

The condition of a structure when it becomes unserviceable 
is called a "limit state". The most important of these limit 
states which must be examined when designing a frame are: 

(a) the ultimate limit state which requires that neither the 
whole structure nor any part of it should collapse under 
foreseable overload (strength requirements); 

(b) the serviceability limit state of deflection which requires 
that the deflection of the structure should not adversely 
affect its appearance and; 

(c) the serviceability limit state of cracking. 

The usual approach is to design a reinforced concrete frame 
on the basis of the most critical limit state and then check 
that the remaining limit state requirements are not violated. In 
general, engineers often try to satisfy the ultimate limit state 
^©quirement) first and then check that the deflections 
of same joints or the sway in the members within acceptable limits. 
They may also check that any crack that may develop is within 






acceptable limits» However, for sway frames, especially forjitall 
slender frames the sway in the columns can often be more critical 
than strengtho For this reason, it is believed that the rational 
approach for the design of such frames will be to consider the 
limit state of deflection first and then proceed to check for 
strength and crack development., 

Clause 2o2o3„l on deflection in CPllO recommends the following 

"The deflection of the structure or any part of the 
structure should not adversely affect the appearance or 
efficiency of the structure. 

The engineer must satisfy himself that deflections 
are not excessive having regard to the requirements of 
the particular structure „»,, 

The effect of lateral deflection should be con¬ 
sidered, particularly for tall slender structures 

The clause does not recommend any specific value as the limit 
of the horizontal deflection for multistorey frames. However, the 
original draft did consider a severe limit of height /500 (Allen, 
1977)o In the absence of a better limit, and as a large sway 
might produce excessive cracks, the limit of height/500 is adopted 
for the design of concrete frames in this thesis, 

6,2 THE PRINCIPLES OF DESIGN 

The design method described in Chapters 2-5 is also applied 
here to generate the values of the required second moments of area 
for the columns of a multistorey sway frame. This is done to 
satisfy the horizontal sway limitations imposed by the service¬ 
ability limit state of deflection. As before, this is conditional 
on that the values of the second moments of area for the beams 
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being known before starting the iteration. The second moment of 
area of each column is then used to obtain the minimum column 
depth required so that the sway does not exceed the limit. 

It was found, when designing several frames, that the actual 
column depth obtained by this method is considerably greater than 
the minimum imposed by the code. 

The code suggests in clause 3,5,8 that if the depth of a column 
is more than its effective height/30 then the sway requirement 
may be satisfied without further calculation. However, the author 
found that this underestimates the column depth needed tor multi¬ 
storey frames and particularly the depth of the columns at the 
lower storeys. 

The 1972, CPI10 also requires that the frame should be designed 
for the serviceability limit state of deflection for different load 
cases. To produce a design which satisfies all the load cases, 
the column sizes obtained, when designing for the first loading 
case, are used as lower bounds for the second loading case, and 
so on. The application of lower bounds obtained when designing 
for a previous load case ensures that each load condition is satis¬ 
fied. 

In the proposed design method, the sections obtained from 
the final deflection limitation design are checked for strength 
requirements. Again the frame has to be checked for all the 
ultimate limit state load cases. For this purpose, an elastic 
analysis using the iteration method given in Chapter 2 is carried 
out for each load case. Once again the axial load effect and the 
axial stiffness (EA/L) are taken into consideration. The maximum 
bending moments, the shear force, and the axial force from all the 
load cases are obtained, and the sections are checked for these 
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ultimate values, 

6,3 DESIGN ASSUMPTIONS 

When iteratively analysing a reinforced concrete frame for 
the ultimate limit state, the relative stiffness of the members 
may be based on any one of the following (CPllO, ,clause 2,4,3,1): 

(a) The concrete section: the entire concrete cross-section, 
ignoring the reinforcement, 

(b) The gross section: The entire cross-section, including the 
reinforcement on the basis of modular ratio. The second 
moment of area for this section is called the uncracked 
second moment of area as all the area of the concrete is 
considered. 

(c) The transformed section: The compression area of the con¬ 
crete cross-section combined with the reinforcement on the 
basis of modular ratio. The second moment of area for this 
section is called the cracked second moment of area as the 
tension area of the concrete cross-section is ignored. 

The work of Beeby and Taylor (1978) shows that the results 
of an analysis vary considerably with the assumption made in 
choosing the properties of the sections. Their analysis of a 
portal frame using different assumptions showed that the maximum 
bending moment was obtained when the cracked second moment of 
area for the beam section was used and the uncracked one for the 
column section. 

In a multistorey frames a column should be designed to resist 
the axial load and the bending moment. Usually the axial force in 
a column dominates the design of the section (except for the top 
one or two storeys where small axial forces may occur), If the 


moment is small and the axial load is comparatively large,tfthen the 
position of the neutral axis would be outside the section. In 
this case the column area would be under compression and the sec¬ 
tion uncracked. For the reasons described above, it is logical 
to assume that: 

i) The relative stiffness of a column section is based on 

assumption (b), and the uncracked second moment of area is 
used, 

ii) For a beam section the bending moment would dominate the 

design of the section and the moment would be large enough 
to make the major parts of the beam crack. Therefore, it is 
reasonable to assume that the relative stiffness of a beam 
section is based on assumption (c), and the cracked second 
moment of area of a rectangular section is used. 

Assumptions (i) and (ii) are used in this chapter for evalu¬ 
ating the relative stiffness for both ultimate limit state and 
serviceability limit state of deflection. In fact CP110:1972 
Appendix A, clause 1, states that greater accuracy may be achieved 
in the deflection calculation if the relative stiffness of the 
cracked beam is used. 

iii) The effects on deflection of temperature, creep,and shrink¬ 
age are ignored in this thesis, 

iv) The bar spacing rules should be satisfied so that no excessive 
cracking occurs. In this respect, it is assumed that the 
reinforcement will be distributed within the section so that 
the distance between two bars does not infringe the maximum 
and minimum limits. Furthermore, the cover for the rein¬ 
forcement is chosen according to the characteristic strength 
of the materials to prevent the influence of weather on 




reinforcement. Minimum concrete covers for the reinforce¬ 
ment are given in clause 3,11,2, and Table 19 of CPI10, In 
the design examples given in this chapter these recommendations 
are taken into consideration, 

v) For practical purposes the overall depth of a beam or a 
column is rounded to the nearest 50 mm, 

6,4 REINFORCEMENT REQUIREMENTS 

Requirements regarding areas of reinforcement in a beam are 
set out in CPI10 as follows: 

(1) Clause 3,11,4,1,. The area of tensile reinforcement should 
not be less than: 0,15% b^d for high yield reinforcement and 
0,25% b d fcr mild steel. Here b is the width of the 
section and d is the effective depth. For a T-section, b^ 
should be taken as the average breadth of the concrete below 
the upper flange, 

(2) Clause 3,11,5, The maximum area of reinforcement in either 
tension or compression should not exceed 4% of the gross 
cross-sectional area of the concrete. 

Requirements regarding areas of reinforcement in the columns 
are defined in clause 3,11,4,1; viz 

(a) The minimum area of main longitudinal bars is 1% of the gross 
area of the column 

(b) The minimum number of bars in a rectangular column is 4 and 
their diameter should not be less than 12 mm. 

The maximum area of longitudinal reinforcement is given in 
clause 3,11,5, CPllO, This should not exceed 6% of the gross cross- 
section in vertically cast columns. 




6o5 DESIGN LOADS 


158 

Unlike the Steel Code of Practice, factored loading is not 
allowed in CP110;1972„ Furthermore, the design loads for the 
serviceability limit states of deflection are not just the 
working loads. 

The characteristic loads are defined in clause 2,3,1 as: 

(a) Characteristic dead load which is the weight of the com¬ 
plete structure, 

(b) Characteristic imposed load which depends on the use of 

the building. These loads are given in CP3 Chapter V, Part 
1, for various buildings, 

(c) Characteristic wind load W^, This is defined and calculated 
in accordance with CP3(1970) Chapter V, Part 2, 

The design load is obtained by multiplying the characteristic 
load by the partial factor of safety which takes account of (1) 
possible overloads; (2) inaccurate assessment of the effects of 
loading and unforeseen stress redistribution within the structure; 
(3) variations in dimensional accuracy; and (4) the limit state 
being considered. The value of y^ varies for the different limit 
states as set out in Tables 6,1 and 6,2, 



TABLE 6,1: VALUES OF y^.SERVICEABILITY LIMIT STATE OF DEFLECTION 
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6o5o1 Design Loads for Serviceability Limit 

For the serviceability limit state of deflection there are 
three load combinations,, "fhe first load combination concerns 
vertical loads only and as sway is mostly due to wind load, the 
first load case is disregarded in this thesis. The second load 
combination gives 1,0 + 1,0 and there is no need to consider 

any change in the load arrangements. The third load combination 
gives 1,0 + 0,8 Qj^ + 0,8 W^, 

6,5,2 Design Load for Ultimate Limit State 

With regard to ultimate limit state, the arrangement of loads 

should cause the most severe forces. Under load combination (1) 

on the loaded spans there should be 1,4 Gj^ + 1,6 Q^, but only 

1,0 G, on the unloaded span. For example, for a three bay frame, 

K 

four arrangements should be considered. Figure 6,1 shows the 
different loading patterns for individual floors. The frame 
should be analysed four times, with the appropriate loading pat¬ 
tern applied on all floors. 


Load Combination 

Dead Load 
Max, Min, 

Imposed Load 
Max, Min, 

Wind Load 

1 - Dead and imposed 
load 

1,4 1,0 

1.6 0 

0 

2 - Dead and wind 
load 

1,4 0,9 

0 

1,4 

3 - Dead, imposed 
and wind load 

1,2 

1,2 

1,2 


TABLE 6,2: VALUES OF ^ULTIMATE LIMIT STATE 






1'60 


Under load combination ( 2 ), the most critical condition will 
have to be considered. This may arise when moments due to 1,4 
on some parts of the structure are added to the wind moments, and 
moment due to 0,9 on other parts of the structure from the 
restoring moment. For the three-bay frame, another four load cases 
should be considered. Figure 6,2 shows the different loading 
patterns for individual floors. It should be noted that the wind 
load for the four arrangements remains the same, i,e. 1.4 Wj^, 

Under load combination (3} a factor of 1.2 on all loads is 
used throughout the structure, with no variations for loaded and 
unloaded spans. This load case is similar to the one considered - 
in steel design and the two codes CPllO and the draft B/20 agree 
on the load factor of 1,2, 

It can be seen that for the three-bay frame, nine load cases 
should be considered, to obtain the maximum internal forces in the 
frame. The number of load cases obtained from the three load 
combinations depends on the number of bays in the frame, but at 
least four load cases are needed for a single bay frame, 

6,6 EQUATION FOR THE MINIMUM COLUMN DEPTH 

The iteration technique described in Chapter 2 is used to 
obtain the second moments of area for the columns so that the 
sway deflections are satisfied. These are then used to obtain a 
lower bound for the overall depth of each colxomn. 

For a rectangular section Figure 6,3, the uncracked second 
moment of area of the section is given by Reynolds and Steedman 
(1976) as: 

I = V b [x^ + (h-x)^] + (a - 1) 

XX J ^ 

[A (d-x)^ + A' (x-d')^] 
s ^ 


( 6 . 1 ) 



FIGURE 6„3; A TYPICAL COLUMN SECTION 


where 

which 

Here i 

and a 

I 

where 


X is the distance of the neutral axis below the top edge, 


is 


= J_ [i bh2 + (a - 1) (A d + A'd')] (6„: 

K 2 e ' ^ s s 

tr 


is the effective area given by: 

'tr 

A = bh + (a - 1) (A^ + A') (6,: 

tr e ^ ^ s s 


is the modular ratio, thus 


a = E /E 
e sc 


E is the modulus of elasticity of steel and 
s 

E is the modulus of elasticity of concrete. Values for 
c 

E are obtained from Clause 2.4.2.2 of CP110 for short 
c 

term loading. In these expressions A^ is the area of tensile 
reinforcements. A' is the area of compression reinforcements 
h is the overall depth of the column section, b is the width 
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of the section, d is the effective depth to tensile rein¬ 
forcement and d’ is the depth to compression reinforcement. 

When designing the column section, it is common practice to 
take the compression reinforcment equal to the tensile reinfor¬ 
cement (i.e, A^=Ap„ In fact, all the design charts for rectan¬ 
gular column cross-section in CP110;Part 2:1972 have been derived 
using A^ ~ ^s* 'this is the case, then equation (6.2) is 
simplified to x = h/2, and equation (6.1) becomes 

I = bhVl2 + 2A (a -1) (h/2 - d')2 (6,5) 

X A o 6 

Equation (6.5) could be written in terms of the percentage of the 


reinforcement p, thus 


'xx 

= bh3/12 +1^ (a^-1) ( ^ - hd* + d•^ 

(6.6) 

where 

p = (100 X 2A^)/(bh) 

(6.7) 

Let 

C = p (a -1)/100 

Jl 0 

(6.8) 


and with > E^ is always positive. Equation (6,6) can be 


written in terms of C^, thus: 

(1 + 3CJ - C.d'bh^ + C d'2 bh - I ^ = 0 (6.9) 

12 1 i 1 

Now v/ith = b (1 + 3C^)/12 (6,10) 

which is always positive, equation (6.9) becomes: 

C,d'b C^d' b 

h3 . -i— h2 * -ij— h - qyq = 0.0 (6.115 

2 2 

which can be written in the form 

^3 ~ C6ol2) 



163 


where 


Aj = -Cp'b/Cj 

(6.13) 

Aj = Cp'2b/C2 

(6.14) 

^ 

(6.15) 

An analytical solution for equation (6„12) is available 

(Spiegel, 

1968), and it depends on the value of the discriminant, 

defined as: 

This is 

D = 

(6.16) 

Q = ( 3 A 2 - a2)/9 

(6.17) 

and R = (9A^A2 - 21- 2a3)/54 

(6.18) 


It can be proved that D in this case is always a positive 
number. For a positive discriminant equation (6„12^ has only one 
real root, and two other complex conjugates„ This means that there 
is only one real solution for equation (6.12) which is 



h = S + T - 1/3 A^ 

(6.19) 

where 

s = 3/p + 

(6.20) 

and 

T = 3v/r - 

(6,21) 


Here R and D are as defined by equations (6.18) and (6,16) above. 

It should be noticed that to be able to calculate h the values of 
p, b and d’ should be known. These are given as data in the compu¬ 
ter program written for this purpose. 





6,7 DESIGN OF A BEAM SECTION 


The second moments of area for a beam section should be 
known before the iteration technique can be applied to obtain the 
column second moments of area. The width b of the beam section is 
first assumed, and then the effective depth d is taken as the maxi¬ 
mum obtained by the following two methods: 

1) The second moment of area I of the beam is obtained from 

Islam (1978) and Anderson and Islam (1979). The plain con¬ 
crete section is then used, only to obtain a preliminary beam 
depth from I = bh2/12 or h = ^/l2I/b. 

2) The effective depth is obtained from the limitation of the 
span/effective depth ratio. The allowable value for span/ 
effective depth ratio depends on: 

a) The span and the support conditions; 

b) The amount of tensile steel and its stress; 

c) The amount of compression steel; and 

d) The type of beam. 

The basic span/effective depth ratios for rectangular beams 
are given in table 8, CPI10, These are based on a beam with 1% 
tension reinforcement with a characteristic strength of 410 N/mm^, 
The limit on deflection is span/250, The values given depend on 

the support conditions and apply to beams up to 10 ra span. The 

basic span/ effective depth ratio from Table 8, CPI10 is 20 for a 

simply supported beam. To obtain a safe preliminary beam depth, 

this ratio is taken to be 20 for the beams of the frame. For a 

2 

steel with a characteristic strength of 250 N/mm , the span/ 
effective depth ratio is taken to be 29 according to Table 10, 


CPllO, 






FIGURE 6.4 A TYPICAL RECTANGULAR BEAM SECTION 


The larger value for the beam depth obtained from (I) or (2) 
above is used with a reinforcement of 1% to obtain the cracked 
second moment of area of the rectangular section. This is then 
used in the iteration technique to find the second moments of area 
of the columns so that the sway in each column is satisfied. 

6.7.1 Desi gn of the Beam Section for Ultimate Limit State 

When cliecking the beam section for strength requirements 
(ultimate limit state), the maximum sagging and hogging bending 
moments are needed. These are obtained from applying all the 
ultimate limit state load cases and then choosing the most severe 
bending moments. These are then used to obtain the reinforcements 
for the beam section. For the maximum hogging moment a rectangular 
cross-section (Figure 6.4) is used with width b, and effective 
dentil d obtained by (1) or (2) above in section 6.7. 


For the 




FIGURE 6.5 A TITICAL T-BEAM SEGTION AT MID-SPAN 

inaxiinuTn sagging moment a T section is used, where the beams are 
integrally cast and support a continuous floor slab. Part of the 
slab adjacent to the beam is thus regarded as the compression 
flange, see (Figure 6.5). In this figure is the effective 
width of the compression flange, b^^^ is the width of the web and 
h^ is the thickness of the floor slab. The effective width 
of- the comnression flange is specified in Clause 3.3.1.2 of CPllO 






For both rectangular and T-sections, the tensile reinforcements 
are considered firsto However, if the concrete in compression is 
overstressed, compression reinforcements is then added to assist 
the concrete to carry the extra loado 

Design charts for singly reinforced beams and doubly rein¬ 
forced beams are given in CP110:Part 2:1972„ All the charts have 
been derived using the assumptions given in Clause 3o3„5ol of 
CPI10:Part 1:1972 for the analysis of cross-sections, A computer 
subroutine was written by MacGinley (1978) to derive these charts. 
This however was limited to rectangular sections, with a constant 
beam depth. For this reason, MacGinley's flow-chart is altered 
in this thesis to design T-sections, and to increase the beam 
depth for cases when the maximum reinforcement in tension and 
compression is exceeded, 

6,8 ULTIMATE LIMIT STATE FOR THE COLUMNS 

Once the columns are designed for deflection limit state, 

the overall calculated depth h^ of a section and the area of 

reinforcement A , must be checked for the ultimate limit state, 
sd 

To do this, it is necessary to calculate the maximum bending 

moment M and the maximum axial load N for all the load 
max ro3.x 

cases. Checking is carried out by recalculating the area of 

reinforcement A for a column section of depth hj to provide the 

section with an ultimate resistance to M^^^ and The area 

A is then adopted if A > A This is provided that the area 
sc sc sa 

of the reinforcement is less than the maximum allowable 6%, If 
not, the computer programme increases the depth hj in increments 
of 50 mm and the area A is recalculated until it is below 6%, 

5 ^ 

The reinforcement for columns resisting moments and axial 
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forces can be either calculated using Clause 3,5,5,3 of CPI 10 or 
obtained using charts given in part 2 of the code. The charts 
are for rectangular parabolic stress distribution for concrete 
and a trilinear stress-strain curve for the steel reinforcements, 
as for beams (Clause 3,3.5,1, CPI10], A computer subroutine is 
written to derive these charts. This'makes use of certain signi¬ 
ficant points on the chart to accelerate the trial and error pro¬ 
cedure needed for the.design. To show how this subroutine works an 
explanation is given first of how the design charts are constructed 
and what properties they have. More information about the use 
of these charts is given in CPllOrPart 2:1972, (Allen 1977) and 
^MacGinley 1978), 

6,8,1 Construction of Design Chart in CPllO 

A symmetrically reinforced column section subjected to the 

ultimate axial load N and the ultimate moment M is con- 

max max 

sidered. The moment is equivalent to the axial load acting at an 
eccentricity e = M /N , Depending on the relative values of 

Tn3.X TTlS-X 

N and M , two cases occur for analysis, 
max max 

i) Compression over the whole section where the neutral 
axis is outside the section, 

ii) The section is divided into two parts by the neutral 
axis, in one part the steel is under tension while in 
the other part both steel and concrete are in compression. 
For a given location of the neutral axis, and for a given 
reinforcement, the strains and stresses in both the concrete and 
the steel can be determined, and from these the values for the 
internal forces can be found. The resultant internal axial force 
and resistance moment can then be evaluated. These should be 
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equal to N and M ,, Consider a section subject to an ultimate 
^ max max 

load N which is less than the capacity of the section under 
max 

pure axial loado This section is able to support ^ 

determinable maximum ultimate moment M ^ when the concrete is at 

max 

its maximum strain and design strength. The stresses in the 
steel depends on the location of the neutral axis. 

It can be noticed from the above that for a given percentage 
of reinforcement the value of the resisting moment and axial 
force for the section depends on the location of the neutral axis. 

This means that successive trials would be required to find the 
solution. To assist the engineer the different trials are readily 
calculated and plotted in the form of a design curve in part 2 of 

CPllOo 

For a selected grade of concrete and reinforcement, the 
design curve for a section with a given percentage of reinforcement 
100 A /bh and its location d/h symmetrically placed about the 

SC 

centre line of the section is formed, Tliis is done by plotting 
values of N/bh against M/bh^ for various positions of the neutral 
axis. It should be noticed that is the total area of rein¬ 
forcement and that this is divided equally between the faces 
parallel to the axis of bending, (e,g, in Figure 6,3 = A' = A^^/2). 

A typical design curve is shown in Figure 6,6 with the stress strain 
curve for the reinforcement shown in Figure 6,7,other curves can 
be constructed for percentages of steel ranging from zero to a 
maximum of 8%, The family of curves forms the design chart for a 
particular combination of materials and a typical chart is shown 

in Figure 6080 

In Figure 6,6, the specific points correspond to those in the 
for steel as shown in Figure 6,7, For example. 


stress-strain curve 





point A in the design curve corresponds to the final yield in the 
stress-strain curve. 

The strain in the steel is obtained from the stress strain 
curve (e.g, for mild steel the strain is 0,00309), While the 
strain in concrete is always assumed to be at its maximum value of 
0,0035 as specified in the code. The depth x of the neutral axis 
is calculated assuming a linear strain distribution as shown in 
Figure (6.9b) . From this: 


0,0035 _ 0,00309 

X 0,8h-x 


or 


X = 0,425h 


Similarly, for point B in the design curve which corresponds to a 
strain of 0,00087 at first yield, the strain diagram shown in 


Figure{6,9 c)gives the position of the neutral axis as: 


0.0035 _ 0,0 0087 

X 0.8h-x 

X = 0,646h 

Point c in the design curve corresponds to zero strain in tensile 
steel (Figure 6,9d), Thus the neutral axis is at the centre of 
the tensile steel with zero stress. For the part CDE of the design 
curve in Figure 6,6 the tensile steel is in compression. For the 
part FABC of the curve the steel is in tension. Point F of the 
curve is for X = h - d. This value of x is the minimum allowable 
in the code. No solutions are given for x less than h-d. The 
design curve is flattened off at part DE to ensure that a minimum 
allowable moment of 0,05 h is automatically included. 

From the above description of the design curve, it could be 

concluded that five significant points can be recognised. These 


are: 
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Cl) Point F where x = h d = d'„ This is the lowest allowable 
value of X 

(2) Point A where x is obtained for final yield of the tensile 
steel 

(3) Point B where x is obtained for first yield of the tension 
steel 

(4) Point C where x = d, and the stress in the tensile steel 
is zero. 

(5) Point D where x is obtained for the section resisting a 
minimum moment of 0.05 Nh. 

The design curve FABCDE is for a specific percentage of 
steel as 4%. Other design curves for the steel percentages of 
0, 1, 3, 5, 6, 7 and 8 may be plotted. These form the design 
chart shown in Figure 6.8, The various zones which depend on the 
location of the neutral axis are marked on the chart. The sig¬ 
nificant points on the design curves are marked, and equivalent 
points are connected by the dotted lines to form these zones. 

6.8.2 Computer Subroutine 

A computer subroutine is written in Fortran to determine the 

amount of reinforcement required in a symmetrically reinforced 

section to resist a given axial force and moment If 

M < 0.05 N h then M is taken as equal to 0.05 N. . A 
max max max max 

process or successive trials making use of the significant points 
on a design chart is used. The process is general, but it would be 
simpler to explain by the following example. 

6.8.3 Example 

Design the reinforcement of a column section subjected to an 

ultimate load of N /bh of 10.5 N/ram^ an ultimate moment of 

max 




2 2 

N/inin o The laaterials are grade 25 concrete and 
mild steel reinforcement. The location of the steel bars is given 
by d/h = 0„8o 

The design chart for the materials specified is shown in 
Figure 6,8, and the subroutine is searching for point R which is 
located as follows: 

(1) Set the percentage of reinforcement at the minimum value of 

1 %, 

(2) Set the depth of the neutral axis to correspond to point A in 
the design chart with x = 0,425 h, 

(3) Calculate the value of the resisting moment and axial force ■ 
Q and W, 

(4) If Q < M and W < N , then increase the percentage of 

II13.X TT13-X 

reinforcement by 0„25%„ 

(4a) Repeat step 4 until Q ^ M or W ^ N , or the percentage 

ms-X Tn3-X 

of steel is starting to become more than 6%, 

(5) If the percentage of reinforcement exceeds 6% and Q is less 

than M and W is less than N , then increase column depth 
max max 

by an increment of 50 mm. Steps 1 to 5 are repeated. 

(6) Depending on the values of Q and W one of these two cases 
occurs first, 

al 0 > M but W < N , This means that point R is above 

^ ^ max max 

the dotted line AA* on the design chart, 

bl 0 < M but W > N , This means that point R is below 

^ ^ max max 

the line AA', For the example under consideration, case (a) 

occurs. In Figure 6,8, as soon as the steel percentage is 

2 

increased from 3,50% to 3,75%, Q/bh will be equal to 3,6 and 
W/bh = 3,8, Thus 0 > M and W < N , This means that 

n/un 0,0, iiioo 

point R for this example is above the line AA' and the depth 
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of the neutral axis is more than 0 o 425 h, 

From the above steps it is only possible to know whether the 
neutral axis is below or above one of the dotted lines shovm in 
Figure 6,80 Repeating the steps with points B and C it is poss¬ 
ible to specify upper and lower boundaries for the depth of the 
neutral axiSo Once this is done, a new trial position of the 
neutral axis is obtained by: 

x^ = (x^ + x ^)/2 ( 6 . 22 ) 

where x^ and are the upper and lower boundaries respectively. 

For point R for instance, x^ = 0.8h, and X 2 = 0.646h„ Thus for 
the first trial x^ = (0.8h + 0.646h)/2 = 0,723h steps 1 to 6 are 
now repeated to determine whether x^ is a new lower or upper 
bound. Equation (6.22) is then applied with this value of x, 
replacing X 2 to calculate a new value for x^ and the process is 
repeated until two successive values of x^ are within a specified 
tolerance. 

The amovint of reinforcement is obtained by starting with 

A /bd as 1% and then increasing it in steps of 0.25% until 
sc 

0 = M and W = N .At this point the actual neutral axis is 
^ max^ max 

known as well as percentage of steel. A large number of columns 
designed by the above computer subroutine and it was found 
to be faster than that used by MacGinley (1978) . Usually five to 
seven trials are needed to determine the actual position of the 
neutral axis. 

5^9 COMPUTER PROGRAM FOR THE DESIGN OF REINFORCED CONCRETE 
FRAMES 

A program for the design of two dimensional, rectangular 





reinforced concrete, sway frames was written in Fortran and run on 
the ICL 1904S computer at the University of Aston, This program 
makes use of the design procedure described in this chapter. The 
program consists of a master segment calling a number of segments, 
the functions of which are similar to that of the steel program 
described in Chapter 4, However, in the program for the design 
of reinforced concrete frames, the beam sections are kept unaltered. 
The ratio between the second moment of area for internal and exter¬ 
nal columns is also kept constant. This is because there is more 
flexibility in the procedure for the design of reinforced concrete 
columns. For instance, the design for deflection limit state is - 
carried out with 1% reinforcement and rather than altering beam 
sections, the ultimate limit state design is obtained by increasing 
this percentage, 

6,9.1 The Program Segments for the Design of 
Reinforced Concrete Frames 

This consists of the master segment calling eight other seg¬ 
ments these are; 

1 - ITERATE 

2 - INTERNAL FORCES as in the steel program 

3 - STABILITY-F 

4 - The function CDEPTH: This calculates the depth of the column 

provided that the second moment of area of the column section 
and the reinforcement ratio are known. This segment makes use 
of equation (6,19), 

5 - The function SMAB: This calculates the cracked second moment 

of area for a beam section provided that the beam depth is known 
1 % reinforcement is assumed, Tne beam depth is calculated as 




described in section 6o7, 

6 - The subroutine U-BEAfl: This designs the reinforcement for a 

beam section for the ultimate limit state making use of 
section 6o7olo 

7 - The subroutine U-COLl)MN: This designs the reinforcement for 

a symmetrically reinforced column section for the ultimate 
limit state making use of section 6„8o2o This is done after 
the ml li m n depth is calculated from the limit state for defl¬ 
ection« 

8 - The subroutine SSSC: This is called by subroutine U-COLUMN 

to calculate the steel stresses from the stress-strain curves. 

The strain in the steel is calculated according to the position 
of the neutral axis. 

6a0 A 15 STOREY REINFORCED CONCRETE FRAME WITH 
THREE UNEQUAL BAYS 

The computer program described was used to design the 15 storey 
3 unequal bay frame shown in Figure 6.10, The frame is to be built 
of reinforced concrete. The characteristic strength of the concrete 
and the steel are 40 N/ram^ and 410 N/mm^ respectively. The applied 
loads on the frame are shown in Figure 6.10. The design loads are 
obtained by making use of sections (6.5, 6.5.1 and 6.5,2). Hence 
two design load cases for deflection limit state, and nine others 
for the ultimate limit state are obtained. 

The limit on the sway is taken to be equal to the height of 
the column/500. The deflection profile used in the design is shown 
in Fio-ure 6.11. The width of the beam and column sections is 
taken equal to 200 mm for all storeys. The flange depth h^ for 
T-beams is taken as 120 mm, while the flange width B^ equals 830 mm. 
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The covers of tensile steel (d’) and compression steel (G) are 
33 mm, and 28 mm respectively. The percentage of the reinforcement 
assumed to resist the sway deflection is 2% for the columns at 
levels 1 to 5„ 1„5% for the levels 6-10, and 1,0% for the columns 
of the top five storeys. The cover for this reinforcement d’ is 
36 mm for all columns, 

A typical storey is represented in Figure 6,12 in which the 
sections at which the reinforcments are calculated are shown. In 
each beam three sections are considered. These are A, B and D, 

For each column one section is considered this is the section with 
the maximum bending moment, which is at one of the column ends. 

The design of the beam sections for the limit state for 
deflection is shown in Table 6,3, All beams have the same section 
and reinforcement. This is because the beam depth obtained from 
the limitation of the span/effective depth ratio dominates the 


Span under 
consideration 

H 

mm 

B 

ram 

Reinforcement 
at A, B and D 

, 2 

A mm 
s 

2 

A' mm 
s 

Left-span beam 
sections 

500 

200 

934 

233 

Mid-span beam 
sections 

300 

200 

534 

133 

Right-span beam 
sections 

400 

200 

734 

183 


TABLE 6,3: 15 STOREY FRAjME - THE DESIGN OF 

THE BEAM SECTIONS FOR THE LIMIT 
STATE OF DEFLECTION 


design. For all the beams, the limitation gives larger beam depths 





than those obtained using Anderson and Islam's equation. The 
percentage of the reinforcement (A^/bd] is assimed to be 1% for all 
the beams. This is the basic percentage referred to in CPllO 
(1972), at which the span/effective depth equals 20 for a simply 
supported beam. The beam sections for the final design which 
satisfy the ultimate state as well as the limit state for deflection 
and cracking remain unchanged. However the steel area is increased 
for most of the beam sections. The steel area for beam sections 
A, B and D of the final design are shown in Table 6.4, It can be 
seen that some reinforcement area are determined by the limit 
state of deflection. These are marked by an asterisk in Table 6,4, 
The design of the column sections is shown in Table 6,5 for 
deflection limit state and the ultimate limit state. It can be 
seen from the table that the ultimate limit state is governing the 
design and most of the column depths obtained from this state are 
more than those obtained from the limit state for deflection. 
Furthermore the reinforcement is also increased considerably. 

It should be noticed, that the column depth obtained for 
deflection limit state is considerably greater than the column 
depth specified in CPllO:1972, The code considers a value of 
effective height/30 to be acceptable. This underestimates the 
column depth needed to satisfy the limit state for deflection. 

For example in this frame, if the effective height is taken equal 
to the actual height (i.e, 3500), then the minimum column depth 
specified in the code is 3500/30 = 150 mm (rounded to the nearest 

50 mm) and except for the top three storeys this is much less than 
the column depth proposed by this new design method. 



6oll A 24 STOREY REINFORCED CONCRETE FRAME 
WITH TWO UNEQUAL BAYS 

The 24 storey - two unequal bay frames shown in Figure 6„13 

was designed using the computer program described in section 6„9„ 

The frame is to be built of reinforced concrete. The characteri- 

2 

Stic strength of the concrete and the steel are 40 N/inra and 410 
2 

N/nrm respectively. The applied load on the frame are shown in 
Figure 6.13, The design loads are obtained as described in 
sections (6.5, 6.5.1 and 6.5.2). This gives two design load cases 
for the limit state for deflection, and five other cases for the 
ultimate limit state. 

The non-linear deflection profile (used in the design) and the 
linear profile are shown in Figure 6.14, in which the sway is 
restricted to h/500. However, for selecting the beam sections, 
Anderson and Islam’s equations are used with sway restricted to 
h/1000. It should be stressed that this latter restriction is 
only applied for the purpose of selecting the beam sections since 
Anderson and Islam's equations underestimate these. 

To calculate the column depths, and the reinforcement of the 
beams and the columns certain sectional properties need to be 
supplied. These are: 

1) The covers to the reinforcement, which are given in Table 6.6 
making use of the recommendation given in CPI10:1972, clause 
3.11.2. In this table d' and G are those shown in Figures 6.3 

and 6o4o 

2) The width of the beam and column sections, which is taken equal 
to 250 mm for storey levels 1-18 inclusive. This is reduced to 


200 mm for the rest of the frame. 





3) 


For T-beam sections the flange depth h^ Csee Figure 6,,5) is 
taken equal to 120 mm for all the storeys. The flange width 
is 880 mm for storey levels 1-18, but is reduced to 830 ram 
for the remaining storeys, 

4) The percentage of reinforcment (lOOA^^/bh) assumed to resist 
the sway of the column. These percentage are given in Table 
6,7, 

A typical storey is represented in Figure 6,15 in which the 
sections at which the reinforcements are calculated are indicated. 
There are three sections in each beam A, B and D^and one section 
in each column. 

The design of the beam sections for the limit state for defl¬ 
ection is shown in Table 6,8, It is noticed that (except for the 
top three storeys) the beam depth designed by Anderson and Islam's 
equations is larger than that obtained from the limitation of the 
span/effective depth ratio. The percentage of reinforcement (A^/bd) 
assumed to be 1% for all the beams. The steel areas obtained 
using this percentage are given in Table 6,8, For the compression 
steel the minimum percentage of 0.25% is used. 

The beam sections for the final design which satisfy the 
ultimate limit state as well as the limit state for deflection and 
cracking remain unchanged. The area of the tensile steel for this 
design is shown in Table 6,9, Sections A, B and D in the table 
refer to those of Figure 6,14, The area of the compression steel 
remains nominal at the minimum of 0,25%, It can be seen that many 
tensile steel areas are determined using the deflection limit state. 

These are marked by asterisks in Table 6,9, 

The design for the column sections for the deflection limit 




state is shown in Table 6,10, This also shows the final design, 

in which the ultimate limit state as well as the limit states for 

deflection and cracking are satisfied. Unlike the 15 storey 

frame, most of the column depths obtained by the limit state for 

deflection are adequate for the ultimate limit state, and no 

increase in the column depth is required. The unchanged column 

depths in the final design are enclosed in brackets in (Table 

6ol0). For some sections the deflection limit state governs the 

design, where the column depth and the steel area are determined by 

this state. For example, column at level 14 (see Table 6.10) 

has a column depth of 450 mm. The area of reinforcment, resisting 

2 

the sway is 2250 mm , which is more than that required for the 

ultimate limit state. This means that the final design should have 

2 

an area of steel equal to 2250 mm to satisfy deflection as well 
as strength requirements. The steel areas determined by the limit 
state for deflection and not by the ultimate limit state are 
marked by asterisks in Table 6.10, However, most of the sections 
need a larger steel area for the ultimate limit state. 

It is noticed, that the column depth obtained for the limit 
state for deflection is considerably greater than the column depth 
specified in CP110:1972, clause 3.5,8, which gives a column depth 
of 150 mm which is adequate only at the top storey. 




















































101.61 mm 



FIGURE 6.11 SWAY DEFLECTION CURVE FOR THE 15 STOREY FRAf4E 
SHOWN IN FIGURE 6.10 





















Steel area 


Storey 

level 


Span 





Left 1214 233 

Mid 854 133 

Right 1159 183 


Left 1282 233 

Mid 847 133 

Right 1115 183 


Left 1329 233 

3 Mid 867 133 

Right 1107 183 


Left 1385 233 

4 Mid 895 133 

Right 1112 183 


Left 1417 233 

Mid 927 133 

Right 1126 183 


Left 

Mid 

Right 


Left 

Mid 

Right 


Left 

Mid 

Right 

Left 

Mid 

Right 

Left 

Mid 

Right 

Left 

Mid 

Right 

Left 

Mid 

Right 


1413 233 

952 133 

1147 183 


1358 233 

981 133 

1151 183 


1270 233 

996 133 

1157 183 

1163 233 

1008 133 

1158 183 

1034 233 

1016 133 

1155 183 

*934 233 

1078 133 

1141 183 

*934 233 

1227 133 

1122 183 


Steel 

Sivea,^ 

Steel 

area2 

Section 

B mm 

Section 

D nun 

A 

A' 

A 

A' 

s 

s 

s 

s 

1459 

233 

2891 

883 

763 

133 

1696 

555 

1117 

183 

2258 

678 

1426 

233 

3089 

1110 

764 

133 

1824 

701 

1096 

183 

2465 

916 

1428 

233 

3093 

1115 

764 

133 

1842 

722 

1097 

183 

2484 

937 

1435 

233 

3037 

1051 

764 

133 

1820 

697 

1102 

183 

2445 

893 

1462 

233 

2925 

921 

765 

133 

1766 

635 

1118 

183 

2346 

799 

1503 

233 

2772 

746 

765 

133 

1696 

554 

1142 

183 

2201 

612 

1557 

233 

2609 

558 

765 

133 

1613 

459 

1174 

183 

2041 

428 

1619 

233 

2444 

369 

765 

133 

1538 

373 

1211 

183 

1870 

231 

1690 

233 

2280 

233 

764 

133 

1467 

291 

1254 

183 

1695 

183 

1773 

233 

2100 

233 

763 

133 

1396 

210 

1306 

183 

1389 

183 

1873 

233 

1783 

233 

760 

133 

1327 

133 

1371 

183 

1120 

183 

1995 

233 

1519 

233 

755 

133 

1263 

133 

1452 

183 

857 

183 


. marks steel area determined by the limit state of deflection 

lABLE 6.4 15 STOREY FR.4HE - THE STEEL AREA FOR BEAM SECTIONS 

OF THE FINAL DESIGN (font) 





































■ 

Steel 

area ^ 

Steel 

area ^ 

Steel 

area ^ 

Storey 

Spa 


Section 

A mm 

Section 

B mm 

Section 

D mm 

level 


A 

A’ 

A 

A' 

A 

A' 




■ 

s 

s 

s 

s 

s 

s 


Left 

*934 

233 

2148 

233 

1248 

233 

13 

Mid 

1 

1330 

133 

747 

133 

1192 

133 

Right 

1084 

183 

1559 

183 

*734 

183 


Left 

*934 

233 

2355 

233 

997 

233 

14 

Mid 


1472 

133 

732 

133 

1260 

133 

Right 

1020 

183 

1713 

183 

*734 

183 


Left 

*934 

233 

2558 

233 

*934 

233 

15 

Mid 


1610 

455 

715 

198 

1373 

133. 

Rig 

ht 

946 

183 

1874 

183 

*734 

183 


TABLE 6.4: 15 STOREY FRAME - THE STEEL AREA FOR BEAM SECTIONS 


OF THE FINAL DESIGN 





















Storey 

Column 

Limit state of 
deflection 

Ultimate limit 
state 

level 

section 

A 2 

h A mm 

mm sc 

2 

h A mm 

mm sc 


Cl 

550 

2200 

600 

3575 


c' 

600 

2400 

900 

6400 

1 

7 

c: 

600 

2400 

750 

5525 


c" 

^4 

500 

2000 

500 

3250 


c, 

550 

2200 

550 

4675 1 

2 

Co 

600 

2400 

850 

7225 

C^ 

600 

2400 

700 

6300 


L 

^4 

500 

2000 

500 

4500 


Cl 

550 

2200 

550 

5700 1 


cj 

600 

2400 

800 

7650 

3 

C^ 

600 

2400 

650 

7500 


"4 

500 

2000 

500 

4750 


c. 

500 

2000 

500 

3825 j 



600 

2400 

700 

6500 

4 

c5 

550 

2200 

600 

4950 


^4 

450 

1800 

■ 450 

3800 


Ci 

400 

1600 

450 

4000 1 


ci 

500 

2000 

650 

6000 

5 

C^ 

450 

1800 

550 

4750 


C^ 

^4 

400 

1600 

400 

4200 


c. 

350 

1050 

400 

3400 1 


ci 

400 

1200 

600 

5500 

6 

c^ 

400 

1200 

500 

4500 


C^ 

^4 

300 

900 

350 

3325 


Cl 

300 

900 

400 

3400 1 


ci 

350 

1050 

550 

5500 

7 


300 

900 

450 

4500 


C^ 

^4 

300 

900 

350 

3325 


c. 

250 

750 

350 

3500 1 


Co 

300 

900 

500 

5000 

8 

ct 

300 

900 

400 

4200 


3 

4 

250 

750 

350 

2625 


c. 

250 

750 

400 

1800 


ci 

250 

750 

450 

4725 

9 

c^ 

250 

750 

350 

4025 


3 

C 

^4 

200 

600 

300 

2850 


TABLE 6.5: 15 STOREY FRAME - DESIGN OF THE COLUMN 


SECTIONS FOR DEFLECTION LIMIT STATE AND 


ULTIMATE LIMIT STATE 


























Limit state of Ultimate limit 

Storey Column deflection state 

level section -^-—-j- 

h A mm h As mm 

mm sc mm c 


300 

3000 

400 

4200 

350 

2800 

300 

2100 

300 

2100 

400 

2800 

300 

2550 

250 

2375 

250 

2125 

300 

3450 

250 

2125 

250 

1500 

200 

2200 

250 

3000 

200 

2100 

200 

1600 

200 

800 

200 

2200 

200 

800 

150 

1800 

150 

300 

150 

825 

150 

300 

150 

300 


TABLE 6.5: 


15 storey ER,/V^'E - DESIGN OF THE COLUMN 
SECTIONS FOR DEFLECTION LIMIT STATE AND 
ULTIMATE LIMIT STATE (Cont) 


































FIGURE 6.14 24 STOREY CONCRETE FRAME - SWAY DEFLECTION PROFILES 




Storey 

level 

Covers for beam 
reinforcement 

Covers for 

COlumn 

reinforcement 

d’ 

G 

d' 

1 to 12 

36 mm 

30 mm 

40 mm 

13 to 24 

33 mm 

28 mm 

36 mm 


TABLE 6.6: 24 STOREY FRAME - REINFORCEMENT 

COVERS 


Storey 

level 

Reinforcement ratio of the 
column sections (100 A^^^/bh) 

1 to 18 

2.0 

19 - 21 

1,5 

22 - 24 

1,0 


TABLE 6,7: 


24 STOREY FRAME-REINFORCEMENT 
RATIOS ASSUMED FOR THE LIMIT 
STATE OF DEFLECTION 
























Reinforcement at 

Storey 

level 

Span 

under 

Beam 

Depth 

A, 

B and D 



Consideration 

h mm 

o 

L 

, 2 



A mm 
s 

mm 

1 - 4 

Left 

750 

1785 

446 


Right 

500 

1160 

290 

5-9 

Left 

700 

1660 

415 


Right 

500 

1160 

290 

10 - 12 

Left 

Right 

650 

450 

1535 

1035 

383 

258 

13 

Left 

650 

1542 

385 

Right 

450 

1042 

260 

14 

Left 

600 

1417 

354 

Right 

450 

1042 

260 

15 

Left 

600 

1417 

354 

Right 

400 

917 

229 

16 - 18 

Left 

i 

550 

1292 

323 

Right 

400 

917 

229 

19 

Left 

500 

1034 

258 

Right 

400 

734 

183 

20 - 21 

Left 

500 

934 

233 

Right 

350 

634 

158 

22 - 24 

Left 

500 

934 

233 

Right 

300 

534 

133 


TABLE 6„8: 24 STOREY FRAiME - THE DESIGN OF THE BEAM 

SECTIONS FOR THE LIMIT STATE FOR 
DEFLECTION 























































Storey 


A A 

^ s^. s . 

at section at section 

A mm 13 2 


at section 
D 



* Marks steel 
deflection. 

TABLE 6.9: 24 ST 

BE.AM 


area determined by the limit state for 

OREY FRAME - THE STEEL AREA FOR 
SECTIONS OF THE FINAL DESIGN 















































































































A 

A 

A 

Storey 

level 

Span 

s . 

at section 

A mm^ 

s . 

at section 

B 

s . 

at section 

D mm^ 

18 

Left 

*1292 

1400 

2272 

Right 

*917 

*917 

1864 

19 

Left 

1090 

1390 

2295 

Right 

*734 

*734 

1750 

20 

Left 

1405 

1534 

2361 

Right 

*634 

661 

1652 

21 

Left 

1359 

1608 

2237 

Right 

*634 

681 

1443 

1 

22 

Left ■ 

1104 

1772 

1925 

Right 

841 

831 

1065 

23 

Left 

*934 

2073 

1398 

Right 

1091 

887 

*534 

24 

Left 

*934 

2446 

957 

Right 

1455 

943 

*534 


* Marks steel area determined by the limit state for 
deflection 


TABLE 6 <,9: 24 STOREY FRAME - THE STEEL AREA FOR BEAM 

SECTIONS OF THE FINAL DESIGN 

































The design for 

the limit state Final design 
Storey Column. . of deflection 
level section———---- 


mm 

, 2 
A mm 
sc 

h 

mm 

, 2 
A mm 
sc 

600 

3000 

750 

10781 

650 

3250 

950 

13656 

500 

2500 

(500) 

6875 

550 

2750 

700 

8750 

600 

3000 

900 

11813 

450 

2250 

(450) 

5906 

500 

2500 

700 

7438 

600 

3000 

850 

11156 

400 

2000 

450 

5063 

500 

2500 

650 

7719 

600 

3000 

800 

11000 

400 

2000 

(400) 

5750 

500 

2500 

650 

6906 

600 

3000 

800 

10500 

400 

2000 

450 

4219 

500 

2500 

600 

6750 

600 

3000 

750 

10313 

400 

2000 

(400) 

5250 

500 

2500 

550 

6875 

600 

3000 

700 

10500 

400 

2000 

(400) 

5000 


500 

2500 

550 

600 

3000 

700 

400 

2000 

(400) 

500 

2500 

(500) 

600 

3000 

650 

400 

2000 

(400) 

500 

2500 

(500) 

600 

3000 

650 

400 

2000 

(400) 


Unchanged coluran depths are enclosed in brackets 

TARTF 6 10- 24 STOREY FRAIC - DESIGN OF THE C0LU1-1N 

TABLE b.iu. limit STATE OF DEFLECTION 

_ ,,, T rUTT CTATP rrnnt . "1 













Storey 

level 

Column 

The design for 
the limit state 
of deflection . 

Final design 

section 

h A mm^ 

mm sc 

h A mm^ 

, mm sc 

11 

cl 

A 

^3 

500 2500 
600 3000 
400 2000 

C500:) 4063 

C 6 OO:) 8525 

C400) 2750 

12 

"cl 

2 

450 2250 
550 2750 
400 2000 

(450) 4219 

600 7500 

(400) 2250 

13 

"cl 

? 

^3 

450 2250 
500 2500 
350 1750 

(450) 3375 

600 6000 
(350) 2844 

14 

Cl 

C 2 

s 

450 2250 
500 2500 
350 1750 

(450) *2250 

(500) 6563 

(350) 2406 

15 

"cl 

^3 

450 2250 
500 2500 
350 1750 

(450) *2250 
(500) 5625 
(350) 1969 

16 

C^ 

2 

"3 

450 2250 
500 2500 
350 1750 

(450) *2250 

(500) 4688 
(350) *1750 

17 

Cl 

C 2 

^3 

400 2000 

450 2250 

350 1750 

(400) 2500 

(450) 4781 

(350) *1750 

18 

Cl 

r 

^^3 

350 1750 
400 2000 
300 1500 

(350) *1750 
(400) 5000 
(300) 1875 

19 

C, 

0 

350 1050 
400 1200 
300 900 

(350) 2100 
(400) 4800 
(300) 1800 

20 

cl 

C^3 

350 1050 

400 1200 

300 900 

— - — ------ 

(350) 2100 
(400) 4000 

(300) 1500 


marks sreei cti ca ^ i 

deflection and not by ultimate limit state 

24 STOREY FRAME - DESIGN OF THE COLUMN 
SECTIONS FOR LIMIT STATE OF DEFLECTION 

AND ULTIMATE LIMIT STATE (Cont.) 


TABLE 6.10: 

















storey Column 
level section 


The design for 
the limit state Final design 
of deflection 


2 2 
A mm h A mm 
sc mm sc 


900 

(300) 

2700 

1050 

(350) 

4025 

750 

(250) 

1875 

500 

300 

1800 

500 

300 

3300 

400 

200 

1000 

400 i 

(200) 

500 

400 

250 

2750 

300 

200 

400 

300 

(150) 

1725 

300 

(150) 

900 

300 

(150) 

*300 


TABLE 6a0: 24 STOREY FRAME - DESIGN OF THE COLUMN 

SECTIONS FOR LIMIT STATE OF DEFLECTION 
AND ULTIMATE LIMIT STATE 




CHAPTER SEVEN 


DESIGN OF SPACE FRAMES 


7,1 INTRODUCTION 


The iteration technique given in Chapter 2 can in fact be 
generalised for the design of space frames. In this chapter, it 
is used to design pin-jointed space frames in which the deflections 
are specified at certain joints in these frames. The aim here is 
to select member areas so that the stress requirements are satis¬ 
fied while the deflection at the specified joints are allowed to 
take place. 

In general the deflection, usually at the middle of certain 
spans, are specified by codes of practice. For instance, the 
British Standard draft B/20 restricts the deflection in the middle 
span of a roof to the span/360, For such a frame, the procedure is 
to calculate the member areas so that the deflection in the middle 
span will be equal to the specified limit. Naturally, the strength 
requirements in some of the members will not be satisfied. The 
cross sectional areas of these are therefore increased and the 
interaction is repeated so that the areas of other members are 
reduced. The process is repeated until the deflection as well as 
the strength requirements are satisfied throughout. 


,2 The Contributions of a Member to the Overall 
Stiffness Matrix 

In a given pin-jointed space frame, let a general member i 
e connected to joint j at its first end and b at its second 
nd. This member is shown in Figure 7.1, in which the arrow on 
he member specifies its positive longitudinal P axis. The system 








coordinates XYZ are also shown in the figure. A right handed 
system is used for the X, Y and Z axes. 


205 


The contributions of member i to the overall stiffness matrix 
[K] of the space frame is given below by equation (7.1) on the 
next page (Majid, 1978). In equation (7.1) {w . w . w ...o.. 

^1 yi ^3 


w , w . -1 J J-J UHV./ VV./V^L.WX 1 . unv^ .A. U X 1 i Cl. X XWCl-VX-J c*. p .L-I-V-'CJ. ^ 

xk 3 k zk 

joints. Each load w has two subscripts, the first denoting 
the direction of the force, and the second, either j or k, denoting 
the joint at the end of the member being considered. The direction 
cosines of the member are 2.^, m^ and n^. It is noticed that only 
the direction cosines of the P axis relative to the X axis are 
involved. In equation(1.7 )a=EA/L is the axial stiffness of the 
member. Where A is the area of the member, L is the length, and 
E is the elastic modules. The column vector {x^ y^ ..... /-^ 

z } is the joint displacements vector. 

It is useful to note that equation (7.1) does not alter if 
the arrow on the member is reversed in Figure 7.1 so that j becomes 
the second end of the member. The direction cosines change their 
sign but the signs of m^, ^^m^ etc do not change. This obser¬ 
vation is useful when deriving the modified stiffness equations 
for the structure, 

7„3 modified stiffness EQUATIONS F OR A GENERAL JOINT IN THE 
PIN JOINTED SPACE FRAME 

The relationship between the external loads W and the vectori- 
ally equivalent joint displacements A, in terms of the system co¬ 
ordinates XYZ, is expressed by the equation 

. (7.; 

w = K A 

and z as three displacements in X, 1 
The vector A has Xj, Xj and z. 


w } is the vector of the external loads applied to the 
ZK 


(7.2) 


and Z directions o 


f a general joint j. Thus, three equations are 





THE CONTRIBUTIONS OF A SPACE FRAME MEMBER TO THE OVERALL STIFFNESS MATRIX 
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needed to define the displacements of each joint. 

In the design problem, the deflection at a joint in a given 

direction may be specified. The corresponding stiffness equation 

is then modified so that it is used to calculate the areas of the 

members connected to that joint. The modified stiffness equation 

for X. is obtained as follows: 

J 

The contributions of a member i to the deflection are cal¬ 
culated from the first row of equations (7,1), thus 


i.ji^.x. + a-.i .m .y- + •'ri -z. - a. £^.x, 

1 pi j 1 pi pi-^l 1 pi pi 3 1 pi k 


a.ji.m.y, -a.£.n.z = w. 
1 pi pi k 1 pi pi k xj 


(7,3) 


Rearranging equation (7,3) for x- we obtain 


X. = r-w . + a.£ .m .y. + a.il .n z - ax 
j xJ 1 pi pi 3 1 pi pi 3 1 pi K 

- a.£ .m .y, - a.£ .n . z )/C-a. £2 ) 

1 pi pi k 1 pi pi k 1 pi 


(7»4) 


Equation (7.4) gives the deflection x^ when one member only is 
connected to the joint. If R members are connected to j then each 
will have the same pattern of contributions, and the modified 
stiffness equations for becomes. 


= l-"xj " iij 

a £ m y - a.£ .n .z )] / ^^i^ni^^ 

i pi pi^k 1 pi pi k ^ 1-1 1 pi 

Similarly, the modified stiffness equation^ for y^ and z^ are. 


(7.5) 


j = [-'V3 " 


(7,6) 
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z. = [-W . + .1, (a.Jl .n .X. + a.m .n .y. - a.H .n . x, 

J ZJ 1=1 1 pi pi J 1 pi pi-';] 1 pi pi k 

- a.m .n .y - a.n^.z )] / [- (a.n^.) 

1 pi pi k 1 pi k^^ *■ 1=1 ^ 1 pi^ 


C7.70 


7„4 DESIGN EQUATION FOR MEMBER AREAS 

If one of the deflections Xj, y^ or is specified, then 
equations (7.5), (7,6) or (7,7) can be used to find the member 
areas. To do this, the axial stiffness of member i is written 
as a factor of the area A^, thus 


a. = A. b. 
1 11 


(7o8) 


where b^ is equal to E/L^. 

In practice, for construction purposes the members are grouped 
together so that the members of a group have the same area A^. 

For a total of n groups with areas .A^^, ,,,», A^^, ,,o,, A^, 




a = .A / A , 
n gn gl 


, a - A / A , ^ ,, o o, 
’ r gt gl 


(7o9) 


Hence, 
be wri 


the axial stiffness of member i given in equation (7,8) can 
tten in terms of the area of the first group of members as: 


^i “ir '^gl ^i 

where the area A- of member i, which is in group r, is equal to 
a. A Each ratio u has two subscripts, the first i denoting 

thl number, and the second x, denoting the group to which 

member i belongs. 

For a specified deflection in a given direction the design 


(7,10) 
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equations are written to calculate the area o£ the meinber in the 
first group. The areas of the other groups are then obtained from 
equation (7,9), 


7,4,1 Calculation of A 


If the deflection Xj in the direction of X axis is specified, 

then equations (7,4), and (7.5) can be used to find the area of 

the members in the first group. Equation (7,4) is rewritten in 

terms of A , and a. as: 
ai ir 


A , = w ./(a. b. £^.x.+a. b. £.m.y. 
gl xy ir 1 pi j ir 1 pi pi j 

+ a. b.il .n .z. -a. b-Jl^.x, - a. b.Jl .m.y, - a. b. .n . z, ) 
ir 1 pi pi j ir 1 pi k ir 1 pi pi k ir i pi pi k 

When R members are connected to joint j then equation (7,11) 
becomes: 

■'gl “ “x/[jil 

- - “lr*’l'pi'"pi>'k - “ir'>l*pi'’pi"k’i 


(7,11) 


(7,12) 


Similarly equations (7.6) and (7.7) become: 


'gl = (“lr'’i*pi'"pl'^i * “lr'’l”pi>'i * “ir'>l"pi'”pi"i 


K n - tn Y - n b V - ct. b.n .m . z, ) 
- a. _b.5, ^,m ,x^ ir i pr k ir i pi pi irJ 


ir 1 pi pi 


(7,13) 


A w fr, b J, n X. + a. b.m .n .y. 

A A = '^zj/i=l pi pi 1 ir 1 pi pi 3 


(7,14) 


bxm_,.n ,y,, - b,n2 z ) 


■ir^'pib - “ir^^pi'pi^k - “ir°i"'pi'Vi^k ' UrU'pi^k 
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It should be noticed, that in equations (7,12), (7,13) and 
(7,14) the area of the members of the first group A , cannot be 

o 

calculated if joint j is unloaded in the direction of the specified 
deflection. However, it is unlikely for a joint to have a critical 
deflection in the unloaded direction. 

The ratios , a„, „,,,, a , ,,,„, a are not knom in the 
first design cycle. For this reason, in this cycle, it is assumed 
that = 02 = ..o = = 1, Obviously, this means 

that for the first design cycle all the member areas are equal and 
the frame is constructed from the same group. The ratios a are 
changed in the subsequent design cycles. 


7,5 DESIGN PROCEDURE 

The design procedure for a pin-jointed space frame is to some 
extent similar to that of a steel sway frame, described in section 
2,8, This consists of the following steps: 

STEP 1 - Specify the joint at which the deflection is expected to 
be critical, 

STEP 2 - Specify the direction (X, Y or Z) and the value of this 
deflection, 

STEP 3 - Set each ct to unity to initiate the first design cycle, 

STEP 4 - Assume the area of the members in the first group, and 

the deflection of all the joints. Infeasible values for 
the deflections such as x^ = " ^j " ^ accep¬ 

table for the purpose of starting the iteration process. 

The area A could have an infeasible value too, but it 

was found that convergence is achieved more rapidly with 

2 

a large A^^^, e,g« 10 mm , 

STEP 5 - Using the modified stiffness equations, calculate the joint 
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STEP 6 


STEP 7 


STEP 8 


STEP 9 


displacements for all the joint up;to that in which the 
deflection is specified. 

Using the appropriate design equation, now calculate the 
first new value of A , and hence calculate the new member 
areas in all the other groups using equation (7,9), 

Continue with the iteration and calculate a new value 
for each unknown joint displacement. At this stage all 
the modified stiffness equations for all the joints have 
been utilized. 

Repeat steps 5, 6 and 7 until convergence is achieved. 

This happens when (the new value of A^^ - the previous . 
value of A^^)/(the New value of A^^) is less than a 
specified tolerance. Notice that the convergence test 
could also be applied to the joint displacements. At 
this stage of the design procedure the first design 
cycle is complete and a new area is provided for each 
member. All the joint deflections are also found. Notice 
that for this first design cycle all the calculated 
member areas have the same value. The member areas and 
the joint deflections calculated fully satisfy the 
stiffness equations and the design obtained in this 
manner is feasible. However, the design process continues 
for the following reason: 

a) To check that while deflection requirements are 
satisfied, the stresses are satisfactory in all the members 

b) To achieve economy by altering the ratios between 

the areas of the groups. 

Check that the deflection of each joint in the frame is 
less than that specified in Steps (Ij and (2) for a 
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given jointo If the maximum deflection 


appears^ [in ja 


joint other than the one specified in step 1, then 
repeat steps 1 to 8 with the deflection limitation 
applied to this new jointo 

STEP 10 “* Using the joint displacements obtained in steps 8 and 9 
calculate the internal forces in the members of the pin 
jointed frame« The force p^ in a member i is expressed 
in terms of the displacements of the joints to which 
it is connected as 

p. = (EA./Lp (-x.lpj -y.mp. -z.Hpj ♦ ♦ y^nip. * 


( 7 . 15 ) 


STEP 11 - Design the member areas to meet the stress requirements» 
The design stresses are those given in British Standard 
* Draft B/20. Notice that two design stresses are used. 

The first is for a member under tension, and the second 
is for compressiono As a result new areas for the mem¬ 
bers in each group are obtained, and the new ratios for 
each group are calculatedo 

STEP 12 - Repeat the iteration process, ioe, steps 5, 6 and 7 until 
convergence is obtained once again. Notice that when 
entering the iteration this time, the current values of 
the joint displacements and member areas are employed. 
These values are more realistic than the initial infea¬ 
sible ones. 

STEP 13 - Repeat steps 10 to 12 until the values of a^, ... , 




become stable. This happens when: 


new value of ct - old value of a | ^ tolerance 


New value of a 
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At this stage the design is complete. The final design obtained 
in this manner satisfies: 

a - The deflection requirement at the critical joint; 
b - The strength requirement, as the stresses in all the members 
are less than or equal to the design stresses; 
c - Design economy is achieved, as this design satisfies the 
deflection requirements, and at the same time the cross 
sectional areas of the members are distributed according to 
their stress requirements. In this manner the member areas 
are distributed more efficiently, 

7,6 COMPUTER PROGRAM 

A computer program based on the design procedure just described 
was written in FORTRAN and run on the ICL 1904S computer. The 
program consists of three segments. These are: 

1 - The MASTER segment: This follows the steps of the design pro¬ 

cedure of Section 7,5, 

2 - The subroutine ITERATE, This calculates the area of the 

members in the first group and the joint displacements, A 
flow-chart of this subroutine is given in Figure 7.2 and it 
will be explained in more detail later, 

3 - The subroutine INTERNAL-FORCES: Tliis calculates the force in 

each member using equation (7,15), A positive force means 
that the member is in tension. 

The data preparation for this program is described in Appendix 
C of this thesis. The data given is similar to that used by the 
PAVO ANALYSIS SYSTEM, Space Structures Research Centre, University 
of Surrey, Obviously, extra data should be supplied to specify 
the joint at which the critical deflection is expected to occur and 
the direction of this deflectioHo 
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7o6,l Subroutine Iterate 

Figure 7,2 shows a flow chart for this subroutine. This 
consists of three loops. The inner loop, starts at level 10 and 
counts the number of members, in the structure. The middle loop 
starts at level 100 and counts the number of joints; and the out¬ 
side loop starts at level 1000 and counts the number of iterations 

NIR needed for each design. The integers I and J are used to count 
the members, and the joints as they are handled, M and N are the 
total number of members and joints respectively. 

For a particular joint J, the contributions to the modified 
stiffness equations for each member connected to the joint are cal¬ 
culated, These are given by equations (7,5), (7.6) and (7,7), 

If the deflection at joint J is specified, in one of the direction, 
X, Y or Z, then the contributions to the appropriate design equa¬ 
tion for each member connected to that joint are found instead. 

The above procedure is carried out for each joint in the space 
frame, calculating each joint deflection, and the area of the 
members A , in the first group. Several iterations will be needed 

gl 

to obtain a complete convergenceo 

7 „ 7 DESIGN EXAMPLE 1 : A SIX MEMBER ISOSTATIC SPACE FRAME 

The isostatic pin-jointed space frame shovm in Figure 7.3 was 
designed by the computer program described above. This frame was 
also designed by hand and the results were compared with those 
obtained by the program. No difference between the two results was 

observed. 

The members were allowed to have any area and were not specifi 
call)' grouped. However, for the first design cycle all the member 

areas were assumed to be equal, thus C6^ ^2 ^ 



The deflection at joint (1) in the positive y direction was specified 
to be lloll mm.. 

Using a small tolerance of 0.001, the area of the members 
obtained by the iteration technique was 402.01 mm^. This is accur¬ 
ate to two decimal places. Thirty five iterations were needed to 
achieve convergence. The result of each iteration is shown by the 
graph in Figure 7.4. The total volume of the frame was 

3 

10951853,8 mm . As a result of this design cycle the joint dis¬ 
placements were also obtained, from which the internal forces in 
the members were calculated using equation (7.15). The forces in 

the members and the design for the stress requirements are given in 

2 , 

Table 7.1. The design stresses used were 0.200 

0.250 kN/mm . These values are taken from Table 6.2.2/5.5 of the 
B/20 draft specification for a slenderness X = 50 and for a 
Robertson constant a = 5.5. Table 7.1 also shows the new ratios 
of the areas of the members. It can be seen from this that the 
area of members 2, 3 and 4 have to be increased, while the rest 
have to be reduced. 

The second design cycle was carried out with these values of 
a, and using the iteration technique. A total of only 9 iteration 
was sufficient to obtain a new set of member areas. This was equal 
to 125.95 mm^ for member 1. Multiplying this value by each a, the 
new member areas were calculated. These are given at the bottom of 
Table 7.1. The total volume for this design which satisfies both 
deflection and stress requirements was 6020203.3 mm^ This is con¬ 
siderably less than the total volume obtained for equal member areas 

in the first design cycle. 

NO further design cycles could be carried out, because the 
frame is isostatic and the member forces are independent of the 
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frame stiffness, 

7,8 DESIGN EXAMPLE 2 - M 8 MEMBER HYPERSTATIC 
SPACE FRAPIE 


To obtain a hyperstatic space frame, two members (7 and 8) 
are added to the structure designed above. Member 7 connects 
joint (5) to (1) and member 8 connects joint (6) to (1), The 
resulting frame, is shown in Figure 7,5, The dimensions, design 
stress, loading and deflection specifications are kept as in 
example 1, The frame is used here for an extensive investigation 
into the proposed design method. The members are kept ungrouped 
and each member is allowed to have a separate area. 

At the first design cycle the member areas were assumed to 
be equal to each other. For a tolerance of 0,01, fifteen iter¬ 
ations were needed to obtain the joint displacements and the member 

2 

areas. These were found to be equal to 187,61 mm , 

For the subsequent design cycles the a ratios were changed 
according to the designs for the stress requirements. The values 
of these ratios are given in Table 7,2,for various design cycles. 
This shows that and converge to zero, which means that for 
the final design, members 5 and 7 are to be removed from the 
structure. Eleven cycles were needed to obtain the final design, 
the member areas are given in Table 7,3, This also shows the 
total volume of the structure. The change in volume at each design 


cycle is represented in Figure 7,6, which shows the reduction in 
the total volume of the frame due to the distribution of the member 
areas according to the stress requirements. The member areas 
needed to satisfy the stress requirements are given in Table 7,4, 
From which, it can be seen that these are less than those needed for 
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the deflection requirements. Thus, for the final design;the member 
areas given in Table 7,3 should be used. Figure 7,7 shows the 
number of iterations per design cycle. This number reduces con¬ 
siderably after the first cycle. 

It can be noticed from the above example, that the hyperstatic 
space frame specified initially was reduced to its basic isostatic 
shape at the final design stage. To examine whether this was 
accidental or not, the frame was redesigned, but this time with 
joint 2 unloaded. All the members and the other loads remain the 
same. Again the members were allowed to have any area and no 
grouping was imposed. After 12 design cycles the final design 
was obtained, in which members 1, 2, 5, 6 and 7 were removed from 
the space frame. The remaining members had areas: 

2 2 2 
A^ = 337,5 mm ; A^ = 198,0 mm , Ag = 90,1 mm 

The total volume of the structure at the first design cycle was 
6726062,7 mm^, while at the final design this was reduced to 

'7 

2930511,0 mm , a reduction of more than 50%, 

From the final design obtained above with joint (2) unloaded, 
the space frame was also seen to be reduced to its basic isostatic 
shape, Tlie members connecting joint (2) to the supports and to 
joint (1) need not be constructed to resist the specified y^^ deflec¬ 
tion at joint (1)- Furthermore, the final design obtained for the 
deflection limitation was found to be satisfactory for the stress 
requirements as well. The areas needed to resist the internal forces 

are: 

2 2 
A = 250,0 mm“; A = 146,6 mm and Ag = 66,7 mm 
3 ^ 

These are less than those given above for the deflection limitation. 
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The same frame with joints 1 and 2 loaded'as in Figure 7„5 
was also designed, but this time grouping the members into three 
groups. The first contained members 1 and 4, member 2, 5 and 6 
were in group 2, while the third group consisted of members 3, 7 
and 8„ 

In the final design stage none of the members were removed. 

The member areas in each group, for the design which satisfied both 
deflection and stress requirements, are: 

2 2 2 
A 1 = 107.9 mm ; A ^ = 196,1 mm A _ = 216,8 mm 
gl g2 g3 

3 

with a volume of 7089646,9 mm , This is considerably greater 
than the design without grouping the members shown in Table 7,3, 

7,9 DESIGN EXAMPLE 3 ^ 

The hyperstatic pin-jointed space frames (a), (b) and (c) 
shown in Figures (7,8), (7,9) and (7,10) were also designed using 
the computer program described in Section 7,6, Tlie space frames 
(a), (b) and (c) were made hyperstatic by adding one, two, and 
three extra members respectively to their basic isostatic frame. 

The joint coordinates, the external applied loads, and the deflec¬ 
tion specifications were kept the same for the three frames. 

These are given in Table 7,5 in which joints, 3, 4, 5 and 6 are 
the supports. Each of the other joints have three degrees of free¬ 
dom. The deflection at joint (7) in the positive y direction was 
specified to be 11,11 mm, Tne members were not grouped. The 
design stresses were taken as in example 1 and 2 above. 

The member areas are given in Table 7,6 for the final designs 
of the three frames. Although, these areas are calculated to 
satisfy deflection requirements it is found that they satisfy stress 
requirements as well. In each case member 3 was removed. Member 



10 was also removed from frames (b) and (c) while member 5 was 
omitted from frame (c)» The members which have to be omitted 
from each frame are denoted by asterisks in Table 7„6„ 

At the final design stage, all the hyperstatic space frames 
specified initially were reduced to their basic isostatic shapes. 
Table 7,6 gives the volumes for the frames from which it can be 
seen that frame (b) has the smallest volume of the three. Curves 
I and II in Figure 7.11 show the change of volume of frames (a) 
and (b), Curve III in Figure 7,12 shows the same information for 
frame (c), 

7^10 DESIGN EXAMPLE 4 - A PYRAMID SHAPE SPACE FRA^'1E 

The four sided pyramid shape space frame shown in Figpre 7,13 
was designed next. Members 9 to 14 are in a plane parallel to the 
Y-Z plane. The joint coordinates and the external loads are given 
in Table 7,7, The deflection of joint (1) at the apex of the 
pyramid in +y direction was specified to be equal to 13,33 mm. 

The member areas at the final cycle are summarised in Table 7, 
The first column shows the member numbers the second gives the 
area of each member obtained to satisfy deflection requirements. 
Each member has a different area, since members were not grouped. 
Columns 3 and 4 of Table 7,8 give the internal forces in the mem¬ 
bers and the member areas obtained by designing for these forces. 
The design stresses used were 0,250 kN/mm^ in tension and 0,200 
kN/mm^ in compression. These were taken from Table 6,2,2/5,5 of 
the B/20 draft specification for A = 50 and a = 5,5, The table 
also gives the values of a at the final design cycle. 

It can be seen from Table 7,8 that the member areas needed to 
satisfy deflection requirements are larger than those needed for 



stress requirements. Thus, in the final, design the member areas 
(which are given in the second column of Table 7,8) needed to 
satisfy the deflection requirement should be used. Notice that 
members 4, 13 and 14 are removed and the final frame is statically 
determinate. 

At the first design cycle the total volume of the frame was 

3 . 3 

4987942 mm , This was reduced to 3322268 mm at the final design, 

a reduction of 50,14%, 

7,11 DESIGN EXAMPLE 5 - A SPACE FRAME WITH 781 MEMBERS 

The practical pin-jointed space frame shown in Figure 7,14 has 
781 members and 242 joints and was designed using the computer 
program described in Section 7,6, In the figure, the structure is 
cut by planes parallel to the X-Y and the Y-Z planes to show how 
the members are connected. Figure 7,15 shows an XZ section at 
level 0 (i„e, Y = 0) , The small circles in the figure represent 
the supports. There are 40 supports and each one is prevented 
from movement in the X, Y and Z directions. The other joints 

have three degrees of freedom each. The distance between two neigh 
bouring joints in XY, ZY or XZ plane is 2000,0 mm. Figure 7,16 
shows the top net of members in plane XZ at level y = 2000.0 mm. 
Various other sections in the XY and YZ planes are shown in 
Figure 7,17, It can be seen that two inclined members are used 
at the edge sections, but only one at other sections. 

The structure is subjected to a dead load of 10 kN at each 
joint and a concentrated imposed load of 140 kN at joint A which 
is at the midspan (see figure 7,15), The vertical deflection at 
A was limited to 20000/360 = 55,55 nm, to satisfy the requirements 
Specified in the draft code B/20o 




For easy construction, the members are grouped into four 

categories as is common practice. The first group consists of 

all the inclined members, the vertical members are all in group 2 

and the horizontal members are in groups 3 and 4„ These are shown 

in figures 7„15 and 7,16 respectively and are for the members in 

the bottom and the top XZ planes. 

The frame is to be built of high yield. This has a tensile 

2 

design strength of 400 N/mm and a compressive design strength 

2 

specified in Table 6,2,2/5,5 of B/20 Draft as 304 N/mm , This 
means that the slenderness X of a member should not exceed 50, 

Figure 7,18 shows the change of volume of the frame at each ■ 
design cycle. Curve (I) shows this change when the frame was 
subject only to deflection constraint at A, Curve (11) gives the 
same information but with the frame subject to stress constraints. 

It can be seen from this figure that, each point on curve (I) 
requires larger volume than a point on curve (II)> which means 
that the deflection requirements govern the design of this frame. 

It should be stressed that, although curve (I) was obtained by 
imposing the deflection requirements only, the designs obtained 
were found satisfactory when stress requirements were checked for 

each and every membero 

Twelve design cycles were required to obtain a design in which 
the ratios a, between the member areas in each group and those in 
the first group, became stable. The values of a for each design 
cycle and the area of the members in group (i) are given in Table 

7,9, A tolerance of 0,005 was adopted. 

It can be seen from Figure 7,18 that the frame has the lowest 
volume at point C obtained in the second cycle of the design process 


The member areas and the total volume for' this design are given by 
the first row of Table 7,10, The corresponding member areas 
required to satisfy stress requirements are given in the second 
row of this table„ These are less than or equal to those required 
for deflection requirements» Thus, the design given in the first 
row is safe and serviceable. 

The same frame was designed, but this time using mild-steel. 

For a slenderness ratio \ not more than 50, Table 6.2.2/5.5 of 

B/20 Draft specifies the tensile and compressive design strengths 

2 2 

for this steel as 250 N/mm , and 200 N/mm respectively. 

Seven design cycles were needed this time to stabilise the 
values of a. Table 7,11 shows the change in these ratios for 
all the design cycles. This table also gives the area of the 
members in the first groupo 

The change in the volume of the structure at each cycle is 
shown in Figure 7.19. Curve (I) shows this change when deflection 
constraints were applied, while curve (II) shows it when stress 
requirements were applied. 

From this figure it can be seen that stress requirements govern 
the design in the case. The design with the lowest volume was 
again obtained in the second design cycle. 


g3 


2 2 
682.2 ram , = 1158.1 mm 


The area of the members at the second cycle are: 

9 2 

A = 497.6 mm ; A „ = 295.7 mm 
gl 

the total volume for this design is 

+. rt-n 1 ■\r wPTP imnosed then the member areas 
deflection requirements only were iraposcu. 


7 3 

119.2639 X 10 mm . When the 


were: 

A^^ = 411.3 mm^; A^2 "" ^ 

1 3 

total volume of 98.198 x 10 mm . 


596.1 mm 


2 


2 

A = 915.5 mm with 

g4 


These areas 


are less than those 
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given abovCo Thus, stress requirements govern the design in this 
case and the first set of areas should be used to obtain a safe and 
serviceable design. 

The number of iterations per design cycle is represented in 
Figure 7,20 for a tolerance of 0,005. 
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remaining two dir¬ 
ections, and add 
to the previous 
terms 


i’lM-JOINTED 


■ACH FRAME, FLOW CHART FOR SUBROUTINE ITERATE 




































Member areas obtained in 

the first design cycle A, = A„ = A, = A = = A, = 402,01 

-i^ ^^2 1 2 3 4 5 6 


Member forces in kN 


= 23,33 = -33,33 ] p 


p^ = 63,33 p^ = 14,14 


Design of the member for I 1 
stress requirements, 
area in mm^ 


253,3 


A 2 = 166,6 


The new ratios between 
the member areas 


= 1,0 


Final design, which 
satisfy strength and 
deflection requirements 
area of the members in 
2 

mm ___ 

Total volume obtained 
in the first cycle 


A^ = 125.954 A 2 = / 


A^ = 342,0 A = 76.3 
4 5 


10951853,8 mm' 


Total volume obtained 
in the second cycle 


6020203,3 mm' 


Percentage reduction 


81,92% 


Table 7,1: 6 MEMBERS ISOSTATIC SPACE FRAME 

Various design values 


* design stress in tension- 0.250 KN/mm ^ 

d«ign stress ,n comprcssionsO.200 KN/mm 










FIGURE 7.5 


8 MEMBER HYPERSTATIC SPACE FRAME 
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2 

A_ = 22,2 mm 

A2 = 112,6 mm^ 

Arr = 337o4 

A, = 197.9 

Total 

1 



4 

volume 

A3 == 0,0 

A, = 180,2 mm^ 

6 

A^ = 0,0 

2 

Ag - 90o5 imn 

4,487 498 
mm^ 


Table 7.3: 8 MEMBER-HYPERSTATIC SPACE FRAME, MEMBER AREAS IN THE FINAL 

DESIGN, BOTH STRESS AND DEFLECTION REQUIREMENTS ARE SATISFIED 


Member 
forces in 
kN 

= -3,30 

= -16.67 

p^ = -49,97 

P 4 = 36o65 

Total volume 

Pg = 0.00 

= 33.33 

p^ = 0,00 

Pg = 16.70 

Member 
areas in i 
mm^ 

A^ = 16.5 

A^ = 83,3 

A3 = 249,8 

A. = 146,6 

4 

3,320,001 mm^ 

A3 = 0.00 

A. = 133o3 

D 

Aj = 0.0 

Ag = 66,8 


Table 7.4: 8 MEMBER-HYPERSTATIC SPACE FRAME MEMBER FORCES, AND MEMBER 

AREAS WHEN ONLY STRESS REQUIREMENTS ARE SATISFIED 


















Member 

number 

Design of 
Frame (a)^ 
area in mm^ 

Design of 
Frame (b)^ 
area in Tran“ 

Design of 
Frame 

area in mm 

1 

394.8 

393.8 

281.3 

2 

69,5 

70,3 

140,6 

3 

0,0* 

0,0* 

0.0* 

4 

33,2 

96,7 

141,7 

5 

60,3 

59,7 

0,0* 

6 

534,8 

534,4 

478,2 

7 

274,1 

274,2 

274,2 

8 

67,8 

114,8 

114,8 

9 

203,1 

105,5 

105,5 

10 

91.1 

0,0* 

0,0* 

11 

- 

98,4 

98,4 

12 

- 

- 

82,0 


Total volume 7554375 ^3 7425500 7450113 nrai' 

in rilin'^ 


* Denote the members which has been removed from the 
final design 


TABLE 7,6: FINAL DESIGNS (IN WHICH DEFLECTION AND STRESS 
REQUIREMENTS ARE SATISFIED) FOR FRAMES (a), 
(b) AND (c) SHOWN IN FIGURES (7,8), (7,9) 

AND r7„10) 
















FIGURE 7„13 THf: PYRAMID SHAPE SPACE FRAME - MEMBER AND JOINT NUMBERS 








Joint 

number 


Coordinates in m 


Load in kN 



X 

Y 

z 

w 

X 

w 

y 

w 

Z 

1 

4,00 

o 

LO 

o 

t—1 

O 

LO 

0 

t—1 

-30,0 

30,0 

5.0 

2 

2,0 

2,25 

2,25 

-5,0 

5,0 

10,0 

3 

0 

o 

2,25 

0,75 

-5,0 

5,0 

5,0 

4 

2.0 

0,75 

2,25 

-10,0 

5,0 

5,0 

5 

2,0 

0,75 

0,75 

-5,0 

5,0 

5,0 

6 

o 

o 

o 

0,0 

3,0 

0,0 

0,0 

0,0 

7 

0,0 

o 

o 

0,0 

0,0 

0,0 

0,0 

8 

0,0 

3,0 

3,0 

0,0 

0,0 

0,0 

9 

0,0 

3,0 

0,0 

0.0 

0,0 

0,0 


Table 7.7: THE PYRAMID SHAPE SPACE FRAME-JOINT 

COORDINATES AND EXTERNAL APPLIED LOADS 




















Member 

number 

Area needed 
for deflec¬ 
tion requi¬ 
rements mm 2 

IntemaJ. 
forces in 
the members 
in kN 

Area needed 
for stress 
requirement 
only in mm^ 

Value 
of a 

1 

114.4 

-19.863 

99.3 

1.0 

2 

330.8 

“5 7 o 536 

287.7 

2.897 

3 

65.2 

00 

T—1 

o 

i-H 

56,6 

0,570 

4 

0 . 0 * 

0 . 0 * 

0 , 0 * 

0 , 0 * 

5 

24.2 

5.248 

21.0 

0,211 

6 

11.7 

-2.015 

10.1 

0.102 

7 

118,4 

-20,612 

103.1 

1.038 

8 

85.8 

18.674 

74.7 

0.752 

9 

55,6 

-9.680 

o 

00 

0.487 

10 

20.0 

4.363 

17.5 

0.176 

11 

45.8 

9.971 

39.9 

0.402 

12 

57.2 

-9.961 

00 

0.502 

13 

0 . 0 * 

0 . 0 * 

0 . 0 * 

0 . 0 * 

14 

0 . 0 * 

0 . 0 * 

0 . 0 * 

0 . 0 * 

15 

81,8 

-14.201 

71.0 

0.715 

16 

276.2 

-48.056 

240.3 

2.420 

17 

22,1 

4.768 

19.1 

0.192 

18 

108.1 

23.529 

94.1 

0,948 

Total volume 
in mm^ 

3 , 322,268 


2 , 888,220 



Table 7o8: PYRAMID SHAPE SPACE FRAME - DESIGN VALUES AT THE 


FINAL CYCLE 



























]0 spans each 2000.0 iran 


h- 

Y = 2000.0 mm 


FIGURE 7,16 


MEMBERS IN GROUP 4 IN THE X-Z PLANE AT \ 


2000.0 mm 
















Design 

cycle 

“l 

“2 

“3 

“4 

Area of the 
members in 
groups (1) 

1 

1,0 

1,0 

1,0 

1.0 

653,86 

2 

1,0 

0,62287 

1,44867 

2,22499 

411,45 

3 

1.0 

0.59418 

1,37098 

2,79649 

389,34 

4 

1,0 

0,58592 

1,30943 

3,28939 

389,96 

5 

1,0 

0,577994 

1.25918 

3,84170 

361,13 

6 

1,0 

0,57470 

1,27180 

4,11295 

360,70 

7 

1,0 

0,57253 

1.28253 

4,38785 

342,56 

8 

1,0 

0.57471 

1,28683 

4,53106 

338,89 

9 

1.0 

0,57585 

1,29028 

4,65067 

337,21 

10 

1.0 

0,57626 

1,29236 

4,75762 

333,89 

11 

1.0 

0,57652 

1,29200 

4,83382 

352A7 

12 

1,0 

0,57644 

1,29209 

4,83555 

332,46 


Table 7,9: THE 781 MEMBER SPACE-FRAME HIGH YIELD STEEL 
-MEMBER AREAS IN THE FIRST GROUP AND VALUES 
OF a FOR EACH DESIGN CYCLE 





Areas in Areas in Areas in Areas in Total 




























Design 

cycle 

“l 

“2 

“3 

“4 

Area of the 
members in 
groups Cl) 
mm^ 

1 

1,0 

1,0 

1,0 

1,0 

653.86 

2 

1,0 

0,6229 

1,4487 

2,2250 

411,45 

3 

1,0 

0.5942 

1.3710 

2,3274 

412,34 

4 

1,0 

0,5861 

1,324 

2.3925 

413,34 

5 

1,0 

0.5828 

1,2916 

2,4370 

412,77 

6 

i 

1,0 

0,5827 

1.2789 

2.4575 

412.06 

1 

/ 

1,0 

1 

] 

1 

_1 

0.5822 

1,2738 

2,4548 

411,75 


Table 7„11: VALUES OF ct FOR THE MILD STEEL SPACE FRA.ME 











127.999J 


Curve (II) designs to satisfy stress 
requirements 


119,2639 


120.2845 


116.2183 


93.198'; 


Curve (I) designs to satisfy 
deflection requirements 



/ 

99.4364 


Design cycle 


FIGURE 7.19 TOE 781 MEMBER SPACE FRAME - MILD STEEL USED. VOLUfTE 
np THF STRUCTURE AT EACH DESIGN CYCLE 






Number of iterations for each design cycle 


Tolerance 0.005 
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CHAPTER EIGHT 

OPTIMUM DESIGN OF RIGIDLY JOINTED STEEL FRAMES 

8.1 INTRODUCTION 

In Chapter 2, once a design was obtained, an extrapolation 
technique was then used to change the beam sections of rigidly 
jointed steel frames. Although this is an approximate method for 
reducing the material cost of the frame, it was found to be satis¬ 
factory from the practical point of view. In this chapter a more 
mathematical approach is presented using a classical method of 
optimisation by calculus. The constraints involved are all equal¬ 
ities and the method of Lagrange multipliers (Hadley 1970, and 
Majid, 1980) can be used to formulate the design problem. Such an 
approach may be more efficient than using linear programming because 
in the latter case each equality requires the introduction of an 
artificial variable which is often difficult to remove from the 
simplex table (Majid, 1974 and Saka, 1975). 

Using the second moment of area of the members as the design 
variables, the objective function Z for a minimum weight design may 
be expressed as: 


Z = f(I) = I L (8.1) 

where M is the total number of members, is the second moment of 
area of member u, which has length L^o If the effects of the axial 
loads in the members are neglected then the number of the stiffness 
constraints can be reduced to 2N, where N is the number of joints 
in the frame« These constraints are: 


a 

°v 


CX) 


= b V 
V 


1 , 


o o o > 


2N 


( 8 . 2 ) 







Each o£ the above equations represents a row in the stiffness 
equations JL = K X, 

To apply the Langrange multiplier technique (see Chapter 1, 
section 1„4) equations (8„2) are first converted to the form: 

= \ [b^ - (X)] V = 1, (8,3) 

where the variable is known as the Lagrange multiplier. The 
constraints are then added to the objective function Z to form 

F (Y. X) = fcn ♦ (8.4) 


F (Y, X) is the unconstrained Lagrangian function which when 
minimum, ensures that Z is also minimum. The difference between 
F(Y, X) and f(I) given by (8,1) is that F(Y, X) has more variables 
than f(I). For example F(Y, X) contains the variables I, as well 
as the stiffness constrains g^ (X) = b^, and the Lagrange multi¬ 
pliers X, Notice that the turning points of F(Y, X) also include 
those of f(I) , 

To obtain the point o£ Absolute Minimum for the new uncon¬ 
strained objective function (8,4), its partial derivatives with 
respect to each of the M+2N variables are set to zero. Thus: 


3F 

9f 

2N 
- z 

A ^ 

91 

31 

V=1 

V 91 

u 

u 


u 

9F 

3f 

2N 


9X 

3X 

" V=1 

A V 3X 

V 

V 


1 

9F 

9X 

= b - 

V 

gvC)^) 

= 0, 


V 


0 , u 


= 0 , 


V 


V 


1 , 


= 1 


1 , 


M (8,5) 


,,, 2N (8,6) 


, 2N (8,7) 


These equations are solved iteratively to obtain the unknown second 
moments of area I, the joint displacements )( and the Lagrangian 



multipliers X, Notice that equations (8.7) simply state that 
~ fact the stiffness equations of the frame. 

These were, on their ovm, used in Chapter 2, to calculate the 
joint displacements and the second moments of area of the columns. 
Thus part of the new iteration process is merely a repeat of the 
process given in Chapter 2. Some of the equations in (8.6) are 
partial derivatives of F with respect to the joint rotations at 
the ends of the columns. These give expressions containing the 
second moments of area Ig of the beams and are thus used to cal¬ 
culate Ig. The rest of equations (8.6) and all the equations (8.5) 
are used to calculate the unknown values of A. 

8.2 A DESIGN EXAMPLE 

To explain the use of equations (8.5) to (8.7) they will be 
formulated for the frame shown in Figure 8.1 and then used to 
demonstrate the steps of the iteration process. The frame was used 
in Chapter 2, section 2.4,4 to explain the derivation of the 
modified stiffness equations. The stiffness equations L = K X 
for this frame were given by equations (2.39). These are:' 






:2S5. 


P = -132 


1 


. h 

X. 

t! 

' ' joint (1)^ 


(2) 

} ^ ^ ^ 
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/ 
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/ 

1 
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/ 

vrrtr 


L 


3660 nun 


. L/2 = 4575 mm 

1 L/2 = 4575 mm 



->4,1^----- 



FIGURE 8.1 SINGLE STOREY - SINGLE BAY FRAME 

where are the rotations of joints 1, 2 and 3 respectively, 

Xt is the horizontal sway and is the vertical sway of joint 1. 

The second moment of area for the columns and the beam are 1^ and 

respectively. The modulus of elasticity E for the steel may be 
2 

taken as 207 kN/mm”, 

The problem is to minimise the objective function Z given by: 


Z = f(I) = 2I^h + I^L (8.9) 

To apply tlie* Lagrange multiplier technique the stiffness equations 
(8.8) are first converted to the form L = S_,i or equation (8.10) below 

In equations (S.IO) each Lagrange multiplier has two suffixes. 

The first denotes the joint number 1, 2 or 3, and the second 
suffix, which is either 1 or 2, denotes the type of deflection 
being considered. One specifies a x'ertical translation (e. g. )'> 
or X.,) while tvvo refers to a joint rotation. 

It should be noticed tliat, equations (8.10) are the matrix 
form of equations (8.3). To obtain the Lagrangian function F(Y, X) 

To.41 each row of 


of eouation 


equations (8.10J is added to the 













2s:? 


objective fimction given by equation ( 8 , 9 ), This results 


in; 


F(Y,A) = 2l^h + l^L + (-Q + 


24EI, 


- 7 ;— 

j z 


6 EI^ 6 EI 

^ + ^09 C -0 -^ X. 


6 EI. 


2 2 T~ 

h ^ h 


2 "2 


4EI^ 4EI 6 EI 2EI. 

+ - ^2 - ""2 - 2 >^1 ^ 1 f-P 

h L/2 (L/2)^ ^ L/2 ^ 


6EI2 24EI2 6EI2 2EI2 

2 "^2 ■" .. ,..3 ^1 ■" .. ,.,2 "^ 3 ^ ■" ^12 ■" l 72 ~ "^2 


(L/2)“ (L/2)" (L/2)' 


8 EI 2 2EI^ 6 EI 

+ - r, + - 1- r_) + A_„ (-0-^ X, 


1 3"^ 32 

L/2 ^ L/2 


h 


r“ 


6 EI 2 2 EI 2 "^^^2 

+ - + - r, + -+ - - r^) 


2-^1 ' "1 
(L/2) L/2 h 


( 8 , 11 ) 


L/2 


The above function includes the turning points of the objec¬ 
tive function Z = f(I) and once minimised we obtain the minimum 
Z, To do so, the partial derivatives with respect to each of 
the variables are set to zero. These derivatives are given in a 
general form by equations (8,5), (8.6),(8,7), Equations (8,5) means 
that the derivatives with respect to I should be equal to zero. 

Thus % 

3F/3I^ = 0,0, 3 F/ 3 I 2 = 0,0 (8,12) 

Equations ( 8 , 6 ) gives: 


3F/3X2 = 0,0, 3F/3r2 = 0,0., 3F/3y^= 0,0 
3F/3r^ = 0,0, 3F/3r2 = 0,0 


( 8 , 12 a) 






Notice that equations (8o7), for which ='0.0' 9F/3X22 = 0,0, 

9F/9Xj^^ = 0,0, 9F/9X^2 = 0,0 and 9F/9X22 “ 0,0, are in fact the 
stiffness equations of the frame. These are given by equations 
8 ,8, and they are used on their own in Chapter 2 (equations 2,40, 
2,41, 2.43, 2,44 and 2,45) to calculate the joint displacements 
and the second moment of area of the columns. 

In the above expressions, X 2 is a specified constant (e,g, 

X 2 = h/300), and in the case of more than one storey frame the 
values of the sway are usually calculated from equation (2,51 ), 
section 2.9. However, a relevant Lagrange multiplier ^ 2 ^ is used, 
and in this case the derivative of F(Y,X) with respect to X 2 should 
be found (Hadley, 1970) although X 2 is a constant. 

The derivatives given by equation (8,12) and (8'^ 12a) are 
calculated below in the same order as in the computer program. 

This is the same order in which the equations are solved during 
the iteration process, 

1 ) 3F/9y^: The differentiation of the fianction F(Y,X) given by 

equation (8.11) with respect to y^^ is: 


6EI2 24EI2 ^^^2 

- + X., - + -T 

(L/2)2 (L/2)3 (L/2) 


(8,13) 


The same result could be obtained from the third column of S in 
equations (8,10), This is because, the elements of this column 
are all coefficients of and remain in the derivative of y^<> 

From equation (8,13), which is relevant to y^ is calculated. 

Hence 


6EI2 6EI2 

'^■ 22 (L/2)“ (L/2)^ 

24EI2 

(L/2)3 


(8,14) 
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2) 9F/9r^: This is obtained from equation (Soil) as: 

2 EI 2 8EI 2EI 

^22 (172-’ * ^ 2 kiT-'f * *32 


The above equation can also be obtained from the fourth column 
of matrix ^ in equations (8„10^„ 

From equations (8<,15) X ^2 which is relevant to the rotation 
r^ is calculated. Rearranging equation 8,15 gives: 


22 


12 


TJi ^ -Tfi 

8EI„ 


ITT 


(8,16) 


3) 9F/9Ij^; This is obtained by differentiating the function 

given by (8,11), thus: 


2h + X 


<24E 


21 


x„ - 


6E 
, 2 


6E 
, 2 


^3^ 


+ X 


22 


6E 

TT 


x^i-^ 


4E 




+ X 32 ( - 


6E 


2 ^2 



0,0 


(8,17) 


The same equation can be obtained by considering the elements 


which involve I, in the matrix S . 

1 —a 


The coefficients of I 2 dis¬ 


appear when differentiating with respect to 

From the above equation the value of X 22 is calculated. Thus 
rearranging for X ^2 gives: 

6E 


24F 

^21 *-73“ ^2 


_ T * 

2 2 


6E . , / 6E 4E ^ , 

2 r ) +- X22 ( 2 ^2 h 3-^ 

h h 


22 


h" 


h 
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4) Differentiating the fionction F(Y,X) or the fifth 


column of matrix S gives: 

—a ® 


6EI .6EI^ 2EI 4EI 4EI 

^21 —2“ ^ Hi -^ ^ H2 - ^ H2 ^ ^ ^ 

h CL/2) L/2 ^ L/2 


(8,19) 


From this equation the value of ^^2 is calculated, hence, rearran¬ 
ging equation (8,19): 


6EI. 


6EI, 


21 - 


Hi 


32 


(L/2) 


H2 


2EI, 


L/2 


4EI 4EI^ 

-f-TT^ ^ it/) 


( 8 , 20 ) 


5) 3F/3x 2: Differentiating the function F(Y,A) with respect to 

x~, or the first column in the matrix S : 

2 ’ —a - 


24EI. 

^21 “3“ 
h 


6EI 6EI 

^22 T2 ^32 ~T~ "" °°° 


( 8 , 21 ) 


from which is: 


6EI^ 

^ ^22 “71“ ^ H 2 


6EI. 


Hi- 


24EI, 


( 8 , 22 ) 


h^ 


6) 3F/3r2: This is obtained using equation (8,11) or from the 

elements of the second row of matrix s thus, 

—a. 


21 


6 EI. 4EI 4EI^ 


6EI, 


h 


11 


(L/2) 


2 ^ ^12 


2EI, 


L/2 


= 0,0 


(8,23) 






■■ ■;*rf 

In the design procedure during the iteration process^ thetvaiues 

^ 21 ^ ^ 22 ^ ^]2 calculated previous to the 

use of the above equationo Thus, this equation is used to cal¬ 
culate I 2 o Rearranging: 


6EI^ 4EI^ 

^21 TT~ " ^22 ~ 

h h 


4E_ , 6 E , 2E 

22 ITT " 11 • 2 ^ S 2 ITT 


80 2.1 Solution of the Derivative Equations by the 
Iteration Technique 


(8.24) 


The modified stiffness equations for the above example were 
solved in section 2.7. Table 2.1 gave the value of the joint dis¬ 
placements and after convergence was obtained. Using these 
values a new iteration solves the derivative equations given 
by (8,14), (8,16), (8,18), (8,20), (8,22) and (8,24), This is 

assuming (as in section 2.7) that X 2 /h = 10,17 mm and I 2 = 2,1345 
8 4 

X 10 mm , The initial values of A^ 2 > ^22’ ^32’ ^21 

to zero. 

The first iteration cycle starts by solving equations (8,14), 

and (8,16) respectively. These give = A ^2 "" Equations 

(8,18) is solved next, which gives A 22 ~ 2,40962 x 10 , For the 

rest of this cycle the new value of A 22 used for the calculation 

of all the other variables. Equations (8,20) and (8,22) gives 

9 

A^ = 0,0, and A = 2,20480 x 10 , The above values of A are now 
0 2 21 

used to calculate the first value for the beam section using equa¬ 
tion (8,24), This gives q = 2,92193 x 10^ mm"^. It can be seen 
that for the first cycle some of the values of A are zero, but 


26 ' 2 ^ '■ . 

these change as the iteration continues. For instance, in the second 

cycle we will find that = 2,75600 x 10^, X^^ = -6o02405 x 10^, 

X 22 = 3,28459 X 10^, X 22 = 3,41365 x 10^, X 2 ^ = 6,12889 x 10^, 

7 4 

and I 2 = 1,51731 x 10 mm , The procedure described is repeated 
until all the variables converge to stable values „ To test this 
convergence a tolerance is introduced, and a variable is said to 
have become stable if; 


New value - old value 
new value 


< The tolerance 


In the next three sections general expressions for the deri¬ 
vative equations will be obtained. 


8.3 THE DERVIATIVES OF THE LAGRANGIAN FUNCTION FOR A 

GENERAL JOINT IN A FRAME 

The derivatives given in equations (8,5) and (8,6) can be 
written in a general form for a typical joint j in a frame. To do 
this, it is necessary to multiply the contributions of a member, to 
the overall stiffness matrix, by the Lagrange variable. Equations 
(2,6) gave these contributions for a beam member connecting joint 
T at its first end to joint D at its second end. These equations 
are altered to: 


First end at T Second end at D 







or 
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(8.25b) 


In the above equations suffix T, and D refer to the first, and 
second end of the member respectively„ Each X has two suffixes„ 

The first denotes the joint number at the end of the member and 
the second suffix which is either 1 or 2 denotes the type of deflec¬ 
tion being considered„ For a beam, one specifies a vertical trans¬ 
lation while two refers to a joint rotation. 

It should be noticed that, in equations (8,25) the matrix 
^ is not symmetrical because each row of the original stiffness 
matrix is multiplied by a different X, For example, the element 
in the first row of ^ in equations (8,25) are multiplied by X^^, 
since this row is used to calculate the vertical deflection y^ 
at the first end of the beam. The elements in the second row are 
multiplied by X^ 2 > because this row is used to calculate the rota¬ 
tion r^o 

The contributions of a column member to the stiffness equations 
{L} = K {X} were given by equation (2,7), These are altered to: 


First 

end 

T 


Second 

end 

D 


first end at second end at 


joint T joint D 





(8,26 


i.e. 


L = S X 


It should be noted, that in equations (8,25) and (8,26) 


the derivatives of the load vectors {V^ 


Vp Mp}' and 


000003000 
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{fip oooooHp Mp} are zero because they are constants. For this 
reason they disappear by differentiation. 

The member and joint specifications were stated in section 
2,3, where a beam joint was distinguished from a column joint. 

In that section, it was also stated that a member is given the 
same number as the joint at its first end. If a beam joint j is 
considered, then three equations are needed to evaluate I.. X., 

y Jl 

and ^his joint. These equations are: 


3F _ n 9F - o 9F 

3Ig. " 9y. - WT 


(8,27) 


Notice that suffix j refers to a specific joint in the frame and 
equations (8,27) from part of the general equations (8.5) and 
(8,6) for the frame. For a column joint equations similar to 
(8,27) are written, thus 


9F n 3F ^ ^ 3F ^ 

— = 3^: = ^ ° 

Cl 1 1 


(8,28) 


In the next two sections, expressions are derived for Ij, 
and done for the modified stiffness equations in 

section 2,4, For a beam joint the differentiation 3F/3Ij need not 
be calculated. This is because in the design procedure the value 
of I for a beam section is specified during a design cycle and 
altered only at the end of that cycle. 


8,4 THE DERIVATIVES OF THE LAGRANGIAN FUNGTION FOR A 


BEAM JOINT 


The specifications for a joint in a beam were shown in Figure 


2,3 and given in Figure 8,2, The stiffness coefficients b^., d^., 
e. and f. for beam B., and b., d., e. and f. for beam B. connected 

in 1 11-L-L X 
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to joint j were given by equations ( 2 , 10 ) and ( 2 , 11 ) respectively. 
The same notation is used here to obtain 3F/3y^. and 3F/3r., 

a) 3F/3y^.: The derivative of Lagrangian function F(X, A.) with 
respect to the veitical deflection y^ is obtained from the first 


Beam B. joint 
--_ 1 _ 


Beam B 4 


FIGURE 8,2 .SPECIFICATIONS FOR A JOINT J IN A BEAM 


column of ^ in equation (8,25) 


The elements 






of this column are all coefficients of y^ and thus when differen¬ 
tial with respect to y.p, they give: 


X . T b . - A . ..d . - A, , b . - A, „d . =0 
ill J2 J kl j k2 J 


(8.29) 


Notice that the first end of beam B^ is connected to joint j and 
its second end is connected to joint k thus when constructing equa¬ 
tion (8.29) from the first column of ^ we have replaced T by j and 
D by k. 

Apart from beam B^, the beam B^ is also connected to joint j. 
For this beam, i replaces T and j replaced D and the elements in 
the third column of ^ in (8,25), differentiation with respect to 
j, give: 


* *i2h * * b2'*i = 


0,0 


(8,30) 




Adding the contributions of B. and Bgiven 
(8o29) and (^ 8 , oO) results in the following equation for 9F/3y.j 

3F/3y. = X h -A d. - A b - A,„d. - A.,b. 

3 Jt 3 32 j kl j k 2 j il 1 


+ A. d. + A. b. + A.„d. = 0 
i2i 3I1 321 


(8,31) 


Hence 


hi' '*h2h * hih * \2h * hih - h2'‘i 


- * b.) 


(8o32) 


Comparing this with equation (2.14) which was used to calculate 
Yy it will be noticed that equation (8.32) can in fact be obtained 
from (2ol4) by replacing each deflection in (2 <,14) by its corres¬ 
ponding A. Thus A^j^ replaces the vertical deflection y^ for joint 
k, Aj ^2 replaces the rotation r^ and so on. 

b) 3F/3rj: The derivative of the Lagrangian function F(X, A) with 
respect to the rotation r^ is obtained as follows; 
i) When differentiating with respect to r^ only the coefficients 
of r^ in the second column of ^ in equation (8.25) remain. All the 
other variables are kept constant and thus their coefficients in 
equation (8,25) disappear during the differentiation. Thus: 


-A.^d. + A.„e. + Ad. + A f- = 0.0 
jl] 32 j kl3 k 23 


(8.33) 


ii) The contributions of the beam to 3F/3rj are calculated 
from the fourth column of ^ in equations (8.25), as joint j is at 


the second end of beam thus 


- "il'^i ^i2^i ^ ^b2®i " °°° 


(8.34) 








iii) The expression for the derivative SF/Sr^ is obtained by adding 
equations (8,33) and (8,34), Thus; 


3F/3r = - X d. + A.„e. + X^d. + X, . - X.,d. 

J Jl J j2 j kl j k2 j il 1 


(8,35) 


+ X.^f. + X d. + A.„e. = 0,0 

i2 1 jl 1 j2 1 


Rearranging the above equation to calculate X^^ gives; 


X.„ = (+ ^^.,d. - id. - A,-£. + A.id. - X.^f. 

j2 jl j kl j k2 j il 1 i2 1 


-A.id.)/(e. + e.) (8,36) 

comparing the above equation with equation (2,17), it is possible 
to see that the above equation can be obtained from (2,17) by 
replacing each deflection by the relevant Lagrange multiplier. 

For example, replaces y. and X^^ replaces y^^ and so on, 

8,5 THE DERIVATIVES OF THE LAGRANGIAN FUNCTION FOR A 
COLUMN JOINT 

The specifications for a general configuration of a column 
joint j were shown in Figure 2,5 and repeated in Figure 8,3, in 
which two beams and and two columns and C^. were connected 
to the joint. Their stiffness coefficients were defined in section 
2,4,1, The derivatives 9F/3I ., 3F/3x. and 3F/3r. are obtained 

^ J J J 

here using the same notationo 

a) For 3F/9I The derivative of the Lagrangian function F(X,A) 
_ ^ J 

with respect to the second moment of area I^j of column Cj is 
obtained from equation (8,26), In this equation I^. appears in 
every element of matrix S^, The coefficients of all the other 
variables disappear when differentiating with respect to I^j, 





FIGURE 8.5 SPECIFICATIONS FOR A GENERAL CONFIGURATION OF A 
COLUMN! JOINT i 


Thus for column C. with joint i at its first end and n at its 
J 


second: 


9F/3I„. = L„. + A., b X. - A.,d^. r. - A b x 
Cj Cl jl C 3 I j jl Cjl 3 jl Cjl n 


- A.,d^,,r_ - A,,d^,TX, + A.^e^.^r, + A,^d„,,,x ^ 


jl Cjl n j2 Cjl j l2 Cjl ] j2 Lj 1 n 


A...r 

32 C 3 I n 


A b .,x. + A -d .,r. + A ,o x 
nl Cjl 3 nl C 3 I 3 nl C 3 1 n 


+ A ,d„.,r 
nl Cl 1 n 


A ^d .,x. + A + A .,d X 

n2 C 3 I 3 n2 C 3 I 3 n2 C 31 n 


" ^.C^CiTn = 


[8„37 


In the above equation the term Lj,. is obtained from differentiating 
the objective function f(I) with respect to This gives the 
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Length of column C^.. Equation (8,37) is rearranged to cal¬ 
culate X.^o Thus 
j2 

h2 = 'hj*hl'^‘‘’cjl''j*'‘cjrr'"cjl=‘n"'‘cjl’'n’ 




(8,38) 


+ ^n2’^^Cjl^j " ^Cjl^j“^Cjl\"®Cjl^n^^/ 


[-d_. 1 X.+e„.Tr. + d_.,x +£„.r 1 
'■ Cjl 3 Cjl 3 Cjl n Cji n^ 


b) 9F/9Xj: The derivative of the Lagrangian function F(X,X) 
with respect to the sway x^, of joint j is obtained as follows: 

i) The contributions of the column C. to 9F/9x. are calculated 
from the first column of the matrix given in equations (8,26), 
since joint j is at the first end of column and, n is at the 
second end,then inequation (8.26) j replaces T and n replaces D, 
thus: 


X.,b^. - X.^d.. - X ^b . - X^^d„. = 0.0 
jl Cj j2 Cj nl Cj n2 Cj 


(8.39) 


ii) Joint j is at the second end of column and replaces D 
in equation (8.26) while joint m which is at its first end replaces 
T.' This means that the contributions of to 9F/9Xj are given 
by the third column of the stiffness matrix, which gives: 


■X .b + X ^d + X.,b + X.„d — 0,0 
ml Cm m2 Cm j 1 Cm ;)2 Cm 


(8.40) 


iii) There are no contributions from the beams and to 
'^F/'ax.. This is because the axial extensions of these beams are 


ignored. 









iv) The expression for 8F/8x^ is thus obtainedib^ adding 
equations (8.39) and (8.40). Hence 


9F/3x. = . - A d - A b - A d„. - A ^b^ 

J JJ- Oj Cj nl Cj n2 Cj ml Cm 


+ A ~d^ + A.^b^ + A.^d = 0.0 

m2 Cm jlCm j2Cm 


(8.41) 


From the above equation the value of A^^ is calculated. This 
gives: 


A., = (+A.„d„. + A ,b-. + A ^d^. + A ,b^ 
jl j2 Cj nl Cj n2 Cj ml Cm 


•A „d„ - A..d„ )/(b^. + b^ ) 

m2 Cm j2 Cm ^ Cj Cm'^ 


(8.42) 


Once again, the above equation can be obtained from equation (2,20) 
by substituting, for each deflection, the relevant Lagrange 
multiplier and ignoring the load. 


c) aF/ar^: The derivative of F(X,A) with respect to the rotation 
rj, of joint j is obtained as follows; 

i) The contributions of the beams and to 3F/3r^ are already 
derived and were given by equation (8.33) and (8.34). 

ii) Column C^ contributes to aF/Sr^ by the coefficients of the 
second column of the metriX S in equation (8.26). Hence 




(8.43) 


Joint j is at the second end of column C^. Hence the fourth 

column of matrix S in equations (8.26) gives the contributions of 
—c 

this member to SF/Br^o Thus 

“^I'^Cra ^m2^Cm ^jl'^Cm ^j2®Cm 


iv) Adding equations (8.33), (8.34), (8.43) and (8.44) together 
gives the final expression for aF/ar^ as: 


(8.44) 




3F/9r. = -X.^d. . X.^a. a X^^d. a X^^^j - »„d. 


* \2*1 * ^Jl^l * -»jldy * X^^a^. a 


C8o45) 


^n7^ri “ ^Tnl^r-xn ■*■ ^ = 0„0 

nz Lj ml Cm m2 Cm jl Cm j2 Cm 
2 

Notice that d^ = -dEI^/L^., e^ = 4EI^./Lj etc, thus equation (8o45) 

contains the second moment of area of the beam B. and is used to 

J 

calculate the value of I for this beam during one step of the 

iteration process„ The use of equation (8.45) to calculate 1^^^. 

is convenient because, apart from ., it contains A(s) and I 

Bj 

for the columns and these will be calculated at stages previous 
to the use of equation (8,45). 

2 

In .equation (8,45) d^ = (- 6 E/Lj)Ij which can be written as 

2 

d. == I.d., where d., = - 6 E/L., Similarly e. = I.e.. where e . 

J 1 Jl Jl 1 1 1 Jl 11 

= 4E/Lj and f^ = '^here f^.^ = 2E/Ly Rearranging equation 

(8,45), we can calculate I . as: 

I . = (A.,d. - - A..d. - ^. 96 - + “ ^-o^rA 

Bl ^ il 1 i2 1 ji 1 j2 1 jl Cj j2 Cj 


-A ^d-. - A -f„. + A -jd - A 9fp_ - 
nl Cj n2 Cj ml Cm m2 Cm j1 Cm 


(8,46) 


^j 2 ®Cm^/'^- 


A.^d.^ 
ll ll 


A.-e.^ 
l2 ll 


A. , d ., 
kl jl 


\ 2 ^jl^ 


Equations (8,32), (8,36), (8.38), (8.42) and (8.46) provide the 
right nuiTibcr of oquations to calculatG the Lagrange multipliers 
and the optimum values for the second moments of area for the 
beamso These are given for a general joint and they are repeated 
for the total number of joints in the frame. The iterative tech¬ 
nique, used to solve the modified stiffness equation, is used here 
once again to obtain the values of A and I^o 
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8.6 OUTLINE OF THE DESIGN PROCEDURE 

The optimum design for a given frame is obtained by a two 
loop iteration,, The first loop solves the modified stiffness 
equations (2.14), (2.17), (2.26), and (2.38) given in Chapter 2, 
section 2.4. The second iteration loop solves the derivative 
equations given above. The procedure is summarised in the flow¬ 
chart shown in Fig 8.4 and a step by step approach would be as 
follows: 

STEP 1 - Select a set of lower bounds for the beam sections. This 
is decided by preventing the failure of each beam under vertical 
loads by the development of a beam mechanism. Under combined 
vertical and wind loads a beam section is not allowed to be smaller 
than its lower bound. 

STEP 2 - Specify the horizontal deflection at each storey level. 

An equation for calculating this is given in section 2.9. 

STEP 3 - Define the ratio of an internal to external column for the 
second moment of area as in equation (2.33). 

STEP 4 - Select an initial set of beam sections using one of the 
methods given in section 2.5. 

STEP 5 - Set the values of the Lagrange multipliers for each joint 
to zero. 

STEP 6 - Solve the modified stiffness equations (2.14), (2.17), 
(2.26) and (2.38) by iteration to obtain the values of the unknown 
sectional property I of the columns and the unknown joint displace¬ 
ments. The method of solving the modified stiffness equations 
was explained in sections 2.6 and 2.7. Notice that the solution 
is carried out until convergence is fully achieved. 

STEP 7 - Using these values of deflections and column sections. 


calculate the Lagrange multiplier for each beam joint from equations 
(8.32) and (8.36). 




Select lower bounds for beam sections under vertical 
loads 


I 

Specify the horizontal deflections 

I 


Define the ratios of the second moment of area of an 
internal to external column 



















STEP_8 - Calculate the Lagrange multipliers forpateolumn joint 
using equations (8,38) and (8,42), 

STEP 9 - Calculate the second moment of area Ig from equation 
(8,46) for each beam associated with its corresponding column, in 
step (8) „ 

STEP 10 - Repeat steps (8) and (9) for all the column joints until 
all the derivative equations are used up, 

STEP 11 - Repeat steps 7 to 10 until convergence is achieved. At 
this stage a new set of beam sections becomes available, 

STEP 12 - With the new beam sections repeat steps 5 to 11 until 
no reduction in any beam section is possible. No beam section is 
allowed to reduce below its lower bound imposed by vertical loading 
conditions, 

It is possible to solve the modified stiffness equations as 
well as the derivative equations in one iteration loop. In this 
case for any iteration cycle, the joint displacements are calcul¬ 
ated first, followed by the calculation of the Lagrange multipliers 
and then I for the beam section. This process is repeated until 

D 

convergence is achieved for the total number of the unknowns. 
However, the procedure of solving all the equations at one go was 
found to be inefficient, as it increases the number of unknowns to 
be calculated at one time, 

8,7 THE COMPUTER PROGRAM 

A computer program was written for the optimum design of 
rigidly jointed plane frames which makes use of the design pro¬ 
cedure described in the previous section. The program consists of 
three segments. These are: 

1 - Tlie MASTER segment; this follows the optimum design flow-chart 
given in Figure 8,4, 



2 - The subroutine ITERATE: as was given in,Chapter 4, 

3 - The subroutine OPTIMUM: This calculates the Lagrange multi^ 
pliers and the optimum, values for the second moments of area of 
the beam sections, making use of the derivative equations„ The 
flow-chart of this subroutine is similar to that of ITERATEo 
Except that for the subroutine OPTIMUM the contributions of a 
member to the derivative equations are calculated at each iter¬ 
ation cycle instead of the contributions to the modified stiffness 
equations„ 

The input data to this program is exactly the same as for 
the design of steel frame given in Appendix A of this thesis, 

8,8 DESIGN EXAMPLE 1 - A FOUR STOREY SINGLE BAY FRAME 

The four storey single bay frame which was shown in Figure 
2,10a is designed here using the computer program described above. 
This frame was designed in section 2,11,1 using the linear 
extrapolation technique and the final reduced cost design was 
shown in Figure 2,18, The optimum design for this frame obtained 
by Saka (1975) using a programming type of optimisation was given 
in Figure 2,19. 

Using the classical type of optimisation as described in 
this chapter the second moment of area for each member is shown in 
Figure 8,5, To calculate the weight of the frame the length L of 
each member is multiplied by the second moment of area I, and the 
sum EIL is calculated. It was found that EIL for the design with 

17 5 

classic optimisation is 1,8940x10 mm , while for the optimum 

13 5 • 

design by Snka this v/as lo8835x10 mni with a difference of 0o55^o 
This small difference may be due to the fact that Saka’s objective 
function minimised the actual volume EAL and not SILo 





Five design cycles were needed to obtain the finali OptiTnum 
design, and the Lagrange multipliers for each joint at the final 
design are given in Table 8.1„ The reduction in the sum HIL for 
each design cycle is represented in Figure 8„6„ This shows that 
the biggest reduction occurs at the second design cycle, where 
the Lagrange multipler are first used„ Table 8„2 shows the 
number of iterations per each design cycle for both the stiffness 
and the derivative iteration techniques. 

The sura ZIL for the reduced weight design which was shown in 

13 5 

Figure 2,18 was 1,903x10 mm with an increase of 0,48% over the 
optimum. This difference is small and for this reason the economy 
by using the linear extrapolation technique may be sufficient for 
practical purposes. 

The computer time used for obtaining the classical optimum 
design was 63 seconds. The time needed for the optimum design 
by Saka was 467 seconds which is more than seven times that of the 
present method, 

8,9 DESIGN EXAMPLE 2 - A FOUR STOREY - TWO EQUAL BAY FRAME 

The four storey two bay frame shown in Figure 8,7 was designed 
using the computer program described in section 8,7, The height 
of each storey and the wind load at each storey level is similar to 
that of the four storey single bay frame designed above. The 
sway deflections used here are shown beside the frame in Figure 
8,7. These were obtained using the procedure described in 
section 2,9 with a = 1/300, 

Six design cycles were needed to obtain the optimum design 
shown in Figure 8,8, in which the second moment of area for each 
member is given. The initial I values for the beam sections were 
taken all equal to 111673 cm"^. The sum ZIL for each design cycle 






is shown in Figure 8.9. This shows that the biggest reduction 
IS obtained at the second design cycle, in which the Lagrange 
multipliers are first introduced. For the optimum design shown 
in Figure 8.5 all the beam sections are reduced to those needed 
by the beam mechanism type of collapse. No further reduction In 
the beam sections is allowed as this violates the strength require- 
ments for vertical loading. 

The same frame was designed by the method used in Chapter 2, 

in which the linear extrapolation technique was used to reduce 

the cost of the frame material. The initial I values for all 

beam sections were taken equal to 111673 cm"^, which are similar to 

those of the initial values for the optimum design. The reduced 

cost design obtained, is shown in Figure 8.10. The total EIL for 

13 5 

this design was 2.0065x10 mm which is 5,68% more than the 
optimum of 1.8955x10^ mm^. 

8,10 DESIGN EXAiMPLE 3 - FOUR STOREY TIVO UNEQUAL - BAY FRAME 

The four storey two unequal bay frame shown in Figure 8,11 was 
first designed by the optimum design computer program described in 
section 8,7, and then by the method used in Chapter 2, The reason 
for this is to provide a comparison of the two designs. The height 
of each storey and the deflection profile used here was the same 
as those used for the two equal-bay frames dealt with above. Both 
designs (the optimum and of Chapter 2) start by assuming the same 

4 

initial set of beam sections, and a high I value of 111673 cm , 

The optimum design for the frame is shown in Figure 8,12 for 
13 5 

which ZIL = 1,53012x10 mm , Five design cycles were needed to 

13 5 

reduce ZIL from its initial value of 7,38128x10 mm to its 


optimum value. Some of the beam sections are reduced to those needed 
by the beam mechanism type of collapse, and marked by asterisks in 





Figure 8„12„ The reduction in 2IL per design cycle is shown in 
Figure 8.13. 

The design by the method proposed in Chapter 2 is shown in 

13 5 

Figure 8.14, for which EIL = 1.63879x10 mm with an increase 
of 7.10%. However, it should be noticed that in this design the 
beam section at each storey is continuous over the two spans, 
which may reduce the erection cost of the frame. 





Second moment of area in cm 



STL = 1,8940x10^'^ mm^"" 


1 CURE 8.5 THE 4 STOREY SINGLE BAY LRAME 
OPTIMUM DESIGN USING GLASS IGA' 
OPTIMISATION 







joint 

number 

xlO^° 

^2^ 

xlO^ 

joint 

number 

xlO^Q 

^2 

xlO^ 

1 

7,713 

2.157 

7 

2.929 

0.885 

2 

1.399 

0.735 

8 

1.798 

2.345 

o 

4.438 

1.798 

9 

7.373 

5.534 

4 

10,717 

3,217 

10 

0,994 

1.602 

5 

0.148 

0.741 

11 

9.353 

11.121 

6 

16.671 

8.188 

12 

2.068 

1.710 


TABLE 8.1 THE 4 STOREY-SINGLE BAY FRAME LAGRANGE 
MULTIPLIERS FOR EACH JOINT IN THE FRAME 


Number of iterations 
for modified 
stiffness equations 


Number of iterations 

for derviatives 28 

equations 


Cycle 

2 

Cycle 

3 

Cycle 

4 

Cycle 

5 

9 

9 

10 

10 

16 

11 

14 

14 


TABLE 8,2 THE 4 STOREY-SINGLE BAY FRAME NUMBER OF ITERATIONS 
PER DESIGN CYCLE 
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EIL = 1.8955x1013 

FIGURE 8,8 OPTIMUM DESIGN OF THE FRAME USING CLASSICAL 
OPTIMISATION 
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CHAPTER NINE 

SUGGESTIONS FOR FURTHER WORK 

The direct design method for deflection limitations presented 
in Chapter 2 to 5 is successful in designing multistorey steel 
frames in which limitations are imposed on their sway. The frames 
designed by this method proved to satisfy the strength requirements 
of the proposed new code. In fact, the design examples of Chapter 
5 show that the frames designed to deflection requirement.remain 
elastic under factorised combined loading, when a non-linear elastic 
plastic analysis is carried out. 

It was found that many existing design methods for sway frames 
which satisfy strength requirements first prove to produce sections 
which do not reduce the storey sway to acceptable limits. This 
indicate that it is more realistic and advantageous to initiate 
the design activity by satisfying the sway requirements first. 

Formulation of the design problem using the modified stiffness 
equation enables the structure to be considered in its entirety, 
without considering each storey individually as in the existing 
design methods (Moy, 1974), (Anderson and Islam, 1979). Hence 
it was possible to obtain an exact design without the need for 
simplifying assumptions, such as the insertion of points of con- 
traflexure at midpoints of the members. The introduction of the 
axial load effect into the design method was made possible using 
stability functionso 

The iteration technique used for solving the non-linear modi¬ 
fied stiffness equations was found to be effective. Hitherto this 
technique has been used only to solve sets of linear equations, but 
it is shown that it may successfully be applied to the solution 





of non-linear equations. In all the design examples solved by this 
technique convergence was always obtained. However, experience 
showed that it is difficult to achieve convergence if the iteration 
process is started with a set of infeasible beam sections. 

In this proposed design procedure, the list of universal 
beam and column sections included only those with the most econo¬ 
mical second moment of area/cost ratios. However, if it is 
required to restrict the depth of a member, it may be necessary to 
select a section that is not among this 'economical' list. In 
this case it is necessary to alter the list accordingly, UB and 
UC are selected for the beams and the columns respectively, but 
the use of UB sections for some column members may sometimes be 
found e conomica 1 , 

The computer program is written for the design of fixed base 
frames. Although most of the multistorey frames in practice are 
designed as such, the program can be modified to include pinned 
base frames. Such modifications can be easily done by changing 
the modified stiffness equations to allow for a hinge rotation at 
the base of the frame. 

The computer time needed for the direct design of the frame 
to deflection limitations is relatively small. The time consumed 
to obtain the final design for a frame, including the time for 
checking the strength requirements, is only twice that needed for 
3 ^ sixigl 0 direct analysis of solving the stiffness equations — 

K X directly. However, for larger frames i,e, those with more 
than 24 storevs, the design time does increase and this is expected 
in such frames. 

The computer storage requirements of the design method proved 
to be so small that it is in fact possible to. use the method to 



design many smaller frames on programable pocket calculators. In 
fact, many of the frames designed in this thesis were obtained by 
desk-top computers using the basic language. This is an important 
advantage of the method and is achieved because the modified stiff¬ 
ness equations of only one joint is operated upon at any given 
stage of the iteration. 

The proposed method was also applied successfully to design 
reinforced concrete sway frames. It was shown, in Chapter 6, that 
the limitations on sway deflections govern the design of some of 
the beams and columns in a frame. The examples show that for 
taller frames, the sway deflections of a frame become more and 
more demanding. 

When designing reinforced concrete frames for deflection 
limitation, the second moment of area of a cracked rectangular 
beam was considered to be constant along the length of the member. 

The proposed method is therefore only approximate in this respect 
as the sagging part of the beam acts as a T-beam, However, it 
is difficult to select a unique second moment of area which represents 
the member accurately throughout its length. 

In Chapter 7, the proposed method was extended to design 
pin-jointed space frames also to satisfy deflection limitation 
first. For this . category of structures, it was found that with 
high yield steel the design is also governed by deflection limit¬ 
ations rather than strength constraints. In fact, the author 
carried out a comparison of prices and found that the large space 
frame in example 5 of Chapter 7, is marginally cheaper when 
manufactured out of high yield steel than mild steel. In this 
practical type frame, it was found tnat only after one cycle of 
design the weight of the frame was reduced to a minimum. Further 
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redistribution o£ the material throughout the frame and in accord¬ 
ance to strength requirements did not prove to be profitable. In 
fact, this realistic example showed that the strength requirements 
of the code played no significant part in selecting the member 
areas. These latter requirements were merely used in the second 
cycle to alter member sizes (i„e, the values of a), but in the 
resulting minimum weight design the deflection requirements proved 
to dictate the member sizes. 

In many small hyperstatic space frames it was found that the 
final design is always statically determinant. This may cause 
misleading conclusions as the structures designed were small and 
grouping members together prevents such an outcome. 

In the present computer program, for the design of pin-jointed 
space frames, the deflection was specified at one joint (e.g.. at 
mid span) and the corresponding modified stiffness equations were 
used to calculate the member areas in the first group. Those 
of other groups were calculated by proportion. It is therefore 
suggested that the program is modified to specify the deflections 
in more than one joint. In that case, proportioning the members 
using stress requirements may become unnecessary. The problem 
arises as to what deflections should be used for these other joints 
and further research is needed in this field. 

In Chapter 2 to 5 a linear extraploation technique was used 
to change the beam sections of rigidly jointed steel frames. 

Although this was found to be satisfactory from a practical point 
of view, the final design is dependent on the initial beam sections. 
It was however, found that these initial sections are difficult to 
be specified satisfactorily by a simple method. Further more, a 
aiven assumed set o£ initial sections can only be improved by the 






extrapolation method without ever obtaining the optimum design. 
This is inspite of the fact that the examples solved gave near 
optimum designs. However, the Lagrange multiplier method, applied 
in Chapter 8, overcame this problem and gives the optimum set of 
beam sections without the need of starting the procedure from any 
particular set of sections. 

In this method of optimisation, the iteration technique 
was used successfully to solve the derivative equations. The aim 
was to reduce the sum, (for all the frame members), of the product 
of the second moment of area and length. It is therefore proposed 
to use the method of Lagrange multiplier to optimise the actual 
weight of the structure. Such an effort may improve the accuracy 
of the method but can also lead to mathematical complications. 
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APPENDIX A 
DATA PREPARATION 
THE DESIGN OF STEEL FRAMES 

The data fox the frame shown in Figure A1 are given below as 
an example„ 
a - Preliminary data 

One card containing: 

(1) The highest joint number. For the frame shown in Figure 
Al, this is equal to 6. 

(2) Total number of joints in the beams (this is equal to 
2 for the frame in Figure Al), 

(3) The yield stress for the steel used in KN/mm (0,250) 

(4) The load factor for vertical loading (1,75) 

2 

(5) Modulus of elasticity in KN/mm (207) 

(6) Tolerance (0,001) 

(7) 1/a for Anderson and Islam's preliminary design (400) 

(8) 1/a allowed in the code (300) 

(9) The load factor for combined loading (1,29) 

b - Beam - joints Data (1 card per joint) 

The specification of a beam-joint was given in Section 2,3, 
The data for beam-joints 1 and 2 of the frame shown in Figure Al 
would appear as follows: 

joint (1): 1 4575o 4 3 -132,0 0,0 

joint (2): 2 4575 6 5 -132,0 0.0 

The data refers to the following: 

(1) The joint number j e,g, (1) for first joint 
The length of beam in mm (4575,) 


(2) 




296 


C3} The joint number following joint j Cwhich ds 4 for beam 1), 

(4) The joint number preceding Cwhich is 3 for joint 1) 

(5j The vertical load acting at joint j in KN C'-132oOj 
(6j The externally applied moment on joint j in KN/mm (0„0j 

c - Column-joints Data (2 cards per joint) 

The specification of a column-joint was given in section 2„3„ 

The data is punched on two cards» For example, the data for joint 
3 of the frame of Figure A1 would appear as follows: 

first card 3 1 0 0 5 10<,17 3 4 

second card 3 4575, 3660o 16,0 0,0 0,0 

The First Card contains the following information: 

(1) The joint number j (e,g, 3) 

(2) The joint number at the second end of beam (which is 
1 for joint 3) 

(3) The joint number at the first end of beam B^ (which is 0 
for joint 3, because B^ does not exist) 

(4) The joint number at the second end of column (which is 
0 for joint 3, because Cj is fixed at its second end) 

(5) The joint number at the first end of Column (which is 
5 for joint 3) 

(6) The sway of the joint in mm (10.17) 

(7) The number of the extreme left joint of the storey of j (3) 

(8) The number of the extreme right joint of the storey of j (4) 

Tlie Second Card, which follows the first one, contains: 

(1) The joint number (e„g„ 3) 

(2) The length of beam B^ in mm (4575,0) 

(3) The length of column C^ in ram (3660,0) 




(4} The ejxternal horizontal load applied at joint j in KN 

C5} The external vertical load applied at joint j in KN (0,0) 

(6) The external applied moment on joint j in KN.inm COoO) 


Notice that a joint is given a number zero if it does not 
exist, or if it is fixed to a support. For example, for joint 3 
the joint niomber at the first end of beam is zero, because 
there is no beam Also the joint number at the second end of 

column Cj is zero, because this end is fixed to the ground. The 
rest of the data for joints 4, 5 and 6 of the frame shown in Figure 
A1 is given in Table Al, 


joint 

4 

Card 1: 

4 0 1 0 6 10,17 3 4 

Card 2: 

4 0,0 3660, 16,0 0,0 0,0 

joint 

5 

Card 1: 

5 2 0 3 0 20,33 5 6 

Card 2: 

5 4575, 3660, 8,0 0.0 o,0 

joint 

6 

Card 1: 

6 0 2 4 0 23,33 5 6 

Card 2: 

6 0,0 3660,0 8,0 0,0 0,0 


TABLE Al DATA FOR COLUMN-JOINTS 4, 5 AND 6 OF THE FRAME OF 

FIGURE Al 


d - Frame Data (two cards) 

The first card contains the number of storeys and the number 
Qp bays in the frame. The second card contains the length of each 
bay in mm starting from the left bay. For the frame shown in Fig-- 
ure Al these two cards are: 2 and 1 on the first card and 9150 on 


the second. 








APPENDIX B 


DATA FOR THE CONCRETE PROGRAM 

It is obvious that the steel universal column and beams are 
not used as data for concrete programo Instead some of the sectional 
properties of the concrete sections are used. These properties dif¬ 
fer from one frame to another and should be given for each individual 
frame. The data consists of: 

(i) Beam Sections Data: (One card per section) Each card contains 
the following information, 

1) The joint number j; 

2) The width of the section b; 

3) The cover of the tensile steel d'; 

4) The cover of the compression steel G; 

5) The flange depth h^ for a T-beam section; 

6) The width of the flange for a T-beam section, 

(ii) Column Sections Data: (One card per section). Each card 
provides three numbers. These are: 

1) The joint nimber, 

2) The cover (of the tensile and compression steel), 

3) The percentage of reinforcement (lOOA^^/bh) assumed to 
resist the sway of the column. This percentage is taken 
equal to 1 at the top storey and increases gradually at 
the lower storeys. 

The width of a column is not given as data because this is 
taken to be equal to the width of the beam. 




APPENDIX C 


DATA PREPARATION 

THE DESIGN OF PIN-JOINTED SPACE FRAMES 

The data for the frame shown in Figure Cl are given below 
as an example, (this frame was shown in Figure 7„3 and designed in 
section 7,7), 
a - Preliminary Data 

One card containing: 

(1) Total number of members in the frame (this is equal to six for 
the frame in Figure Cl), 

(2) Total number of joints including supports. For the frame shown 
this is equal to six, 

(3) Total number of supports (equal to four for the frame), 

(4) Total number of joints at which deflection is specified. This 
is equal to one, 

2 

(5) Modulus of elasticity in KN/mm , (200) 

2 

(6) Design stress for members in tension in KN/mm , (0,250) 

(7) Design stress for members in compression in KN/mm'^, (0.200) 

(8) Tolerance, (0.001) 

(9) Total number of groups, (For the frame in Figure Cl this is 

equal to six, as members were not grouped in this example). 

b “ Section Data (one card per group) 

Each card contains the group number and the member areas for 
the group in ram^ s.g, for the first group the card appears as 
follows: 

1 100 

If the ratios between the member areas in the groups are to be 
taken initially as equal to one, then the area of the members in 




joint (5) 
CO.3.0) 
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the remaining groups should be initially set to 100. Thus 

2 100 

3 100 

4 100 

5 100 

6 100 

c - Member Data [one card per member) 


The data for member 1 of the frame of Fig. Cl would appear 
as follows: 

15 2 1 

the data refers to the following: 

[1) the member number e.g. [1) for first member 

[2) the joint number at the first end (which is 5 for member 1) 

(3) the joint number at the second end (which is 2 for member 1) 

(4) the number of the group in which the member belongs (e.g. 1) 
The remaininr data fro members 2, 3, 4, 5 and b of the frame 

siiown in Fieure Cl, are given in Table Cl. 






Member 2 

2 

4 

2 

2 

Member 3 

3 

4 

1 

5 

Member 4 

4 

3 

1 

4 

Member 5 

5 

1 

2 

5 

Member 6 

6 

6 

2 

6 


TABLE Cl DATA FOR MEMBERS 2, 3, 4, 5, AND 6 OF THE FRAME 
OF FIGURE Cl 

d - Joint Data (one card per joint) 

The data for joint 1 of the frame of Figure Cl would appear 
as follows: 

1 4000 0 3000 lOoO 30„0 10„0 

the data refers to the following: 

(1) the joint number (e,go 1) 

(2) the joint coordinate X in mm (4000) 

(3) the joint coordinate Y in mm (0) 

(4) the joint coordinate Z in mm (3000) 

(5) the externally applied load to the joint in the X direction in 
in KN (10.0) 

(6) The externally applied load to the joint in the Y direction in 
KN (30.0) 

(7) the externally applied load to the joint in the Z direction in 
KN (lOoO) 

The rest of the data for joints 2, 3, 4, 5 and 6 of the frame, 
shoivn in Figure Cl, are given in Table C2„ 




joint 2 

4000 3000 0 10,0 20,0 10,0 

joint 3 

0 0 3000 0,0 0,0 0,0 

joint 4 

0 3000 3000 0,0 0,0 0,0 

joint 5 

0 3000 0 0,0 0,0 0,0 

joint 6 

0 0 0 0,0 0,0 0,0 

i 


TABLE C2 DATA FOR JOINTS 2, 3, 4, 5 AND 6 OF THE FRAME 
OF FIGURE Cl 

e - Restraint Data (one card per support) 

Each card contains: 

(1) The number of the joint which is restrained in some way 

(2) Degree of freedom in X, Y and Z directions„ (one if the 
joint is constrained in this direction and 0 if it is not)„ 

For the frame shown in Figure Cl joints 3, 4, 5 and 6 are 

restrained in the X, Y and Z directions, thus the data for these 
joints are: 

3 111 

4 111 

5 111 

6 111 

f - Deflection specification (only for the joints at which deflection 

is specified - one card per each such 
joint) 

For the frame of Figure Cl there is only one card„ This would 
appear as: 

1 OoO lloll OoO 0 1 0 

The above data refers to the following: 

(1) The joint number at which deflection is specified. 





Deflections x, 7 and z for this joint (the value of the 
specified deflection in mm, or 0„0 if the deflection in the 
corresponding direction is not specified)„ 

Direction specification X, Y or Z (1 for the specified direc¬ 
tion and 0 for the rest)„ 
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